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SECTION  I 
IHTRCIDUCTICK 

The  principal  object  of  this  research  is  to  determine  the  optimum 
drive  logic  configuration  for  a specific  moving  base  simulator,  the  Large 
Amplitude  Multi-mode  Aerospace  Research  Simulator  (LAMARS),  located  at 
WPAFB.  The  problem  of  "optimizing"  simulator  motion  reproduction  is  a 
general  one  which  need  not  be  application  specific.  Therefore,  any  method 
for  determining  the  optimum  LAMARS  drive  logic  configuration  is  really  a 
general  procedure,  applicable  for  any  motion  base.  This  provides  the 
motivation  to  develop  a "general"  tool  which  can  then  be  used  for  LAMARS 
optimization. 

It  is  important  to  appreciate  the  basic  concepts  of  motion  simulation 
from  the  start.  The  pilot  in  any  vehicle  is  able  to  sense  motion  as  the 
result  of  the  angular  velocity  and  specific  force  environment  motion  pro- 
vides (visual,  aural,  thermal,  etc.,  cues  aside).  The  motion  simulator 
is  controlled  by  a computerized  drive  logic  which  commands  the  simulator 
actuators  and,  in  turn,  simulator  motion.  Put  simply,  the  drive  logic 
need  only  compute  the  appropriate  commands,  based  on  the  particular  vehicle's 
equations  of  motion,  to  cause  the  simulator  to  reproduce  the  actual  motions. 
In  this  way,  the  pilot  will  be  stimulated  by  the  same  angular  velocities 
and  specific  forces  as  in  actual  flight,  and  motion  fidelity  of  simulation 
is  perfect. 

Obviously,  the  problem  is  not  that  simple.  This  is  primarily  because 
the  simulator  motion  base  is  usually  severely  constrained  to  move  within 
some  relatively  small  volume.  The  challenge  of  drive  logic  design  is  to 
present  high  fidelity  motion  cues  within  the  simulator  constraints.  Know- 
ledge of  the  human  sensory  processes  and  apparatus  is  of  some  help  here 
since  it  is  only  necessary  to  reproduce  motions  which  can  be  sensed.  But 
the  sensory  apparatus  bandpass  is  relatively  large,  so  the  attenuation 
obtained  on  the  basis  of  this  consideration  alone  is  insufficient.  At  this 
point  the  central  question  then  concerns  the  relationship  between  recovered 
motion  and  perceived  fidelity  (Ref.  l). 
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The  princ ipal  cooqponents  of  the  drive  logic  are  the  washout  filters, 
the  limiting  scheme  and  the  coordinate  transformations.  The  washout 
filters  are  the  main  object  of  the  optimization  process  presented  here. 

A washout  filter  is  simply  a device  for  systematically  attenuating  actual 
motion.  Washout  filter  properties  can  be  linear,  nonlinear,  or  time- 
varying.  Linear  and  nonlinear  filters  are  considered  in  this  report. 

Linear  washout  filters  are  characterized  by  their  gain  and  effective 
washout  break  frequency.  These  are  the  two  parameters  which  are  used  for 
the  optimization  or  "tuning"  process. 

The  first  step  in  evaluating  the  LAMARS  drive  logic  is  to  develop  a 
set  of  typical  aircraft  maneuver  scenarios.  These  scenarios  are  the  basis 
for  generating  representative  inputs  used  when  evaluating  the  washout  fil- 
ters. An  approximate  method  was  developed  to  construct  abstracts  of  the 
maneuvers  from  strip  chart  data.  The  method  uses  kinematic  equations  to 
compute  all  required  variables.  Eight  different  scenarios  were  generated. 
The  procedure  is  described  in  Section  II. 

The  next  step  is  to  design  a computer  simulation  to  evaluate  the  IAMARS 
drive  logic.  In  addition,  this  simulation  includes  models  for  the  human 
sensory  processes  and  any  sensory  thresholds  which  are  operative.  The 
inputs  to  this  drive  logic/sensory  model  simulation  are  the  angular  velo- 
city and  specific  force  time  histories  as  developed  in  the  scenario  genera- 
tion phase.  The  outputs  include  both  actual  (in  flight)  and  simulated  (in 
simulator)  quantities  at  several  points  in  the  model.  This  modeling  process 
is  described  in  Section  III. 

Once  the  various  motion  quantities  are  available  from  the  drive  logic/ 
sensory  model,  it  is  desirable  to  have  a method  for  comparing  the  actual 
and  sensed  motion  in  order  to  evaluate  the  simulation  fidelity.  To  this 
end,  a set  of  measures  and  a resulting  criterion  were  developed  to  examine 
and  optimize  the  drive  logic.  Using  these  measures,  the  scenarios  were 
evaluated  to  choose  one  which  could  represent  most  of  the  others  in  the 
later  optimization,  fhen,  using  this  representative  scenario,  the  gains 
and  break  frequencies  of  the  washouts  were  varied  in  order  to  illuminate 
the  sensitivity  of  the  measures  to  these  parameters.  This  parameter 
variation  exposed  the  critical  axes  in  terms  of  optimization,  and  also 
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provided  a functional  relationship  between  gain  and  break  frequency  for 
the  washouts  corresponding  to  given  simulator  motion  constraints.  This 
development  and  use  of  the  measures  and  criterion  are  described  in  Sec- 
tion IV. 

At  this  point,  a piloted  simulation  is  necessary  to  solidify  the  rela- 
tionship between  the  break  frequency  of  a washout  and  the  resulting  cri- 
terion value  for  each  critical  axis.  Piloted  simulation  is  used  to  expose 
this  relationship  in  terms  of  a pilot  preference  parameter  which  is  an 
integral  element  of  the  criterion.  Also,  the  simulator  experiment  deter- 
mines the  effect  of  the  g-suit,  a sensory  aid  often  vised  in  simulating 
high-g  environments.  The  experimental  plan  and  the  results  of  the  experi- 
ment are  presented  in  Section  V. 

The  results  of  the  simulator  experiment  are  then  applied  to  optimize 
the  critical  axes,  providing  a plausible  set  of  optimum  gains  and  break 
frequencies  for  the  LAMARS  linear  drive  logic  washouts.  These  optimum 
parameters  are  for  a given  drive  logic  configuration.  In  order  to  fur- 
ther examine  the  optimization  scheme,  two  nonlinear  washout  filters  were 
included  in  the  drive  logic  simulation.  Again,  a set  of  optimum  para- 
meters is  obtained  for  each  configuration.  This  optimization  process  is 
reported  on  in  Section  VI. 

The  final  results  and  conclusions  on  the  linear  washout  scheme  optimi- 
zation, as  well  as  the  g-suit  evaluation  and  the  investigation  of  the  two 
nonlinear  washout  schemes,  are  presented  in  Section  VII.  Suggestions  for 
further  research  are  found  in  Section  VIII.  Four  appendices  are  included: 
Appendix  A provides  the  kinematic  equations  used  in  the  scenario  genera- 
tion; Appendix  B presents  the  nonlinear  "subliminal"  washout  scheme  developed 
in  this  research  program  for  possible  LAMARS  application;  Appendix  C con- 
tains some  of  the  raw  fidelity  measure  data  generated  in  this  program;  and 
Appendix  D gives  a procedure  for  comparing  and  extrapolating  results  for 
different  linear  drive  logic  configurations.  Finally,  the  documentation 
of  the  software  used  in  support  of  this  study  is  found  in  the  Software 
User's  Guide,  presented  in  a separate  report  (Ref.  2). 
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SECTION  II 


SCENARIO  GENERATION 


A.  INTRODUCTION 


In  order  to  perform  optimization  of  drive  logic,  it  is  necessary  to  have 
representative  inputs  which  will  expose  the  critical  drive  logic  parameters. 
Once  the  drive  logic  has  been  "tuned"  for  a particular  representative  input 
scenario,  it  is  assumed  to  be  optimum  not  only  for  that  input,  but  also 
for  any  scenario  which  falls  within  the  representative  scenario's  particular 
motion  envelope.  Thus,  the  ultimate  purpose  served  by  this  scenario  genera- 
tion is  to  provide  a range  of  input  scenarios  from  which  to  choose  one  "rep- 
resentative" scenario.  The  optimization  is  then  performed  using  the  chosen 
scenario  as  the  input  to  the  drive  logic.  Note  the  implicit  assumption 
that  drive  logic  optimization  for  this  "representative"  scenario  will  hold 
for  all  other  scenarios  within  the  "representative"  scenario  motion  enve- 
lope. This,  in  turn,  implies  the  following  restrictions  on  the  optimiza- 
tion process:  the  process  by  which  the  drive  logic  is  optimized  should  be 
independent  of  the  drive  logic  inputs.  This  restriction  is  not  observed 
in  the  strict  sense  by  the  optimization  procedure  stated  later  on.  This 
is  an  important  point,  and  will  be  discussed  later  in  Section  IV. B. 


Eight  scenarios  were  generated  in  this  effort;  these  included  three 
handling  qualities  tasks  (HQDT),  one  air-to-ground  weapon  delivery  task, 
and  four  air-to-air  gunnery  tasks.  Table  1 presents  a brief  description 
of  each  scenario,  the  identification  number  assigned  and  its  length  in 
seconds. 


The  process  of  generating  each  scenario  started  with  abstracting  each 
scenario  in  terms  of  three  quantities. 

p — roll  acceleration  (rad/sec^) 

q — pitch  acceleration  (rad/sec^) 


• 
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Scenario 

Working 

Scene 

Length 

Description 

Notation 

Number 

of  Run 

Roll  into  4 g constant  altitude 
turn  at  about  1 g/sec 

HQPT1 

1 

10  sec 

Roll  into  4 g constant  altitude 
turn,  unload,  and  reverse  to 
opposite  4 g constant  altitude 
turn 

HQDT2 

2 

l6  sec 

Wind-up  turn  from  a constant 

2 g turn  to  a constant  6 g 
turn  at  l/3  g/sec 

HQPT4 

3 

20  sec 

Air-to-ground  delivery  from 
lefthand  racetrack,  45  deg 
dive  angle 

A/G 

4 

120  sec 

Air-to-air,  low  closure  rate, 
low  angle  off 

AA1 

5 

64  sec 

Air-to-air,  high  closure  rate, 
high  angle  off 

AA2 

6 

17  sec 

Air-to-air,  high  closure  rate, 
low  angle  off 

AA3 

7 

10  sec 

Air-to-air,  head-on  pass 

AA4 

8 

15  sec 

The  required  results  are  the  quantities  which  serve  as  typical  inputs  to 
the  drive  logic: 

PA  - in-flight  roll  rate  (rad/sec) 

qA  - in-flight  pitch  rate  (rad/sec) 

rA  - in-flight  yaw  rate  (rad/sec) 

Up0  - longitudinal  specific  force  (g)  'j 

I At  pilot's 

Vpo  — lateral  specific  force  (g)  > location  in 

I the  aircraft 

Wpo  — normal  specific  force  (g)  / 


The  f>,  4 and  Vp  input  quantities  are  used  to  drive  a computer  program  which 
uses  kinematic  equations  of  motion  to  produce  the  pA,  qA,  rA,  Upo,  Vpo  and 
*po  inputs  to  the  drive  logic.  The  kinematic  equations  are  presented  in 
Appendix  A,  and  a complete  description  of  the  computer  program  which  imple- 
ments the  scenario  generation  process  can  be  found  in  the  User's  Guide 
(Ref.  2). 

The  input  for  the  scenario  generation  consists  of  an  array  of  discrete 
p,  q and  Vp  values  which  specify  discrete  time  intervals  and  constant  values 
for  these  variables  during  the  scenario.  While  it  might  seem  odd  to  have 
chosen  $,  4 and  Vp  as  inputs , these  variables  have  been  selected  because 
they  are  relatively  easy  to  specify  for  most  flying  tasks  of  interest,  and 
their  integrals  are  smooth  curves.  Since  the  integrals  of  the  inputs  (e.g. , 
pA  and  qA)  are  used  as  inputs  to  the  drive  logic,  it  Is  assured  that  the 
drive  logic  inputs  are  smooth  and  continuous  functions. 

Two  methods  for  abstracting  the  p,  q and  Vp  inputs  were  used.  The 
first  method,  employed  for  the  air-to-ground  and  air-to-air  scenarios, 
consists  of  developing  piecewise  constant  time  histories  for  f>,  q and 
such  that  the  resulting  airplane  variables  produced  by  the  kinematic 
equations  mimic  a strip  chart  recording  of  those  same  variables.  The 
strip  charts  were  produced  by  the  Air  Force  by  flying  each  scenario  of 
interest  in  a simulator.  Table  2 presents  a list  of  the  variables 
recorded  on  the  strip  chart  for  each  scenario  and  an  indication  of  which 
of  those  variables  were  reproduced  using  the  kinematic  equations. 

The  second  method,  used  for  generating  the  HQDT  scenarios,  employed 
a strategy  similar  to  the  one  described  above.  Since  strip  chart  record- 
ings were  not  available  for  the  HQDT  scenarios,  however,  it  was  necessary 
to  construct  the  inputs  to  the  kinematic  equations  solely  on  the  basis  of 
verbal  descriptions  of  the  scenarios.  For  these  cases,  only  $ and  Vp  were 
extracted,  while  q was  computed  to  assure  that  the  airplane  was  maintained 
in  constant  altitude  flight  (see  Appendix  A).  The  variables  reproduced 
for  the  HQDT  scenarios  are  also  shown  in  Table  2. 

The  following  figures  provide  plots  of  the  resulting  airplane  varia- 
bles. Table  3 lists  the  variables  which  are  plotted,  a description  of 
those  variables,  and  the  corresponding  label  used  on  the  plots. 
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TABLE  2.  STRIP  CHART  AND  SCENARIO  VARIABlfS 


VARIA- 

BLES 

EXIST  ON  STRIP  CHART  FOR  SCENARIOS: 

REPRODUCED  FOR  SCENARIOS: 

12  5 
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TABLE  3.  PLOTTED  SCENARIO  GENERATION  VARIABLES 


PLOT 

LABEL 


DESCRIPTION 

Body  axis  pitch  acceleration  (deg/sec2) 

Body  axis  pitch  rate  (deg/sec) 

Pitch  attitude  rate  (deg/sec) 

Pitch  attitude  (deg) 

Longitudinal  specific  force  (g) 

Body  axis  roll  acceleration  (deg/sec2) 

Body  axis  roll  rate  (deg/sec) 

Bank  angle  rate  (deg/sec) 

Bank  angle  (deg) 

Lateral  specific  force  (g) 

Inertial  speed  (ft/sec) 

Body  axis  yaw  rate  (deg/ sec) 

Heading  rate  (deg/ sec) 

Heading  (deg) 

Normal  specific  force  (g) 

Inertial  speed  rate  of  change '(ft/ sec2) 

Derivative  of  inertial  speed  rate  of 
change  (ft/sec3) 

Angle  of  attack  (deg) 

Angle  of  attack  rate  (deg/sec) 


VARIABLE 


WPOD 


VTDD 


ALFA 


ALFAD 


* input  for  scenario  generation  (except  for  HQDT  scenarios) 
■*  input  for  scenario  generation 
■*  input  for  drive  logic  simulation 


Figure  1.  (Concluded) 
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Figure  2.  Scenario  2,  HQDT2 
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Figure  4.  Scenario  4,  A/G 
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Figure  4.  (Concluded) 


Figure  5.  Scenario  5,  AA1 
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Figure  5.  (Concluded) 
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Figure  6.  (Continued) 
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The  generation  of  scenario  data  consisting  of  three  angular  velocity 
inputs  (pA,  rA)  and  three  specific  force  inputs  (upo>  Vpot  *po)  is 
necessary  in  order  to  provide  inputs  for  the  drive  logic  sinxlation.  This 
section  describes  the  model  used  in  that  simulation. 


This  model  consists  of  three  parts:  the  drive  logic  Itself,  a model 
of  the  human  vestibular  system  (assuned  here  to  be  the  primary  motion  sens- 
ing apparatus),  and  the  various  sensory  thresholds  which  are  presumed  to 
be  operative.  In  addition,  a secondary  motion  cueing  device,  the  g-suit, 
is  modeled  in  connection  with  noxmal  specific  force  perception.  Although 
the  motion  base  model  simulates  the  LAMARS  drive  logic,  the  technique 
employed  is  generally  applicable.  In  order  to  apply  this  research  to  a 
different  motion  base,  it  would  only  be  necessary  to  change  the  drive  logic 
model.  Since  the  approach  to  the  overall  model  is  a modular  one,  similar 
changes  could  be  made  to  incorporate  other  sensory  models  as  well. 


The  following  sections  discuss  the  drive  logic  model,  the  perceptual 
model  (including  vestibular  sensors  and  threshold  effects)  and  the  g-suit 
model.  The  final  section  presents  the  overall  block  diagrams  for  the  models 
for  each  of  the  six  axes.  This  section  also  includes  an  explanation  of  the 
symbols  used  in  the  remainder  of  this  report. 


B.  DRIVE  LOGIC  MODEL 
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the  hard  stops  imposed  by  the  physical  travel  limits  of  the  simulator. 

The  latter  is  accomplished  by  the  combination  of  washout  filters  and  limit- 
ing computations.  Limiting  computations  are  not  considered  in  this  research 
program. 

There  are  two  distinct  methods  for  computing  the  drive  commands  for  each 
axis:  the  crossfeed  scheme  (Ref.  3)  and  the  crosaproduct  scheme  (Ref.  4). 
These  are  illustrated  in  Fig.  9,  which  presents  simplified  linearized  block 
diagrams  for  roll- sway  axis  implementation  for  each  scheme.  The  crossfeed 
scheme  was  designed  by  Bray  for  implementation  on  the  Flight  Simulator  for 
Advanced  Aircraft  (FSAA)  motion  base.  The  crossproduct  scheme  was  designed 
by  Schmidt  and  Conrad;  it  forms  the  basis  for  the  LAMARS  drive  logic. 

The  purpose  of  these  two  different  washout  implementations  is  the  same. 

In  order  to  provide  an  accurate  set  of  motion  cues  to  the  simulator  pilot, 
it  is  necessary  to  coordinate  roll  and  sway  axis  motions  and  pitch  and 
surge  axis  motions.  The  coordination  causes  an  angular  command  in  roll  or 
pitch  to  also  induce  a sway  or  surge  acceleration  of  the  motion  base  as  well. 
This  assures  that  no  spurious  lateral  or  longitudinal  specific  force  cues 
result  from  the  gravitation  vector  when  the  simulator  cab  is  tilted  from 
vertical.  To  accomplish  this  coordination,  the  cab  angular  displacement, 
cp^,  for  example,  must  be  crosscoupled  into  the  lateral  specific  force  path. 

In  order  to  simulate  low  frequency  components  of  lateral  or  longitudinal 
specific  force,  components  of  gravitational  force  must  be  used.  These 
components  are  obtained  by  tilting  the  simulator  cab.  The  low  frequency 
components  of  lateral  specific  force,  for  example,  must  be  crosscoupled 
into  the  cab  tilt  (<pm)  path  to  accomplish  this.  These  are  the  purposes  of 
the  two  paths  labeled  "residual  tilt  path"  and  "coordinating  path"  in  the 
block  diagrams  of  Fig.  9* 

Note,  however,  that  the  implementation  of  these  coordinating  paths  is 
different  for  the  two  schemes.  In  the  crossproduct  scheme  the  residual 
tilt  and  coordinating  paths  form  a closed  loop.  There  are  no  closed  loops 
in  the  crossfeed  scheme,  which  provides  the  necessary  coordination  via 
open- loop  crossfeeds.  The  following  observations  can  thus  be  made; 

• The  closed- loop  coordination  of  the  crossproduct  scheme 
provides  an  effective  second-order  washout  filter  in 
the  angular  velocity  axis. 
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• The  croasproduct  scheme  is  designed  to  produce  per- 
fect recovery  of  specific  force:  the  transfer  function 
between  Vp0  and  Vp0p  is  always  1 . 

• The  crossfeed  scheme  provides  a means  for  introducing 
deliberate  miscoordination.  This  is  accomplished  by 
varying  the  coordinating  path  gain,  Kq>y. 

The  crossproduct  scheme  is  the  basis  for  the  LAMARS  drive  logic.  The 
washout  filter  parameters  allow  the  simulator  to  be  "tuned"  to  provide  the 
most  realistic  motion  cues  available  within  the  displacement,  velocity  and 
acceleration  constraints  of  the  motion  base.  The  hey  parameters  of  the 
linear,  second-order  washouts  used  in  the  LAMARS  scheme  are  the  gain  Kp 
and  break  frequency  c%.  Variation  of  the  gain  of  the  washout  filter  pro- 
vides a method  for  attenuating  the  input  signal,  while  variation  of  the 
break  frequency  determines  the  motion  washout  rate  (Ref.  5) • 

In  the  course  of  this  research,  two  additional  nonlinear  washout  filters 
were  investigated.  Each  of  these  adds  a nonlinear  computation  to  the  under- 
lying linear  LAMARS  drive  logic.  The  first,  the  "Parrish"  washout  filter 
(Refs.  6,  7),  was  developed  to  eliminate  the  so-called  false  cue  encountered 
in  simulating  pulse- type  angular  motions.  This  technique  uses  a cost  func- 
tion to  vary  the  gain  of  the  washout  filter  automatically.  It  provides  an 
"optimum"  gain  for  the  filter  at  each  point  in  time,  based  on  a number  of 
time-varying  parameters  (such  as  angular  velocity  and  angular  displacement). 
Figure  10a  presents  a block  diagram  of  the  Parrish  filter  as  implemented  for 
an  angular  velocity  axis.  The  investigation  of  this  technique  examined  its 
use  in  the  roll  axis  only,  since  results  from  Parrish's  manned  simulation 
(Ref.  8)  indicate  that  it  is  most  effective  in  that  axis.  Detailed  infor- 
mation on  the  Parrish  washout  filter  is  in  Refs.  6-9  • 

The  second  scheme,  the  "subliminal"  washout,  was  developed  in  this 
research  program  (Refs.  10  and  11).  It,  too,  builds  on  the  underlying 
linear  LAMARS  drive  logic.  This  scheme  employs  two  nonlinear  functions 
which  increase  the  effective  washout  rate  of  the  drive  logic  whenever  the 
input  is  below  a specified  "indifference  threshold"  level.  It  has  been 
observed  by  several  researchers  that  a level  of  angular  velocity  or  specific 
force  cues  exist,  such  that  cues  which  fall  below  this  level  are  not  sensed 
by  a pilot  under  normal  workload  conditions.  This  "indifference"  threshold 
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is  exploited  in  the  subliminal  scheme.  It  provides  justification  for 
deliberately  introducing  translational  miscoordination  to  reduce  the  sway 
displacement  required.  Figure  10b  presents  a block  diagram  for  imple- 
mentation of  this  subliminal  scheme  for  an  angular  velocity  axis.  The 
subliminal,  scheme  was  included  in  the  sway  axis  only,  in  this  investigation, 
although  it  can  be  applied  for  both  sway  and  surge  translational  axes. 

Figure  11  depicts  the  roll- sway  axes  for  the  LAMARS  drive  logic  with  both 
nonlinear  schemes  included.  Further  details  of  the  subliminal,  washout  are 
presented  in  Appendix  B. 

C.  PERCEPTUAL  MODEL 

One  set  of  outputs  from  the  simulation  is  the  drive  logic  commands  to 
the  physical  motion  base  actuators.  These  consist  of  position  commands  to 
any  or  all  of  the  simulator's  axes.  A second  set  of  outputs  consists  of 
the  motion  cues  presumed  to  be  sensed  by  the  pilot.  These  include  angular 
velocities  and  specific  forces  as  appropriate  to  the  six  motion  axes.  In 
order  to  include  the  pilot  in  the  overall  simulation  of  the  LAMARS,  it  is 
necessary  to  model  the  sensory  apparatus  which  enables  the  pilot  to  perceive 
his  orientation.  This  perceptual  submodel  is  the  subject  of  this  section. 

The  vestibular  apparatus  in  the  inner  ear  is  the  primary  sensor  of 
angular  velocities  and  specific  forces  to  be  modeled  here  (Refs.  12-17) - 
The  vestibular  system  is  comprised  of  the  semi-circular  canals  and  the  oto- 
liths. Each  inner  ear  contains  three  canals,  oriented  in  an  approximately 
orthogonal  axis  system.  These  organs  are  the  angular  velocity  sensors  and 
are  modeled,  according  to  their  physical  properties,  as  heavily  damped 
angular  accelerometers.  A schematic  drawing  of  a set  of  canals  and  the 
accompanying  mathematical  model  relating  sensed  angular  velocity  to  input 
angular  velocity  is  given  in  Fig.  12. 

The  otoliths  are  the  other  sensory  organs  in  the  vestibular  system 
(Ref.  18).  Each  inner  ear  contains  two  otolith  organs:  the  saccule  is 
located  in  an  approximately  vertical  plane,  while  the  utricle  is  located 
in  a plane  inclined  approximately  25  deg  from  the  head  upright  horizontal 
plane.  The  otoliths  are  the  principle  specific  force  sensors  in  the  vestibu- 
lar system.  The  physical  properties  of  the  otoliths  are  less  well  defined 
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Figure  1 1 . IAMARS  Drive  Logic  with  Nonlinear  Washouts 
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Figure  12.  Semi-Circular  Canal  Model 
(After  Qrmsby,  Ref.  1^) 


than  the  properties  of  the  semi-circular  canals.  The  mathematical  model 
used  here,  however,  accounts  for  a major  portion  of  the  experimentally 
validated  observations  of  otolithic  function.  A schematic  drawing  of 
the  otoliths  is  shown  in  Fig.  13.  The  mathematical  model  used  to  compute 
sensed  specific  force  from  input  specific  force  is  also  included  in  Fig.  1 3- 

Threshold  effects  are  important  in  modeling  sensory  processes.  Given 
a threshold  effect,  it  is  assumed  that  motions  whose  magnitudes  fall  below 
the  threshold  level  will  not  be  perceived  by  the  pilot.  There  are  two 
types  of  thresholds  modeled  here.  The  first  is  a sensory  threshold;  it 
arises  from  physical  limitations  of  the  organ  itself.  This  threshold  type 
is  evident  for  the  semi-circular  canals  — they  are  unable  to  sense  angular 
velocities  of  magnitude  less  than  .035  rad/sec  (Ref.  15). 

The  second  type  of  threshold  is  the  "indifference"  threshold,  mentioned 
previously  in  connection  with  the  subliminal  washout  scheme.  This  threshold 
type  is  evident  for  specific  force  sensing.  Under  normal  workload,  pilots 
are  unable  to  detect  specific  forces  of  magnitude  less  than  approximately 
0.1  g (Ref.  20).  It  is  important  to  note  the  phrase  "under  nonnal  workload'' 
in  discussions  of  indifference  thresholds.  If  a subject  is  asked  to  con- 
centrate on  determining  when  the  sensation  of  specific  force  begins,  the 
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Figure  1j.  Otolith  Model 
(After  Ormsby,  Ref.  14) 

sensory  threshold  is  found  to  be  approximately  0.005  g (e.g.,  Ref.  15)  — 
the  otoliths  are  very  sensitive  to  stimulation.  If  the  subject  is  given 
a task  to  perform  in  addition  to  indicating  the  onset  of  a specific  force, 
this  "indifference"  threshold  is  found  to  be  much  higher  than  the  sensory 
threshold  — approximately  0.1  g,  as  stated  above.  In  the  simulation  of 
actual  flight  scenarios,  it  is  assumed  that  the  indifference  threshold 
will  be  operative,  since  the  pilot's  primary  task  will  not  be  mere  detec- 
tion of  the  motion  cues,  but  rather  an  actual  flying  task.  The  indiffer- 
ence threshold  may  be  visualized  not  as  a sensory  limitation  but  as  the 
result  of  an  information  processing  allocation  decision  made  in  the  central 
nervous  system  to  weight  the  primary  tasks  associated  with  flying  the  air- 
craft more  heavily  than  the  task  associated  with  specific  force  sensing. 
Although  this  allocation  process  is  not  at  all  well  understood,  the  result- 
ing threshold  effect  is  easily  modeled. 
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The  perceptual  model  Includes  the  sensory  threshold  for  angular  velo- 
city and  the  indifference  threshold  for  specific  force  in  the  appropriate 
paths.  Figure  lU  depicts  the  model  used  to  implement  these  thresholds  as 
a mathematical  function.  Also,  the  threshold  level  for  each  code  is 
indicated. 


D.  O-SUIT  MQD8L 


A second  sensory  modality  is  included  in  the  simulation  of  normal 
specific  force  perception.  This  is  necessary  because  motion  simulator 
travel  limits  severely  constrain  reproduction  of  normal  specific  force 
cues.  Although  visual  scenes  are  often  very  compelling  in  terns  of  simu- 
lating angular  and  translational  velocities,  it  is  not  possible  to  induce 
a sensation  of  8*s  visually.  This  is  because  the  primary  g- force  sensory 


o)  Presumed  Form  of  Indifference  Threshold 
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modalities  are  vestibular,  tactile  and  proprioceptive.  One  approach  to 
this  problem  of  simulating  g-force  sensations  in  an  essentially  static 
environment  is  by  means  of  tactile  stimulation  using  a g-suit  (Ref.  19) • 

The  g-suit  (or  pressure  suit)  was  designed  originally  for  use  in  high-g 
environments,  in  order  to  reduce  the  often  debilitating  head-to-toe  hydro- 
static blood  pressure  gradients  encountered  in  high-g  aircraft  maneuvers. 
Different  use  is  made  of  the  g-suit  in  the  simulator  context.  Here,  pres- 
sure is  applied  to  the  pilot  vie  the  suit  in  order  to  provide  tactile 
stimulation  similar  to  that  encountered  in  flying  as  normal  specific  force 
(g's)  vary.  The  g-suit  is  driven  in  accordance  with  the  in-flight  normal 
specific  force  generated  by  the  aircraft.  The  objective  of  g-suit  use  is 
to  induce  a more  realistic  perception  of  the  normal  specific  force  envi- 
ronment . 


The  model  used  for  the  g-suit  pressure  response  is  a lag  with  break 
frequency  at  1.0  rad/sec.  This  model  is  driven  by  the  actual  (in-flight) 
normal  specific  force,  and  produces  an  additional  sensed  quantity  which  is 
integrated  with  the  output  from  the  vestibular  perceptual  portion  of  the 
model.  In  this  way,  it  is  possible  to  examine  the  effect  of  the  g-suit  on 
the  total  perception  of  normal  specific  force.  The  model  for  the  g-suit 
is  such  that  the  combined  output  from  the  suit  and  the  vestibular  sensors 
yields  full  recovery  of  the  normal  specific  force  cue  for  the  in-flight 
situation.  The  manner  of  integrating  vestibular  and  g-suit  induced  per- 
ception of  normal  specific  force  is  shown  in  block  diagram  form  in  the 
next  subsection. 


E.  THE  TOTAL  MODEL 


The  overall  model  used  in  this  research  in  order  to  evaluate  the  LAMARS 
drive  logic  consists  of  a simulation  of  the  drive  logic  itself,  and  a 
simulation  of  the  major  sensory  channels  known  to  be  active  in  the  simu- 
lator pilot.  For  each  axis,  the  model  combines  the  drive  logic  equations 
for  that  axis,  the  appropriate  vestibular  sensor  for  that  axis  (semi- 
circular canal  model  for  angular  velocity  axes,  otolith  model  for  specific 
force  axes)  and  the  particular  threshold  in  effect  for  that  axis  (sensory 
threshold  for  angular  velocity  axes,  indifference  threshold  for  specific 


! 


I 

I I 

i i 


/ 


force  axes).  In  addition,  a model  for  perception  of  normal  specific  force 
via  tactile  stimulation  provided  by  a g-suit  is  included  in  the  heave  axis. 
It  is  presumed  that  the  g-suit  provides  a more  compelling  normal  specific 
force  simulation  than  do  the  visual  and  motion  cues  alone.  Integration 
of  the  g-suit-induced  cue  is  based  on  this  presumption. 

The  figures  which  follow  illustrate,  in  block  diagram  form,  the  total 
linear  model  for  each  pair  of  axes.  Table  4 presents  a list  of  the  major 
variables  for  each  axis,  and  Table  5 provides  a list  of  the  parameters  to 
be  set  and  their  nominal  (F-17  configuration  from  Ref.  l)  values.  There 
are  several  notes  which  accompany  Figs.  15-17: 

1 ) LAMARS  has  no  longitudinal  travel  capability;  there- 
fore, no  coefficient  values  relating  to  x,  k or  x are 
shown. 

2)  Roll- sway  axes  show  the  linear  washout  case  only;  the 
nonlinear  cases  are  shown  in  Fig.  1 1 . 

3)  Yaw  and  heave  axes  are  uncoupled. 

4)  No  perceptual  or  indifference  thresholds  are  thought 

to  be  applicable  in  the  heave  axis,  so  none  are  included. 

5)  An  equalization  is  postulated  in  the  heave  axis  so  that 
the  total  in-flight  perception  computed  by  summing  the 
g-suit  tactile  sensation  and  the  vestibular  sensation 
is  equal  to  the  in-flight  normal  specific  force. 

Further  details  on  the  actual  software  implementation  of  the  model  is  given 
in  the  Software  User's  Guide. 
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TABLE  4. 

LIST  OF  VARIABLES  FOR  MODEL  BLOCK  DIAGRAMS 

* 

VARIABLE 

> 

NAME 

DESCRIPTION 

: ■ 

- 

f\ 

tJA 

Input  pitch  velocity  (actual) 

\ 

EX  J 

<3AF 

Recovered  pitch  velocity  (simulated) 

' 

®M 

Commanded  pitch  angle 

F 

<le 

Scaled  error  between  <ja  and  q/yp 

1 

<1AP 

Actual  pitch  velocity  after  perceptual,  processing 

| 

^AFP 

Simulated  pitch  velocity  after  perceptual  processing 

cipe 

Scaled  error  between  q/^p  and  q^pp 

1 

^AI 

Actual  pitch  velocity  with  threshold  imposed  after 
perceptual  processing 

clAFI 

Simulated  pitch  velocity  with  threshold  imposed  after 
perceptual  processing 

r 

S : 

Qle 

Scaled  error  between  qAI  and  q^pj 

^po 

Input  longitudinal  specific  force  (actual) 

UpoF 

Recovered  longitudinal  specific  force  (simulated) 

r' 

i 

Ue 

Scaled  error  between  Up0  and  UpQp 

"poP 

Actual  longitudinal  specific  force  after  perceptual 
processing 

\ / 

"poFP 

Simulated  longitudinal  specific  force  after  perceptual 
processing 

Upe 

Scaled  error  between  Up0p  and  Up0pp 

"pol 

Actual  longitudinal  specific  force  with  threshold  imposed 
after  perceptual  processing 

* 

"poFI 

Simulated  longitudinal  specific  force  with  threshold 
imposed  after  perceptual  processing 

"ie 

Scaled  error  between  Up0j  and  UpQpj 

Pa 

Input  roll  velocity  (actual) 

Paf 

Recovered  roll  velocity  (simulated) 

<Pm 

Commanded  roll  angle 

Pe 

Scaled  error  between  pA  and  p^p 

1 

Pap 

Actual  roll  velocity  after  perceptual  processing 

Pafp 

Simulated  roll  velocity  after  perceptual  processing 

1 
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L... 
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Ppe 

PAI 

PAFI 

Pie 

■frpo 

^poF 

% 

▼e 

*poP 

▼poFP 

▼pe 

▼pol 

^poFI 

*Ie 

rA 

rAF 

♦m 

re 

rAP 

rAFP 

rpe 

rAI 

rAFI 

rIe 


TABLE  4.  (CONTINUED) 

Scaled  error  between  p^p  and  p^pp 

Actual  roll  velocity  with  threshold  imposed  after 
perceptual  processing 

Simulated  roll  velocity  with  threshold  imposed  after 
perceptual  processing 

Scaled  error  between  pAI  and  Pafi 

Input  lateral  specific  force  (actual) 

Recovered  lateral  specific  force  (simulated) 

Commanded  lateral  position 

Scaled  error  between  VpC  and  Vp0p 

Actual  lateral  specific  force  after  perceptual  pro- 
cessing 

Simulated  lateral  specific  force  after  perceptual 
processing 

Scaled  error  between  VpoP  and  VpoFP 

Actual  lateral  specific  force  with  thresholds  imposed 
after  perceptual  processing 

Simulated  lateral  specific  force  with  thresholds  imposed 
after  perceptual  processing 

Scaled  error  between  VpGj  and  'frpoFx 

Input  yaw  velocity  (actual) 

Recovered  yaw  velocity  (simulated) 

Commanded  yaw  angle 

Scaled  error  between  rA  and  r^p 

Actual  yaw  velocity  after  perceptual  processing 

Simulated  yaw  velocity  after  perceptual  processing 

Scaled  error  between  r^p  and  r^pp 

Actual  yaw  velocity  with  thresholds  imposed  after  per- 
ceptual processing 

Simulated  yaw  velocity  with  thresholds  imposed  after 
perceptual  processing 

Scaled  error  between  rAj  and  r^pj 
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TABLE  4.  (CONCLUDED) 


Input  normal  specific  force  (actual) 

Recovered,  normal  specific  force  (simulated) 

Commanded  vertical  position 
Scaled  error  between  Wp0  and  WpoF 

Actual  norma]  specific  force  after  perceptual  processing 

Simulated  normal  specific  force  after  perceptual  pro- 
cessing 

Scaled  error  between  WpQp  and  Wp0pp 

Actual  normal  specific  force  after  equalization 

Actual  normal  specific  force  alter  g-suit  processing 

Actual  total  perception  of  normal  specific  force 

Simulated  normal  specific  force  after  equalization 

Scaled  error  in  total  perception  of  normal  specific 
force  without  g-suit 

Scaled  error  in  total  perception  of  normal  specific 
force  with  g-suit 

Simulated  total  perception  of  normal  specific  force 


TABLE  5* 

LIST  OF  MODEL  PARAMETERS  AND 

1 

a 

THEIR  NOMINAL  VALUES 

:■ 

Parameter 

Value 

Units 

Kq 

0.5 

f 

K2P 

1.0 

F 

0.5 

rad/sec 

“5P 

0.2 

rad/sec 

<“4p 

0.2 

rad/sec 

r 

bp 

0.7 

Kp 

0.05 

K2R 

0.2 

•j 

% 

1.0 

3 

i 

<“p 

1.57 

rad/sec 

<»1R 

0.266 

rad/sec 

1 

<a2R 

0.65 

rad/sec 

“SR 

0.65 

rad/sec 

i 

“4r 

1.0 

rad/sec 

r 

ciy/(2\/c^ ) 

0.7 

c2y 

0.0025 

bp 

0.7 

j 

Kr 

0.5 

i 

* k2z 

0.15 

1 

(Dp 

1.0 

rad/sec 

C1Z 

0.05 

rad/sec 

• 

C2Z 

0.0 

(rad/sec)2 

y 

C3Z 

1.4 

rad/sec 

c4z 

1.0 

(rad/sec)2 

C5Z 

0.3 

rad/sec 

i 

< 

%sr 

0.0 

rad/sec 

^R 

0.0 

rad/sec 

®6P 

0.0 

rad/sec 

[ 
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SECTION  IV 


MEASURES,  CRITERIA  AND  PARAMETRIC  STUDY 


A.  INTRODUCTION 

Once  scenario  data  for  typical  flying  tasks  aryl  drive  logic  and  perceptual 
models  are  in  hand,  it  is  possible  to  compute  time'  histories  for  all  in-flight 
and  simulation  variables  of  interest.  In  order  to  interpret  this  data  it  is 
necessary  to  design  a method  by  which  the  resulting  time  histories  of  the 
various  variables  (listed  in  Table  4)  can  be  evaluated.  Evaluation  is  with 
respect  to  effect  upon  simulator  motion  fidelity  and  required  acceleration, 
velocity  and  displacement  as  drive  logic  parameters  are  varied.  The  method 
used  in  this  research  employs  a set  of  "measures"  computed  from  the  drive 
logic  simulation  variables.  These  measures  are  of  two  types.  Those  of  the 
first  type  are  constituents  of  a motion  fidelity  criterion.  The  constituents 
are  a dissection  of  a normalized  mean  square  error  according  to  contributory 
sources  believed  to  be  significant.  These  measures  then  re-combine  to 
form  a criterion  for  assigning  a meaningful  "value"  to  a particular  simu- 
lator configuration  (i.e.,  given  set  of  drive  logic  parameters).  This  value 
is  an  index  of  fidelity  for  simulator  motion  reproduction.  By  using  this 
criterion,  it  is  possible  to  compare  drive  logic  configurations  and  thereby 
choose  an  "optimum"  configuration. 

The  second  type  of  measure  represents  peak  normalized  acceleration, 
velocity  and  displacement  required  of  the  simulator  hardware  as  a function 
of  key  drive  logic  parameters  for  a given  scenario.  These  measures  are 
then  used  to  determine  those  combinations  of  drive  logic  parameter  values 
which  will  not  result  in  exceedence  of  the  physical  motion  capability  of  a 
given  simulator  (e.g.,  LAMARS). 

The  next  sub- section  begins  by  motivating  the  choice  of  measures  and  the 
criterion.  A list  of  the  computed  measures  for  each  axis  is  presented. 

Then,  the  computed  peak  measures  for  the  eight  scenario  inputs  to  the  nomi- 
nal LAMARS  drive  logic  are  compared.  This  results  in  the  choice  of  a "repre- 
sentative" scenario  for  use  in  the  remainder  of  the  study.  Using  the 


representative  scenario,  a parametric  study  with  respect  to  the  key  drive 
logic  parameters  is  undertaken  (see  Table  5) • Measures  for  sets  of  syste- 
matic parameter  variations  are  computed.  Comparison  of  the  measures  exposes 
the  important  sensitivities  of  the  motion  fidelity  criteria  to  the  drive 
logic  parameters,  and  provides  direction  for  a simulator  experiment  to  deter- 
mine one  remaining  key  parameter.  This  key  parameter  is  the  marginal  utility 
of  dynamic  motion  reproduction  errors  with  respect  to  motion  scaling  errors 
as  implied  by  pilot  subjective  choice  of  optimum  motion  reproduction.  A 
value  for  this  parameter  completes  the  motion  fidelity  criterion  develop- 
ment. The  criterion  is  thus  ready  for  application  as  the  basis  for  drive 
logic  optimization. 

B.  GENERAL  EXPLANATION  OF  MEASURES  AND  CRITERIA 

Suppose  a typical  motion  variable  upon  which  a measure  is  to  be  based 
is  x.  The  actual,  in-flight  motion  variable  is  xa,  while  the  simulated 
motion  variable  is  xs.  Let  the  scale  factor  applied  to  the  xa  variable 
to  reduce  the  level  of  simulated  motion  be  K,  while  the  scale  factor  for 
other  motion  variables  which  couple  with  the  xs  response  is  K.  The  purpose 
of  the  proposed  measure  is  to  provide  a meaningful  basis  for  comparing  the 
xa  and  xB  variables  produced  by  different  drive  logic  configurations. 

Consider  one  method  of  comparing  the  simulated  motion  with  the  actual 
motion:  the  error  variable  composed  of  the  difference  between  actual  and 
simulated  motion. 


e = xa  - xs  = (l  - K)xa  + (Kxa  - xs) 

The  algebraic  expression  on  the  far  right  side  of  the  equation  separates 
this  error  into  two  distinct  parts:  the  first  term  (l  — K)xa  is  the  component 
of  the  error  arising  from  scaling,  while  the  second  term  (Kxa  — xB)  is  the 
component  of  the  error  arising  from  the  washout  dynamics. 

The  process  of  separating  the  various  components  of  the  error  can  be 
taken  one  step  further.  Consider  the  mean  of  the  error  squared,  normal- 
ized by  the  mean  of  the  actual  motion  variable  squared,  x^.  The  normalization 


assures  that  this  measure  will  be  fairly  independent  of  the  particular  maneu- 
ver. This  is  a desirable  feature  which  was  identified  earlier  in  Section  II.  A. 
The  resulting  calculation  illustrates  this  approach  to  criterion  development: 

=5  n 

The  three  terms  on  the  far  right  of  this  equation  are  the  constituent  measures. 

2 

The  first  measure,  (l  — K)  , represents  the  effect  of  motion  scaling  error. 

The  second  measure  represents  the  effect  of  dynamic  error  in  reproducing  the 
scaled  motion.  The  third  measure  corrects  for  the  combined  effects  of  scaling 
and  dynamic  error. 

This  J0  criterion  and  its  component  measures  constitute  a plausible  method 
for  application  to  uncoupled  axes  such  as  heave  or  yaw.  In  fact,  it  will  be 
seen  that  this  very  criterion  is  used  for  those  axes,  with  one  important  change 
discussed  shortly.  A problem  exists,  however,  when  attempting  to  apply  this 
J0  criterion  to  the  coupled  roll-sway  or  pitch-surge  axis  sets,  since  it  is 
not  possible  to  distinguish  between  sources  of  infidelity  arising  from  dynamics 
and  those  arising  from  coupling  between  pairs  of  axes.  In  order  to  distinguish 
between  these  separate  sources,  the  J0  equation  is  reformulated  in  an  equiva- 
lent, but  more  elaborated  form,  as  follows: 


(Kxa  - xB)~ 


J,  = (1  - K)2  + — ~ — + 2(1  — K) 

4 


Kx|  - (xaxfl)_ 


K=0 


(Kxa  - X|f  - (Kxft  - x*)-=Q 
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4-2(1  - K) 


(xaxs ) k-o  “ (*•*•) 
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This  J-|  criterion  now  includes  measures  which  isolate  the  coupling 
effects  arising  from  the  two  inputs.  This  is  accomplished  by  evaluating 
the  measures  with  one  input  "shut  off"  (when  K ■ 0),  and  then  evaluating 
the  Incremental  effect  on  the  measure  with  that  same  input  "turned  on. " 
Again,  the  first  measure  accounts  for  the  effect  of  motion-scaling  error. 
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The  second  measure  represents  the  effect  of  dynamic  error  in  reproducing 
Bcaled  motion,  in  the  absence  of  coupling  between  simulator  axes.  The  third 
measure  corrects  for  the  combined  effects  of  scaling  and  washout  dynamic 
error  in  the  absence  of  coupling  between  the  simulator  axes.  The  fourth 
measure  represents  the  incremental  effect  of  coupling  between  the  simulator 
axes  on  the  dynamic  error  in  reproducing  scaled  motion.  The  fifth  measure 
represents  the  incremental  effect  of  coupling  between  simulator  axes  on  the 
combined  effects  of  scaling  and  washout  dynamic  error. 


Two  criteria  forms  have  been  developed  for  the  purpose  of  comparing 
various  drive  logic  schemes.  The  JQ  criterion  is  useful  for  uncoupled 
axes,  while  the  «T|  criterion  permits  the  inter-axis  coupling  effects  to 
be  taken  into  account.  One  further  refinement  requires  the  measures  to 
be  normalized  by  K or  K2  while  at  the  same  time  allowing  for  a weighting 
factor  k where  k = aK.  The  following  criteria  are  those  used  for  the 
remainder  of  this  study: 
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Notice  now  that  a can  be  treated  as  a marginal  utility  parameter  which 
reflects  the  pilot’s  subjective  perception  of  the  importance  of  dynamic  washout- 
induced  error  relative  to  scaling- induced  error.  This  weighting  factor,  a, 
permits  the  effects  of  dynamic  washout  error  to  be  expressed  in  units  which 
are  equivalent,  in  the  pilot’s  subjective  judgment,  to  the  motion  scaling 
error.  As  will* be  seen  in  the  next  section,  it  is  the  purpose  of  the  simulator 
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experiment  to  determine  this  pilot-preference  weighting  factor,  a,  in 
the  criterion  for  each  axis. 

The  JQ  and  J-|  criteria  have  been  specially  constructed  to  have  the 
following  characteristics: 

• JQ  = 0 or  J-j  = 0 implies  ideal  motion  recovery  (l.e., 
in-flight  and  simulated  motions  are  identical) 

• JQ  i 1 or  J|  > 1 implies  that  the  simulated  motion  has 
ceased  to  make  a positive  contribution  to  the  overall 
simulation  fidelity. 

Thus,  once  a has  been  selected  via  the  piloted  simulation  experiments,  the 
objective  of  the  drive  logic  optimization  is  to  minimize  JQ  or  Ji  by  select- 
ing  the  appropriate  scaling  gain  and  washout  break  frequency.  The  criterion 
is  applied  separately  for  each  axis.  Optimization  with  respect  to  these 
criteria  ultimately  results  in  the  determination  of  drive  logic  gains  and 
break  frequencies  for  each  axis. 

C.  ALTERNATIVE  CRITERIA  AND  MEASURES 

Three  alternative  criteria  for  motion  fidelity  in  each  axis  are  con- 
sidered. These  are  based  upon  the  normalized  mean  of  the  squared  errors 
following  drive  logic  processing  (e.g.,  [pA~pAp]g/^  and  [VpQ 
in  Fig.  i6);  following  sensory  processing  (e.g.,  fP^p  ~ PafP^/Pap  and 
[vpop  - VpoppJ^/vpop  in  Fig.  l6);  and  following  threshold  processing  (e.g., 

[PAI  " PAFI^/PAI  8111(1  f^poi  ~ W2Af0l  in  Fig.  16).  Constituent  measures 
must  be  computed  for  each  of  the  three  alternate  criteria  for  each  axis. 

Table  6 presents  a list  of  the  measures  which  are  computed  for  each  axis 
and  their  appropriate  units.  Basically,  there  are  three  sets  of  measures 
based  on  mean  square  errors  for  each  axis:  those  computed  after  the  drive 
logic  processing,  those  computed  after  sensory  processing,  and  those  computed 
after  threshold  processing.  For  the  coupled  axes  (pitch-surge  and  roll-sway), 
there  are  six  measures  in  each  set:  the  mean  of  the  squared  actual  motion, 
the  four  bracketed  measures  which  comprise  Ji,  and  the  mean  of  the  squared 
simulated  motion.  For  the  uncoupled  (yaw  and  heave)  axes,  there  are  only 
two  bracketed  measures  which  comprise  the  JQ  criterion,  so  the  total  number 
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MEASURES  COMPUTED  BY  AXIS 


used  to  represent  KzK^t  for 


M 


/ 


of  measures  per  set  is  four.  For  the  heave  axis,  there  are  additional 
measures  (GS1,  GS2,  GS5)  which  are  similar  to  the  PNL1,  PNL3,  and  FNI5 
confutations,  but  do  not  include  the  effect  of  the  g-suit. 

D.  CHOOSING  A REPRESENTATIVE  SCENARIO 

Besides  the  measures  which  comprise  the  J0  and  J1  criteria,  peak  value 
measures  are  also  confuted.  These  measures  are  the  normalized  peak  values 
of  commanded  acceleration,  yelocity  and  displacement  (refer  to  Table  6) 
which  we  the  drive  logic  outputs.  These  measures,  when  compared  to  the 
appropriate  hardware  limits  of  the  particular  motion  base,  provide  a basis 
for  choosing  the  scaling  gain  and  washout  break  frequency  combinations 
required  in  order  to  fit  the  reproduced  motions  to  the  capabilities  of  a 
given  simulator.  These  peak  value  measures  are  also  used  to  select  a repre- 
sentative scenario. 

It  was  stated  in  the  scenario  generation  discussion  that  the  purpose  of 
generating  eight  scenarios  was  to  provide  a data  base  from  which  to  choose 
the  scenario  which  best  enveloped  the  remaining  scenarios.  This  "representa- 
tive" scenario  would  then  be  employed  throughout  the  parametric  study/drive 
logic  optimization  phases.  This  is  a valid  approach,  since  the  criteria 
and  their  constituent  measures  were  deliberately  developed  to  be  fairly 
independent  of  the  input.  Thus,  optimization  based  on  the  inputs  from  a 
"representative"  scenario  can  be  said  to  be  typical  for  all  inputs  which 
fall  within  those  "representative"  limits. 

In  order  to  choose  a representative  scenario,  each  of  the  eight  scenarios 
was  used  as  the  input  to  the  drive  logic/sensor  simulation.  The  resulting 
peak  measures  for  each  axis  are  listed  as  the  Bft.CC,  FVEL  and  PDIS  entries 
in  Tables  C-1,  C-2,  and  C-3  of  Appendix  C.  Note  that  the  drive  logic  para- 
meters used  in  each  case  are  those  listed  as  the  "nominal"  parameters  in 
Table  5,  so  the  drive  logic  used  to  select  the  representative  scenarios  is 


that  currently  used  for  IAMARS 


The  measures  of  interest  here  are  the  normalized  peak  values  of  velocity 
and  displacement  for  the  angular  axes,  and  the  peak  values  of  acceleration, 
velocity  and  displacement  for  the  translational  axes.  These  peak  values 
are  normalized  by  the  scaling  gains  having  primary  influence  on  that  peak. 
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Table  7 presents  a summary  list  of  the  appropriate  IAMARS  motion  base  limits 
for  comparison  with  the  peak  values  for  each  scenario.  In  order  to  compare 
the  scenarios,  some  of  the  peak  values  are  plotted  versus  scenario.  These 
plots  are  found  in  Pig.  18  (refer  to  Table  1 for  interpretation  of  the 
scenario  numbers).  The  equation  used  to  compute  the  normalized  LAMARS  hard- 
ware limit  is  included  on  each  plot. 

As  can  be  seen  from  Table  7 or  the  plots,  especially  the  heave  and  pitch 
axis  plots.  Scenarios  6 and  8 clearly  do  not  "fit"  on  the  simulator  when  the 
nominal  drive  logic  gains  are  used.  Also,  it  is  evident  that,  for  the  most 
part.  Scenario  4 (the  air-to-ground  weapon  delivery  task)  envelopes  the 

Jj 

requirements  for  Scenarios  1 through  5 and  7.  On  the  basis  of  this  obser- 
vation, Scenario  4 is  designated  the  "representative"  scenario  (Scenarios  6 
and  8 are  presumed  to  be  inappropriate  for  simulation  for  convenience  in  pre- 
senting subsequent  methods  and  results).  It  will  be  used  as  the  drive  logic 
input  for  the  remainder  of  this  research.  Once  the  drive  logic  is  optimized 
for  Scenario  4,  it  is  assumed  that  the  same  optimization  applies  for  Scenarios 
1 through  5 and  7 as  well. 


Having  chosen  the  representative  scenario,  it  is  now  possible  to  identify 
the  upper  bounds  on  the  scaling  gain  for  each  axis.  This  is  accomplished  by 
determining  that  scaling  gain  value  which  when  multiplied  by  the  normalized 
peak  value  for  Scenario  4 results  in  a product  which  is  equal  to  the  IAMARS 
hardware  limit.  This  determination  of  the  maximum  scaling  gain  assumes 
the  other  drive  logic  parameters  are  constant.  Table  8 presents  the  results 
of  these  computations.  These  results  may  be  regarded  as  a first  tentative 
step  toward  optimizing  the  drive  logic.  However,  it  is  usually  desirable 
to  optimize  with  respect  to  the  washout  break  frequencies  as  well  as  the 
scaling  gains.  A method  for  determining  those  combinations  of  scaling  gain 
and  washout  break  frequency  which  satisfy  the  simulator  hardware  limits  is 
discussed  in  the  next  subsection. 

E.  PARAMETRIC  STUDY 

Three  goals  must  be  achieved  in  this  drive  logic  parametric  study.  First, 
a range  of  drive  logic  parameters  must  be  established  for  the  IAMBS  simulator 
experiments  at  Vought.  In  particular,  the  range  of  drive  logic  parameters 
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TABLE  8.  COMPUTED  UPPER  BOUNDS  ON  SCALING  GAINS  (POR  NOMINAL 
WASHOUT  BREAK  FREQUENCIES  OF  TABLE  5 AND  SCENARIO  4) 
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i 

: 

!, 

I 
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Pitch* 

Kq 

£ 0.14 

Roll 

Kp 

< 0.81 

Sway 

Kp 

Ky  S.  0. 14 

Yawt 

Kr 

< 1.0 

Heave 

Kz 

k2z  <0.11 

* LAMARS  pitch  limits  have  been  reduced  to  account  for 
the  derotation  required  in  connection  with  heave 
motion. 

t The  upper  bound  of  unity  on  the  scaling  gains  is 
operative  here. 


which  can  be  investigated  (i.e.,  which  will  fit)  on  that  simulator  must  be 
defined.  ("Fitting"  the  simulator  refers  to  commanding  motions  which  do 
not  exceed  the  displacement,  velocity  and  acceleration  limits  of  the  motion 
base).  The  second  goal  is  to  identify  the  key  simulator  axes  for  experi- 
mental investigation.  The  third  goal  is  to  obtain  the  drive  logic  parametric 
effects  on  the  measures.  This  will  provide  the  data  base  required  in  the 
next  section  for  refining  the  fidelity  criteria  and  for  optimizing  the  LAMARS 
drive  logic.  The  latter  data  is  also  used  to  define  the  range  of  drive  logic 
parameters  which  result  in  commanded  motions  which  will  fit  within  the  IAMARS 
limits . 

This  parametric  study  is  not  intended  to  be  exhaustive.  Rather  it  is 
structured  to  highlight  the  most  significant  effects  of  the  drive  logic 
parameters  on  the  measures  and  criteria.  To  this  end,  certain  break  fre- 
quencies of  the  LAMARS  drive  logic  are  chosen  to  covary.  For  example,  in 
the  pitch-surge  axes  the  following  relationship  is  used  to  relate  u^,  o^p 
and  o^p: 

uq  = 1.4o>5p  = 1 .4  o^p 
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Similar  relationships  are  used  for  the  other  axes  (refer  to  Table  9).  These 
relationships  limit  the  number  of  parameters  varied  to  the  essential  minimum 
number.  Obviously,  such  relationships  are  somewhat  arbitrary  and  it  is 
certainly  possible  to  vary  each  parameter  separately.  Such  a study  is  out- 
side the  bounds  of  this  investigation,  however.  Furthermore,  it  is  possible 
to  approximate  the  effects  of  changes  in  these  secondary  parameters  from 
values  specified  in  the  given  relationships  by  defining  "effective  values" 
of  the  essential  parameters.  One  procedure  for  defining  values  for  the 
effective  lead  frequency  and  effective  delay  parameters  is  described  in 
Appendix  D. 


Table  5 lists  the  parameters  available  for  variation  in  the  existing 


drive  logic.  The  parametric  study  uses  a linear  model  of  the  IAMARS  drive 
logic  configuration . Table  9 describes  the  parameter  variation  matrix. 

Comparing  the  two  tables  provides  a means  for  understanding  which  parameters 
other  than  scaling  gains  are  thought  to  be  critical  for  each  axis.  Scaling 
gain  effects  are  treated  via  normalization  with  respect  to  the  scaling  gains 
since  the  drive  logic  model  is  linear. 

The  parameter  variations  result  in  tables  of  measures.  Tables  C-^  through 
C-l4  in  Appendix  C.  Again,  it  is  useful  to  compare  the  peak  values  of  accelera- 
tion, velocity  and  displacement  for  each  axis  with  the  hardware  constraints 
of  the  IAMARS  motion  base  fbr  that  same  axis.  Table  10  makes  this  comparison. 

These  data  are  used  to  identify  drive  logic  parameter  combinations 
which  result  in  simulated  motions  which  will  "fit"  on  the  motion  base.  As  in 
the  previous  section,  the  concept  of  "fitting"  on  the  simulator  is  based  on 
the  normalized  IAMARS  limits.  And  again,  the  upper  bound  on  the  gain  necessary 
to  make  the  motion  fit  on  the  motion  base  can  be  computed  for  each  value  of 
break  frequency,  B(  y If  we  assume  that  the  optimum  occurs  for  the  scaling 
gain  at  its  upper  bound  a unique  relationship  between  scaling  gain  and  break 
frequency  for  motion  is  defined.  Existence  of  such  a relationship  greatly 
simplifies  the  optimization  task.  It  is  then  sufficient  to  optimize  with 
respect  to  either  break  frequency  or  gain,  rather  than  with  respect  to  both. 

Scaling  gain  versus  break  frequency  plots  for  heave,  roll  and  sway  of 


IAMARS  are  shown  in  Figs.  19,  20  and  21,  respectively.  These  plots  are 
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TABLE  10 

PEAK  VALUES  FOR  PARAMETER  VARIATION 
VS.  BREAK  FREQUENCY 


LAMARS 

BF» O.O675 

0.125 

rad/aec 

0.25 

rad/aec 

0.5 

rad/aec 

1.0 

rad/aec 

2.0 

rad/aec 

UNITS 

rad/aec 

Pitch 

De rotated  Dlapl.  (deg) 

± 4.76 

Velocity  (deg/aec) 

60.0 

Acceleration  (deg/sec2) 

K gp  * OeO 

±400.0 

De rotated  Dlapl.  (deg) 

20.36 

19.41 

11.7 

5.51 

2.12 

0.97 

Velocity  (deg/aec) 

K2P  ■ 0.5 

5.58 

5.^3 

5.0 

3.30 

2.16 

1.09 

De rotated  Dlapl.  (deg) 

21 .06 

22.01 

16.0 

12.58 

10.5 

10.42 

Velocity  (deg/aec) 

K2P  * 1 .0 

5.67 

5.65 

5.73 

4.89 

3.96 

4.51 

De  rotated  Dlapl.  (deg) 

22.79 

24. 6l 

20.34 

18.26 

19.51 

20.40 

Velocity  (deg/aec) 

5.75 

5.87 

6.47 

7-07 

7.14 

8.94 

Yaw 

Derotated  Dlapl.  (deg) 

± 4.76 

13.5 

7.28 

3.97 

1.96 

0.93 

0.35 

Velocity  (deg/aec) 

50.0 

1.94 

2.18 

1.84 

1.16 

0.72 

0.37 

Acceleration  (deg/sec2) 

±200.0 

Heave 

Diaplacement  (ft) 

± 10.0 

300.0 

210.0 

91.5 

24.0 

6.0 

1 .41 

Velocity  (ft/aec) 

13.0 

42.0 

34.5 

28.5 

8.7 

4.5 

1.5 

Acceleration  (ft/aec2) 

± 96.6 

14. U 

15.0 

0.05 

8.4 

4.8 

3.9 

Roll  (Kp  - 0.5) 

Displacement  (deg) 

± 25.0 

Velocity  (deg/aec) 

60.0 

Acceleration  (deg/aec)2 

± 460 

3 0,0 

Displacement  (deg) 

42.88 

43.6 

41.63 

26.05 

11.05 

5.33 

Velocity  (deg/sec) 

kZR  " 0.5 

27.9 

28.33 

24.01 

21 .03 

16.31 

7.88 

Displacement  (deg) 

42.8 

43.8 

41.5 

24.9 

9.37 

4.6 

Velocity  (deg/aec) 

Krg  = 1.0 

27.9 

28.3 

23.7 

20.7 

14.7 

5.1 

Displacement  (deg) 

42.9 

43.7 

4i.3 

24.0 

9.0 

4.6 

Velocity  (deg/sec) 

27.9 

28.3 

23.5 

20.5 

13.4 

4.3 

i 
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TABLE  10.  (CONCLUDED) 


LAMARS 

BP- 0.0675 

0.125 

rad/sec 

0.5 

rad/sec 

1.0 

rad/sec 

LIMITS 

rad/sec 

Sway 

Displacement  (ft) 

± 10.0 

Velocity  (ft /sec) 

10.0 

Acceleration  (ft/sec2) 

± 51.52 

Kg,  * 0.0,  C2y  - 0.0025 

Displacement  (ft) 

900.0 

460.0 

195.0 

75.0 

21.5 

Velocity  (ft /sec) 

85.O 

90.0 

48.5 

27.5 

11.0 

Acceleration  (ft/sec2) 

26.9 

23.6 

24.16 

13.7 

5.84 

K gp  — 0.0,  Cgy  * 0.0225 

Displacement  (ft) 

365.0 

540.0 

180.O 

65.0 

20.0 

Velocity  (ft /sec) 

80.0 

65.O 

55.0 

30.0 

10.5 

Acceleration  (ft /sec2) 

26.4 

26.0 

21.7 

10.7 

6.33 

1C gp  * 0.0,  Cgy  * 0.25 

Displacement  (ft) 

80.0 

70.0 

60.0 

35.0 

12.0 

Velocity  (ft/sec) 

28.4 

30.2 

29.8 

12.9 

6.36 

Acceleration  (ft/sec2) 

13.3 

11.5 

10.8 

9.9 

6.35 

Kjj,  = 0.5,  Cgy  * 0.0225 

Displacement  (ft) 

365.0 

345.0 

180.0 

65.O 

22.0 

Velocity  (ft /sec) 

80.0 

65.0 

55.0 

30.5 

11.5 

Acceleration  (ft /sec2) 

25.0 

26.5 

22.0 

10.7 

6.4 

Kg,  = 1.0,  Cgy  » 0.0025 

Displacement  (ft) 

900.0 

470.0 

200.0 

85.O 

26.0 

Velocity  (ft/sec) 

85.0 

90.0 

50.0 

29.0 

12.5 

Acceleration  (ft /sec2) 

25.6 

22.6 

2U.1 

14.2 

7.4 

Kg,  * 1.0,  Cgy  * 0.0225 

Displacement  (ft) 

365.0 

345.0 

185.0 

70.0 

24.0 

Velocity  (ft /sec) 

80.0 

65.0 

55-0 

31.0 

12.5 

Acceleration  (ft/sec2) 

24.4 

26.9 

21.8 

11.6 

7.1 

Kg,  » 1.0,  Cgy  > 0.25 

Displacement  (ft) 

80.0 

70.0 

60.0 

37.5 

14.0 

Velocity  (ft /sec) 

27.1 

29.0 

24.5 

13.7 

7.2 

Acceleration  (ft/sec2) 

12.2 

10.8 

10.7 

9.7 

7.9 

2.0 

rad/sec 


5.5 

3.6 
3.14 

5.5 

3.4 
3.** 

3.  *5 
2.33 
3.12 

6.5 
4.0 

4.0 

8.0 

5.0 
5.2 

8.0 

4.9 

5.2 

5.0 

3.4 

4.7 


Computed  Stray  Scaling  vs.  Roll  Washout  Break  Frequency 


derived  from  the  PACC,  PVEL  and  PDIS  entries  in  Tables  C-9,  and  Tables  11-1 4 
and  the  IAMARS  motion  limits  listed  in  Table  7.  Notice  that  the  sway  axis 
is  also  sensitive  to  variations  in  C2y  (as  expected),  so  Fig.  21  is  actually 
a family  of  curves.  Table  11  provides  a set  of  relationships  which  define 
the  scaling  gains  for  each  axis.  These  relationships  define  the  gains  for 
the  subsequent  simulator  experiment. 

There  is  one  important  unknown  left  in  the  optimization  process  — the 
pilot  marginal  utility  parameter,  a.  The  marginal  utility  can  be  determined 
toy  a series  of  piloted  simulation  experiments.  These  are  run  for  each 
axis  using  several  combinations  of  gain  and  break  frequency  satisfying  the 
relationships  given  in  Table  11.  Pilot  opinion  ratings  and  commentary  can 
be  used  to  establish  the  best  among  the  several  combinations  and  a relative 
ranking  of  the  combinations.  An  estimate  of  the  marginal  utility  parameter  a, 
for  the  critical  axis  may  then  be  inferred  by  means  of  the  following  process. 
The  criterion  (JQ  or  J1 ) is  evaluated  as  a function  of  a for  each  gain  and 
break  frequency  combination  used  in  the  simulation  experiment.  The  value  of  a 
which  results  in  the  minimum  criterion  value  for  the  pilot-identified  "best" 
combination  and  also  results  in  the  correct  relative  ranking  according  to 
criterion  value,  is  accepted  as  the  "true"  marginal  utility  parameter  value. 
The  results  of  such  an  experiment  will  be  presented  in  the  next  section. 

Once  the  criterion  and  associated  marginal  utility  parameter  value,  a, 
for  each  axis  are  available,  this  motion  fidelity  criterion  is  ready  for 
application. 

The  next  section  presents  the  simulator  experiments  and  results,  which 
are  used  in  the  manner  indicated  above  to  discover  appropriate  values  for 
the  marginal  utility  parameter,  a,  for  each  axis. 

TABLE  1 1 . RELATIONSHIPS  WHICH  DEFINE  LAMBS  EXPERIMENT  SCALING 
GAINS  FOR  EACH  AXIS  AS  A FUNCTION  OF  WASHOUT  BREAK  FREQUENCY 


Heave 

Kz 

- Bj2/1* 

Sway 

- 2^ 

Roll 

*P 

= 0.286  Bp 

Pitch 

Kq 

= 0.286  Bq 

Yaw 

Kr 

= 0.286  Br 

1 
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SECTION  V 


SIMULATOR  EXPERIMENTS 


A.  INTRODUCTION  AND  SUMMARY 

A piloted  simulation  experimental  series  run  on  the  Vought  Corporation 
LAMBS  moving  base  simulator  was  necessary  to  guide  development  of  the  drive 
logic  optimization  procedure.  This  provides: 

• A data  base  including  pilot  ratings  and  rankings  with 
systematic  variations  in  drive  logic  parameters. 

• A basis  for  selecting  key  parameters  in  motion  fidelity 
criteria. 

• A systematic  survey  of  critical  issues  in  moving  base 
simulation  for  moderate  maneuver  amplitude  flying  tasks, 
particularly:  g-suit  effectiveness  and  roll  sway  coordi- 
nation. 

The  first  point  is  particularly  important  because  there  is  no  extant  data 
base  for  moderate  (or  large)  maneuver  amplitude  flying  tasks  wherein  sys- 
tematic variations  of  drive  logic  parameters  have  been  explored.  Such  a 
data  base  is  required  in  order  to  formulate  an  orderly  drive  logic  optimi- 
zation procedure. 

Key  experimental  results  obtained  for  three  subject  pilots  and  three 
air-to-ground  delivery  scenarios  are; 

• G-suit  use  is  essential  to  perception  of  load  factor 
for  positive  g maneuvering. 

• Load  factor  representation  via  simulator  heave  motion 
is  so  distorted  by  washout  effects  that  maximum  sub- 
jective fidelity  results  for  zero  scaling  gain  (i.e., 
no  motion)  for  this  axis. 

• Pull  simulator  sway  coordination  for  cab  roll  angle 
is  required  at  high  frequencies. 

• Small  roll  axis  scaling  gains  and  associated  small 
values  of  roll  washout  break  frequencies  result  in 
maximum  subjective  fidelity.  Furthermore,  some  small 
amount  of  roll  axis  motion  is  required  to  ameliorate 
vertigo  when  using  limited  field-of-view  visual  dis- 
plays. 
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• Combinations  of  sway  translation  and  residual  cab  tilt 
in  roll  should  be  such  that  lateral  specific  force 
response  increases  monotonically  for  a step  lateral 
specific  force  command  in  order  to  maximize  perceived 
motion  fidelity. 

• Subjective  evaluations  of  yaw  and  pitch  rate  motion 
fidelity  are  not  particularly  sensitive  to  drive  logic 
parameter  variations  over  the  range  of  practical  interest 
for  large  amplitude,  five  degree-of- freedom  motion  bases. 

Additionally,  the  data  required  to  infer  values  of  the  marginal  utility  para- 
meters (a's)  relating  scaling  and  dynamic  components  of  motion  reproduction 
error  were  obtained.  Actual  values  for  these  parameters  are  developed  in 
the  next  section  of  the  report. 

B.  LARGE  AMPLITUDE  MOTION  SYSTEM,  COCKPIT,  AND 

SIMULATED  AIRCRAFT 

The  motion  system  is  a five  degree-of-freedom  configuration  (vertical, 
lateral,  roll,  yaw,  pitch)  that  is  driven  by  hydraulically  powered  linear 
actuators.  The  simulator  cockpit  is  supported  at  the  forward  end  of  a rigid 
twenty-four  foot  beam.  Beam  angular  motion  through  an  aft  mounted  two -axis 
gimbal  provides  cockpit  linear  heave  and  lateral  accelerations.  Cockpit 
angular  motion  in  pitch,  roll  and  yaw  is  provided  by  the  three-axis  forward 
gimbal  system. 

The  hydraulic  system  is  a 500  to  3,200  psi,  360  gpm  continuous,  560  gpm 
transient  supply  containing  three  120  gpm  electro- hydraulic  servo  controlled 
variable  displacement  pumps  and  three  250  hp  motors. 

A summary  of  the  motion  system  performance  characteristics  is  presented 
below: 


Axis 

± Displacement 

± Velocity 

± Acceleration 

Roll 

15  deg 

3k  deg/ sec 

59  rad/sec2 

Yaw 

22.5  deg 

48  deg/sec 

29  rad/sec2 

Pitch 

35  deg 

77  deg/ sec 

68  rad/sec2 

Vertical 

(Heave) 

10  ft 

25  ft/sec 

3.4g  (from  1g) 

Lateral 

10  ft 

13  ft/sec 

1.5g 
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The  LAMBS  cockpit  is  a single  place  general  purpose  design.  Its  layout  and 
general  arrangement  are  representative  of  the  A-7.  The  instrument  display 
consists  of  conventional  aircraft  instruments.  Cockpit  controls  include  a 
conventional  center  stick,  rudder  pedals  and  throttle  lever.  Variable  pri- 
mary flight  control  feel  forces  are  provided  by  an  electrohydraulic  force 
servo  system  longitudinally  and  laterally  and  a mechanical  system  direction- 
ally. 

The  external  visual  scene  is  a 36  deg  by  48  deg  (60  deg  diagonal)  field- 
of-view  presentation  on  a 1000  line  monochrome  television  display  in  the 
forward  windscreen  position.  A 4 000: 1 scale  terrain  board  with  limited 
relief  provides  a VFR  visual  scene  up  to  8000  ft  altitude.  A gunnery  range 
feature  (in  relief)  is  included  on  the  terrain  board.  A computer-driven 
g-suit  feature  is  available.  A-7  aircraft  dynamics,  valid  at  all  attitudes 
and  angle  of  attack  and  including  speed  brakes  and  stability  augmentation, 
are  simulated. 

C.  SIMULATION  FLYING  TASKS 

Three  air-to-ground  delivery  scenarios  were  used;  one  curvilinear  and 
two  conventional  with  15,  30  and  4-5  deg  dive  angles,  respectively  (it  was 
not  possible  to  use  pop-up  deliveries  because  the  visual  display  cannot 
provide  a horizon  during  the  pop-up).  These  air-to-ground  flying  tasks 
were  set  up  by  Bill  Neely,  Capt.  USAF  to  provide  a maximum  face  validity 
within  simulation  environment  constraints.  The  three  scenario  plan  views 
are  shown  in  Fig.  22.  Initial  conditions  and  g levels  are  indicated.  Each 
scenario  is  terminated  with  a 4 g pull-off  and  left  climbing  turn,  in  formal 
runs,  each  scenario  was  flown  three  times  in  order  to  obtain  an  evaluation 
of  a single  motion  drive  logic  configuration.  (All  scenarios  were  not 
necessarily  flown  to  evaluate  a single  drive  logic  configuration,  however.) 
The  fixed  sight  depression  angle  for  these  flying  tasks  was  appropriate 
to  rocket  delivery. 

Differences  in  flying  technique  for  Air  Force  and  Navy  pilots  is  evi- 
dent in  conventional  delivery  scenarios.  Air  Force  pilots  tend  to  roll-in 
in  the  dive  plane,  and  then  execute  an  unloaded  roll-out  on  to  final.  The 
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T Indicates  Target 


V = 430  kt 
H = 7700  ft 
¥ = 360  deg 


a)  Curvilinear,  15  deg  Dive  Angie 


\ V = 350  kt 
f H=  7700ft 
I ¥ = 270  deg 


b)  Conventional , 30  deg  Dive  Angie 


V = 350  kt 
H = 7700  ft 
¥ = 270  deg 


c)  Conventional , 45  deg  Dive  Angle 


Figure  22.  Air-to-Ground  Delivery  Scenarios 


Ufa  ■ a*#  Aa 


Navy  pilot  tends  to  roll-in  above  the  dive  plane,  pull  down  to  the  dive 
plane,  and  then  execute  a loaded  roll-out  onto  final. 


D.  EXPERIMENTAL  PLAN 

1.  Experimental.  Objectives 

Experimental  objectives  are  to  obtain  sufficient  data  for  preliminary 
resolution  of  the  following  issues. 

• G-suit  effectiveness  in  augmenting  perceived  load  factor 
simulation  fidelity. 

• Marginal  utility  of  motion  amplitude  fidelity  relative 
to  motion  dynamic  fidelity  (the  a parameters  of  the 
previous  section). 

• Procedure  for  determining  combinations  of  drive  logic 
parameters  which  cause  simulated  motion  for  a given 
flying  task  to  fit  on  a given  simulator. 

• Scaling  and  washout  requirements  for  simulator 
roll-sway  coordination. 

2.  Motion  Base  Drive  Logic 

The  general  form  of  the  LAMBS  drive  logic  used  in  this  experimental  pro- 

| 

gram  is  shown  in  Fig.  23.  Features  added  to  the  LAMBS  drive  logic  expressly 
for  this  program  include  the 

• roll-sway  coordination  path  having  parameters  k,  £c, 
and  o\>. 

• roll  axis  residual  tilt  path  having  parameters  Kyr 
and  Ty 

• provision  for  less  than  100  percent  derotation  in  yaw 
per  unit  sway  via  the  parameter  k^  (required  because 
full  yaw  gimbal  range  is  required  to  derotate  the  full 
sway  displacement  of  the  cab). 

Under  ideal  circumstance  of  full  coordination  at  all  frequencies,  kc  s 1.0 
and  03c  = 0.  Under  ideal  circumstances  of  full  lateral  specific  force  recovery, 

Kyr  = 1.0,  Ty  = [2^BRy{nBRY  + 2^yO)yj-1  and  Ky  = 1.0. 


88 


Motion  Base  Drive  Logic 


Drive  logic  parameters  used  in  the  experimental  series  for  each  drive 
logic  configuration  are  listed  in  a later  table. 

Drive  logic  scaling  and  dynamic  response  was  verified  by  comparing  strip 
chart  response  recordings  with  responses  on  transparent  overlays  which  were 
computed  by  a completely  independent  method.  This  comparison  was  made  for 
step  inputs  applied  to  each  of  the  six  drive  logic  input  points  while  record- 
ing each  of  the  five  drive  logic  output  servo  commands.  Dynamic  response 
fidelity  of  the  motion  base  hardware  itself  had  been  verified  at  an  earlier 
time  by  Vought  Corporation  personnel. 

3.  Experiment  Series 

Three  major  series  of  experiments  were  run. 

• Fixed-base  series.  All  drive  logic  configurations  were 
tested  to  determine  that  no  flying  task  would  result 

in  encounters  with  the  motion  base  displacement,  velocity 
or  acceleration  limits.  This  series  of  runs  was  also 
used  by  the  three  pilot  subjects  to  practice  the  flying 
tasks. 

• Basic  series  with  motion.  The  "Reference"  drive  logic 
configuration  was  used  for  all  runs.  However,  runs 
were  made  g-Suit  on/off,  roll-sway  coordination  on/off 
and  heave  motion  on/off. 

• Final  series  with  motion.  Parametric  variations  of  the 
drive  logic  were  made.  Separate  subseries  were  run  for 
investigation  of  the  heave  axis,  roll-sway  coordination 
and  roll  axis  parametric  effects. 

The  "Reference"  drive  logic  configuration  parameters  are  listed  in  Tables  12 
through  l4.  The  three  subseries  of  the  Final  Series  are  described  in  detail 
by  the  next  three  paragraphs. 

Heave  Preference  Series  — The  Heave  Preference  (HP)  Series  of  drive 
logic  configurations  are  structured  to  investigate  the  tradeoff  between 
normal  specific  force  scaling  gain  and  washout  break  frequency.  Peak  heave 
displacement  required  is  held  constant  in  this  investigation.  This  constraint 
on  peak  heave  displacement  is  satisfied  by  Kg  = const)  for  a fourth 

order  washout  circuit.  Over  the  limited  range  available  for  meaningful  experi- 
mental investigation  (0.05  i Kg  £ O.25),  the  constraint  relationship  is 


TABLE  12 


sec  I sec 


Sbrz  5z 


0.707 


S — FIXED  PARAMETERS  — ALL  CONFIGURATIONS 


Kp  I ar  = o>r  Kq 


0.707 


sec  sec  I ft/sec 


5y  ^c  he 

ft/sec2 
'eg 


0.707  O.562 


de 


ft/sec 


0.890  2.0  0.890  2 
(deg/g 

28.65) 


* Fraction  of  full  yaw  derotation 


TABLE  13.  LAMBS  — VARIABLE  PARAMETERS  — HEAVE  AXIS 


DRIVE  LOGIC 
CONFIGURATION 

“z 

“>BRZ  - Ofc 
rad/sec 

Reference  Case 

0.10 

1.0 

Heave  Preference 
(with  and  without 
g-Suit) 

(Roll  and  Sway  i 
settings) 

ixes  at  Reference  Case 

Reference  Case 

0.10 

1.00 

HP  0.15 

0.15 

1.15 

HP  0.20 

0.20 

1.30 

HP  0.25 

0.25 

1.45 

further  approximated  as  tug  = 3(Kg  — O.l)  + 1.0  rad/sec.  The  drive  logic 
configuration  parameters  used  are  listed  in  Table  13. 

Coordination  Preference  Series  — The  Coordination  Preference  (CP) 

Series  of  drive  logic  configurations  are  structured  to  investigate  the 
tradeoff  between  high  frequency  coordination  scaling  gain  and  coordination 
washout  break  frequency.  Peak  sway  displacement  required  and  roll  rate 
scaling  gain  and  washout  break  frequency  are  held  constant.  The  constraint 
on  peak  sway  displacement  is  satisfied  by  kc  = 2.0  for  a second  order 
washout  circuit.  The  drive  logic  configuration  parameters  used  are  listed 
in  Table  lU. 

Roll  Preference  Series  — The  Roll  Preference  (RP)  Series  drive  logic 
configurations  are  structured  to  investigate  the  tradeoff  between  roll  rate 
scaling  gain  and  washout  break  frequency.  Peak  cab  roll  displacement  required 
is  held  constant  at  that  constant  value  which  allows  the  sway  Mis  displace- 
ment to  remain  within  the  motion  base  limit  for  all  roll  axis  configurations 

TABLE  14.  LAMBS  — VARIABLE  PARAMETERS  — ROLL  AND  SWAY  AXES 


DRIVE  LOGIC 
CONFIGURATION 

KP 

Op  = 0>p 

rad/ sec 

o>bry  = a* 
rad:/  sec 

*c 

03c 

rad/ sec 

Reference  Case 

0.2 

0.7 

0.625 

1.0 

0.5 

Coordination 

Preference 

(Heave  axis 

1 1 l 

at  Reference  Case  settings) 

1 t 1 

CP  0.25 
cp  0.50 

Reference  Case 

Roll  Preference 
(Washed  Out, 
Full  Coordina- 
tion) 

RP  0.05 

RP  0.15 

Reference  Case 
RP  0.30 


0.7  0.625  0.25  0.125 

0.7  0.625  0.50  0.25 

0.7  0.625  1.00  0.5 

(Heave  axis  at  Reference  Case  settings) 


0.05 

0.175 

O.625 

1.00 

0.5 

0.15 

0.525 

0.625 

1.00 

0.5 

0.20 

0.70 

0.625 

1.00 

0.5 

0.30 

1 .05 

0.625 

1 .00 

0.5 

T 


specified.  Coordination  high  frequency  scaling  gain  is  1.0  and  the  coordina- 
tion washout  break  frequency  is  0.5  rad/sec  for  all  roll  washout  configurations. 
The  constraint  on  peak  cab  roll  displacement  is  satisfied  by  Kp  = cop/3.5  for 
a second  order  washout  circuit.  The  drive  logic  configuration  parameters 
are  listed  in  Table  l4. 

Drive  logic  parameter  values  for  the  pitch  and  roll  residual  tilt  paths 
and  for  the  second  order  yaw  rate  washout  circuit  were  fixed  and  not  subject 
to  investigation. 

4.  Subjective  Ratings  and  Protocols 

A sample  pilot  rating  form  is  reproduced  in  Table  15.  One  rating  form 
was  completed  for  each  group  of  runs  using  a single  drive  logic  configuration. 
Typically  all  three  flying  tasks  were  flown  to  evaluate  a single  drive  logic 
configuration. 

The  purpose  of  each  question  on  the  pilot  rating  form  is  summarized 
below. 

Question  1 . — Obtain  indication  of  special  simulator  flying  technique. 
Requires  pilot  to  identify  limitations  imposed  by  simulator. 

Question  2.  — Basic  source  of  pilot  rating  data  for  simulated  motion 
quality/fidelity  vis-a-vis  actual  flight.  Separate  ratings  are  solicited 
for  the  pilot  body-oriented  three  angular  velocity  and  three  specific  force 
axes.  A "one"  rating  corresponds  to  "same  as  actual  flight".  A "five"  rat- 
ing corresponds  to  "aborting  the  flying  task  because  of  nausea  and/or  dis- 
orientation". A "three"  rating  corresponds  to  a "net  no  value"  for  the  motion 
cue  presented.  "Net  no  value"  indicates  either  that  the  positive  value  fac- 
tors offset  the  negative  value  factors  or  zero  value  of  all  factors. 

Question  3«  — A crosscheck  on  ratings  of  Question  2. 

Question  4.  — Separation  of  combined  visual/motion  presentation  factors 
from  motion  presentation  factors  alone. 

Question  3.  — Opportunity  to  interpret  motion  infidelity  as  a disturbance- 
like effect. 
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□ Check  if  you  have  recorded  comments  supporting  above  ratings  on  tape. 

3.  Circle  the  axes  for  which  "no  motion"  would  be  preferred  to  the  motion  provided. 


ROIL 


LATERAL 
G FORCE 


NORMAL 
G FORCE 


PITCH 


LONGI- 
TUDINAL G 
FORCE 


YAW 


U.  Are  the  visual  scene  and  motion  sensations  consistent?  Yes  , No 

Ibssibly  . Describe  any  inconsistencies  


3-  Do  you  detect  turbulence  or  other  disturbances?  Yes 
describe  characteristics  of  the  disturbances.  


_»  *> 


If  Yes. 


6 . If  the  G-Sult  drive  is  active,  does  it  add  to,  or  detract  from  realism  of 
Q-force  sensations? 


ADDS  TO , NO  EFFECT 

Cbmments  on  G- force  sensations 


H)  TRACTS  FROM 


7.  Please  tape  record  other  00 intents  on  task,  fidelity  of  A- 7 representation, 
visual  fidelity,  etc. 





9k 


I ; 


/ 


Quest Ion  6.  — G-suit  effectiveness  evaluation. 

Question  7.  — Separation  of  simulated  motion  factors  from  other  simu- 
lation factors.  Collect  serendipitous  comments . Attract  nonrelevant  comments 
which  might  otherwise  be  made  in  responding  to  Questions  1 through  6. 

Pilot  rating  forms  were  filled  out  in  their  entirety  for  the  Basic  Series 
of  experiments.  In  the  Final  Series,  only  Questions  2 and  3 were  filled 
out. 

During  all  Series,  pilots  were  instructed  to  restrict  their  flying  to 
the  three  prescribed  tasks. 

No  particular  order  was  observed  in  performing  the  three  flying  tasks 
during  the  Basic  Series  of  runs.  Replications  were  performed  as  the  pilots 
deemed  necessary  to  formulate  their  evaluations.  All  evaluations  were 
delivered  in  a single  session. 

In  the  Final  Series,  flying  tasks  A and  C were  flown  alternately.  Three 
replications  for  each  task-drive  logic  configuration  combination  were  flown. 
Strip  chart  and  digital  magnetic  tape  recorded  data  were  collected  on  the 
second  replication.  If  that  trial  was  "blown",  additional  records  were 
collected  on  the  third  replication.  Pilot  comments  were  recorded  on  rating 
forms  and  voice  tape  at  the  conclusion  of  runs  for  each  drive  logic  configura- 
tion. Orders  of  presentation  for  drive  logic  configurations  in  the  Final 
Subseries  are  listed  in  Table  16.  After  completing  each  configuration  the 
pilot  was  reminded  of  the  motion  axis  under  investigation.  He  was  then 
prompted  to  make  a forced  declaration  of  preference  for  the  motion  (in  that/ 
axis)  of  the  last  configuration  or  next-to-last  configuration.  He  was  a^dso 
asked  to  order  the  configurations  in  the  Subseries  according  to  his  subjective 
preference  for  the  motion  provided  by  each. 

5.  Pilot  Subject  Background 

Backgrounds  for  the  three  participating  pilots  are  summarized  in  Table  17. 
All  three  are  test  pilot  school  graduates. 


/ 


95 


>felteite-iri4.-awa 


TABLE  16.  ORDERS  OF  PRESENTATION  FOR  DRIVE  LOGIC  CONFIGURATIONS 
IN  FINAL  SUBSERIES  (G- SUIT- ON  ONLY) 

Heave  Preference  (HP) 

Reference 
HP.  25 
Reference 

Coordination  Preference  (CP) 

Reference 
CP.  50 
Reference 

Roll  Preference  (RP) 

RP.05 

RP.15 

Reference 

RP.30 


E.  DATA  AND  INTERPRETATION 

1 . Fixed- Base  Series 

Monitoring  of  the  drive  logic  commands  to  the  motion  base  servos  showed 
validity  of  the  procedure  for  constraining  peak  values  of  selected  commands 
to  be  approximately  constant.  High  face  validity  of  the  three  flying  tasks 
was  confirmed. 

2.  Basic  Series 

A summary  of  pilot  rating  data  for  the  "Reference"  drive  logic  configu- 
ration and  various  special  cases  of  that  configuration  are  given  in  Table  18. 
No  motion  was  identified  as  preferable  to  the  motion  provided  for  12  out  of 
15  ratings  for  the  "lateral  g- force"  presentation;  for  2 out  of  15  ratings 
for  the  "roll"  presentation  and  for  1 out  of  l4  ratings  for  the  "yaw"  presen- 
tation. This  indicates  that  numerical  ratings  for  the  roll  axis  motion  tend 
to  be  biased  downward.  Numerical  ratings  for  other  axes  tend  to  be  confirmed. 
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TABLE  l8.  PILOT  MOTION  FIDELITY  RATING  DATA  FOR  BASIC  SERIES 
(Reference  drive  logic  configuration  and 
special  cases  as  noted) 


MOTION  IS: 

1 

REALISTIC 

2 

USEFUL 

3 

DISTRACTING 

4 

DISORIENTING 

5 


LATERAL  NORMAL 
ROLL  G FORCE  G FORCE 


LONGITUDINAL 

PITCH  G FORCE  YAW 


I 


2 


3 


4 


5 


Legend: 

O Reference  Cose , Special  Cases: 

□ Normal  Specific  Force  Off, 

V G-Suit  Off,  A Coordination  Off 

. 

TABLE  19.  SUNMARY  OF  PILOT  PREFERENCE  COMMENTS  FOR  BASIC  SERIES 

SUBJECT 

PREFERENCE 

PREFERENCE 

ED 

RH 

BN 

ISSUE 

G-Suit 

ON 

ON 

ON 

Normal  Specific 

OK  ON  or  OFF 

OFF 

OFF 

Force 

Coordination 

ON* 

ON* 

ON* 

* Preference  is  slight 


Pilot  comments  indicate  g-suit  induced  cues  have  very  significant  value. 
This  was  especially  apparent  in  comments  delivered  when  the  subjects  were 
deprived  of  this  cue.  G-suit  cues  have  only  slight  value  when  maneuvering 
in  a nearly  unloaded  condition. 

Pilot  comments  indicate  that  "lateral  g force"  presentation  is  slightly 
improved  by  the  presence  of  coordination,  but  that  the  "lateral  g force" 
presentation  is  poor  with  or  without  coordination. 

Pilot  comments  indicate  that  the  g-suit  induced  cues  considerably  domi- 
nate "normal  g force"  in  perception  of  load  factor.  These  separate  cues  are 
not  perceived  in  an  integrated  manner.  "Normal  g force"  presentation  has 
a nm«n  positive  effect  at  load  factor  onset.  "Normal  g force"  presentation 
has  a moderately  negative  effect  when  a sustained  load  factor  is  pulled  on 
as  for  rapid  roll-in  to  a moderate  g turn,  or  as  in  pulling  off  the  target. 
Pilots  describe  this  negative  effect  as  a "pitch  hop".  (This  is  the  result 
of  rapid  washout  of  the  normal  specific  force  cue . ) Pilots  tend  to  comment 
on  the  "normal  g force"  cue  presentation  in  the  "pitch"  evaluation  in  dis- 
tinction to  the  "normal  g force"  evaluation  column  of  the  rating  form.  Elimi- 
nation of  the  "normal  g force"  cue  presentation  resulted  in  generally  favorable 
pilot  comments.  Most  reflected  improvement  because  of  absence  of  the  "pitch 
hop"  while  seme  reflected  degradation  because  of  absence  of  the  load  factor 
onset  cue.  Pilot  preferences  for  the  special  case  conditions  have  been 
summarized  from  the  comments  in  Table  19« 

Interpretations : 

• Fidelity  of  the  lateral  specific  force  presentation  is 
poor. 

• Roll  axis  fidelity  ratings  are  downgraded  because  of 
poor  lateral  specific  force  fidelity. 

• Roll- sway  coordination  has  some  influence  on  perceived 
lateral  specific  force  fidelity. 

• G-suit  induced  cues  have  a strong  positive  effect  upon 
perception  of  load  factor  in  moderate  amplitude  positive 
g maneuvers. 

• Normal  specific  force  cues  may  actually  detract  from 
simulation  fidelity  for  moderate  amplitude  maneuvers. 
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3.  Final  Series 
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a.  Heave  Preference  Series 

The  Basic  Series  experiments  suggested  that  "normal  g force"  onset  cues 
might  have  significant  value.  This  hypothesis  was  tested  by  comparing  data 
for  a drive  logic  configuration  (HP.25)  which  emphasized  onset  cues,  with 
that  for  the  Reference  configuration.  The  experimental  order  of  presentation 
is  given  in  Table  16.  The  subjects  were  informed  only  as  to  which  motion 
axis  was  the  focus  of  this  investigation.  Subjects  were  asked  to  rank  three 
configurations  by  preference.  Resulting  numerical  ratings  are  summarized 
in  Table  20  and  rankings  are  given  in  Table  21 . Numerical  rating  and  rank- 
ing for  the  "pitch"  and  "normal  g force"  presentation  evaluations  both  show 
a clear  preference  for  the  Reference  configuration.  It  is  remarkable  that 
2 of  the  3 subjects  identified  the  first  and  third  configurations  as  being 
the  same  even  though  there  was  no  reason  for  them  to  expect  two  presentations 
of  the  same  configurations. 

Time  histories  in  Fig.  2k  show  the  increased  emphasis  given  normal  specific 
force  (azM)  onset  by  configuration  HP.25  relative  to  that  for  the  Reference 
configuration. 

\ 

Interpretation: 

• No  normal  specific  force  motion  is  preferable  to  that  which 
could  be  provided  within  the  constraints  of  these  experiments. 

This  result  is  based  on  both  the  Heave  Preference  and  Basic 
Series  findings.  The  rapid  washout  which  must  be  used  in 
normal  specific  force  drive  logic  channels  is  responsible 

for  this  finding. 

• Normal  specific  force  onset  cues  have  some  positive  value, 
but  this  is  largely  overshadowed  by  the  negative  effects 
of  rapid  washouts  for  abrupt  maneuvers  of  moderate  g level. 

b.  Coordination  Preference  Series 

The  Basic  Series  experiments  showed  washed  out  coordination  was  preferred 
slightly  to  no  coordination.  The  purpose  of  the  Coordination  Series  is  to 
determine  if  there  is  a favorable  tradeoff  between  reduced  high  frequency 
coordination  path  gain  and  reduced  coordination  washout  break  frequency. 
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TABLE  20.  PILOT  MOTION  FIDELITY  RATING  DATA  FOR 
HEAVE  PREFERENCE  SERIES 


. 

> 

I 


4 
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MOTION  IS: 

1 

REALISTIC 

2 

USEFUL 

3 

DISTRACTING 

4 

DISORIENTING 

5 


LATERAL  NORMAL  LONGITUDINAL 


ROLL  G FORCE  G FORCE  PITCH  G FORCE  YAW 


I 


2 


3 


4 
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Configurations:  © Reference,  Gf  HP-25 


TABLE  21 . PILOT  PREFERENCE  RANKINGS  FOR  HEAVE 


RANK 

preference 

ED 

SERIES  CONFIGURATIONS 

SUBJECT 

RH 

BN 

1 

1st  Ref. 

Last  Ref. 

Last  Ref. 

2 

Last  Ref. 

1st  Ref. 

1st  Ref. 

5 

HP.  25 

HP. 25 

HP.  25 

I 

This  hypothesis  was  tested  by  comparing  data  for  a drive  logic  configuration 
(CP. 50)  and  the  Reference  configuration.  Configuration  CP.50  reduces  both 
high  frequency  gain  and  washout  break  frequency  to  one-half  the  values  in 
the  coordination  path  for  the  Reference  configuration.  Ifoe  experimental 
order  of  presentation  is  given  in  Table  16.  The  subjects  were  informed 
only  as  to  which  motion  axis  was  the  focus  of  this  investigation.  Subjects 
were  asked  to  rank  three  configurations  by  preference.  Resulting  numerical 
ratings  are  summarized  in  Table  22  and  rankings  are  given  in  Table  23. 

Time  histories  in  Fig.  25  show  that  the  level  of  spurious  lateral  specific 
force  is  increased  for  configuration  CP.50  with  respect  to  the  Reference  con- 
figuration. Fhasing  of  the  presented  lateral  specific  force,  ayM,  with  the 
actual  in-flight  quantity,  ay^,  is  seen  to  be  good  for  both  drive  logic 
configurations.  Since  the  input  scaling  gain  on  lateral  specific  force 
commands  is  Ky  = 0.5,  the  ayiN  and  ayM  traces  are  directly  comparable  for 
the  recording  sensitivities  used. 

Pilot  numerical  rating  and  ranking  data  and  pilot  comments  generally 
indicate  a slight  preference  for  the  configuration  providing  full  coordina- 
tion at  high  frequencies  (i.e.,  k^  * 1.0).  Two  of  the  three  subjects  indicated 
spurious  "lateral,  g forces"  were  least  for  this  configuration.  The  third  sub- 
ject comments  indicated  his  preference  for  configuration  CP.50  over  the  Ref- 
erence configuration  was  slight. 

Interpretation : 

• There  is  at  least  slight  pilot  preference  for  full  roll- 
sway  coordination  at  high  frequencies. 

• Experimental,  conclusions  are  weak  because  the  maximum 
level  of  low  frequency  coordination  available  provided 
marginal,  or  inadequate  roll-sway  axis  motion  fidelity. 

c.  Roll  Preference  Series 

The  purpose  of  the  Roll  Preference  Series  is  to  eaqslore  the  tradeoff 
between  roll  washout  scaling  gain  and  break  frequency.  These  results  are 
also  key  in  establishing  the  marginal  utility  (op)  parameter  value.  Hie 
order  of  presentation  of  the  configurations  in  the  experiment  is  given  in 
Table  16.  The  subjects  were  informed  only  as  to  which  motion  axis  was  the 
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TABLE  22.  PILOT  MOTION  FIDELITY  RATING  DATA  FOR 
COORDINATION  PREFERENCE  SERIES 


MOTION  IS: 

1 

REALISTIC 

2 

USEFUL 

3 

DISTRACTING 

4 

DISORIENTING 

5 


LATERAL 
G FORCE 


NORMAL 
G FORCE 


PITCH 


LONGITUDINAL 
G FORCE 


Of 

o 

0 

O 

O 

Of 

Configurations;  ©Reference,  gTcP.  50 


TABLE  23.  PILOT  PREFERENCE  RANKINGS  FOR  COORDINATION 
PREFERENCE  SERIES 


ED 

1 st  Ref. 
Last  Ref. 
CP.  50 


SUBJECT 

RH 

CP.  50 

1st  Ref. 

Last  Ref. 


BN 

1st  Ref. 
Last  Ref. 
CP.  50 


• > -T*  ' 

_ . • 


I 


Reference 
0 ®®  ® 
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focus  of  this  investigation.  Subjects  were  asked  to  rank  four  configurations 
by  preference.  Resulting  numerical  ratings  are  summarized  in  Table  24,  and 
rankings  are  given  in  Table  25.  Numerical  rating  and  ranking  for  the  "roll" 
and  "lateral  g force"  presentations  show  configuration  RP.30  to  be  the 
least  preferred.  Preferences  otherwise  tend  to  favor  the  small  scaling 
gain/low  break  frequency  combinations,  RP.05  and  RP.15.  (Configuration 
RP.18  is  the  result  of  an  error  in  setting  up  the  RP.05  configuration  for 
subject  RH. ) The  setting  for  Op  was  1.05  rad/sec  instead  of  0.175  rad/sec. 
op  was  set  correctly  to  0.175  rad/sec.  The  effective  lead  break  frequency 
is  1.225  rad/sec  which  is  between  that  for  the  Reference  configuration 
(l.4o  rad/sec)  and  the  RP. 15  configuration  (l.05  rad/sec).  The  effective 
lead  break  frequency  is  computed  using  the  Appendix  D method. 

Time  histories  in  Fig.  26  show  the  increased  roll  washout  rate  effects 
as  one  proceeds  through  configurations  RP.05,  RP.15,  Reference  and  RP.30. 
Notice  also  that  peak  cab  roll  displacement  is  approximately  constant  for 
all  configurations. 

Interpretation : 

• Second  order  roll  washout  break  frequencies  up  to  1.0  rad/sec 
are  acceptable.  Only  at  1.0  rad/sec  do  roll  washout  effects 
begin  to  show  up  in  pilot  comments. 

• There  is  an  approximately  even  tradeoff  between  roll  wash- 
out scaling  gain  and  break  frequency  (below  1.0  rad/sec) 
when  the  relationship  between  the  two  parameters  is  deter- 
mined by  holding  peak  cab  roll  angle  constant. 

F.  SIMULATION  SQUAWKS 

Three  residual  problems  existed  throughout  this  simulation  and  could  not 
be  eliminated: 

• A- 7 equations  of  motion  produce  strong  proverse  yaw  while 
in  fact  the  actual  aircraft  has  an  adverse  yaw  character- 
istic at  these  flight  conditions. 

• The  motion  base  produces  a "shuddering"  motion  in  sway 
when  the  sign  of  cab  roll  rate  changes. 

• Occasional,  very  small,  sway  axis  glitches  give  subjects 
the  impression  of  small  rudder  kick  inputs. 
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TABLE  2k.  PILOT  MOTION  FIDELITY  RATING  DATA 
FOR  ROLL  PREFERENCE  SERIES 


MOTION  IS: 

ROL 

REALISTIC 

2 — 

A O 

USEFUL  © ^ 

3 — 

©a 

DISTRACTING 

DISORIENTING  © 
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G FORCE  G FORCE 
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PITCH  G FORCE 


OAO  V I 


0 

© 0©V 
©A 


OA0DV 
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Configurations:  ©Reference,  A RP.05,  0 RP.I5,  B RP.I8,  V RP.30 


TABLE  25.  PILOT  PREFERENCE  RANKINGS  FOR  ROLL  PREFERENCE 
SERIES  CONFIGURATIONS 

SUBJECT 


RANK 

ED 

RH 

BN 

1 

RP.15 

j RP.15 1 

RP.05 

2 

( RP.05  ) 

I RP.18J 

REF. 
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( RP.30  ) 

REF. 

RP.15 
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• Sensation  of  sitting  on  top  of  a rotating  drum,  i.e., 
significantly  above  the  pilot's  location  in  the  actual 
aircraft. 

The  first  comment  is  easily  explained  because  the  pilot  is  indeed  tilted 
toward  the  inside  of  the  "turn".  To  the  extent  that  simulator  cab  tilt  is 
not  coordinated  with  simulator  cab  sway  acceleration,  the  pilot  is  exposed 
to  a spurious  lateral  specific  force,  (—  [gcfty  — Y^] ) which  is  exactly  that 
required  to  produce  the  sensation  of  falling  to  the  inside  of  the  turn. 

The  second  comment  arises  from  the  same  aspect  of  lateral  specific  force 
reproduction,  but  requires  interpretation  from  a dynamic  rather  than  a static 
viewpoint.  Consider  the  following  block  diagram  which  shows  components  of 
lateral  specific  force  arising  from  cab  tilt,  ideal  sway  coordination  and 
the  missing  component  of  sway  coordination  which  must  be  endured  as  a prac- 
tical matter  because  of  limited  simulator  sway  displacement  capability. 

The  dashed  path  in  Fig.  27  represents  the  kinematic  effect  of  locating  the 
pilot's  seat  above  the  cab  roll  axis  when  (— iz)>0.  Consider  the,  (— iz)  value 
giving  equivalence  between  the  "missing  coordination  component"  and  the 

"sitting  on  top  of  a drum"  paths  at  an  arbitrary  frequency,  <u.  This  is 

—2 

obtained  by  evaluating  the  frequency  response  for  s times  the  transfer 
function  for  the  "missing  coordination  component"  path  at  the  frequency,  cu. 

gVl  + (2^/tuc)2 

(-iz)  = 


^[l  “ (o>/o>c)2]2  + [?^0)/03c]2 


Actual 

Coordination 


Perfect 

Coordination 


Missing  Component 
of  Sway  Coordination 


(s+u>c/2£) 
s*  + 2C‘»cS+«»! 


Figure  27.  Block  Diagram  for  Coordination  Effects 


The  values  of  (— iz)  which  result  clearly  justify  the  subjective  impression 
of  "sitting  on  top  of  a drum"  at  low  and  mid  frequencies. 


H.  CONCLUSIONS 


G-suit  cue  dominates  any  motion  which  can  be  provided  by 
the  simulator.  Subjects  strongly  prefer  to  use  the  g-suit 


• G-suit  pressure  schedule  could  be  improved  for  better 
unloaded  flight  sensation.  1 g pressure  could  be 
increased,  gradient  could  be  increased. 


• G-suit  cue  is  not  integrated  with  motion  cue  in  per- 
ception. G-suit  cue  essentially  "replaces"  z axis 
motion  cue  roles  in  the  flying  task. 

• Near-zero  g and  negative  .r  simulation  sensations  are 
missing. 

• Pilot  preference  in  a z axis  motion  scaling  versus 
washout  rate  tradeoff  (with  simulator  heave  axis  dis- 
placement held  constant)  tends  to  favor  zero  seeding 
gain  and  zero  washout  rate.  That  is,  no  z-axis  motion 
is  preferred.  This  result  presumes  the  g-suit  is  used. 

• Spurious  z-axis  motion  washout  effects  are  quite 
apparent  for  moderate  amplitude  maneuvers  such  as 
are  used  in  air-to-ground  flying  tasks  (e.g.,  tight 
turn  entry,  pull-off).  Subjectively  these  effects 
sire  identified  as  a "pitch  hop,"  a rapid  washout  of 
the  normal  specific  force. 


Roll  Axis 

• Roll  rate  scaling  and  washout  rate  do  not  appear  to  have 
a direct  sensitive  influence  upon  angular  motion  percep- 
tion. Effects  are  secondary  because  of  the  influence 

of  these  two  characteristics  upon  the  degree  of  coordina- 
tion which  can  be  achieved. 

• Roll  rate  scaling  gains  over  the  range  0.05  to  0.50  (with 
washout  rate  set  to  give  constant  peak  cab  roll  displace- 
ment, and  with  constant  coordinating  crossfeed  character- 
istics) did  not  result  in  pilot  comments  indicating  per- 
ceptual sensitivity  to  roll  rate  scaling. 

• Second-order  roll  rate  washout  break  frequencies  ranged 
from  0.175  to  1.05  rad/sec  in  a manner  consistent  with 
the  roll  rate  scaling  gains  to  achieve  constant  peak 
cab  roll  displacement.  Only  at  the  most  rapid  washout 
rate  (1.O5  rad/sec)  did  roll  washout  effects  begin  to 
show  up  in  pilot  comments.  Washout  break  frequencies 
up  to  1.0  rad/sec  therefore  seem  acceptable.  Break 
frequencies  should  be  limited  to  this  value  by  appro- 
priate use  of  scaling  gain. 


Lateral  Specific  Force  Axis 


• Moderate  amplitude  maneuvers  having  rapid  onset  charac- 
teristics require  full  roll- sway  simulator  coordination 
at  high  frequencies. 
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A second-order  washout  break  frequency  of  0.5  rad/sec 
was  found  to  be  marginally  acceptable  for  the  roll- 
sway  coordination  path.  A still  lower  value  is 
desirable.  Lack  of  roll-sway  coordination  at  lower 
frequencies  tends  to  give  the  pilot  the  impression 
that  he  is  sitting  further  above  the  aircraft  roll 
axis  than  is  actually  the  case. 


Simulator  lateral  g response  to  lateral  g step  commands 
should  be  monotonically  increasing  with  time  in  order 
to  avoid  dynamic  cue  infidelity  for  rudder  inputs.  A 
residual  tilt  wash- in  first-order  break  frequency  of 
0.5  rad/sec  results  in  acceptable  tilt  rates. 


Longitudinal  Specific  Force  Axis 


Residual  tilt  representation  of  longitudinal  g forces 
was  adequate  as  long  as  there  was  no  rapid,  large  appli- 
cation of  thrust  or  speed  brakes.  Rapid,  large  applica- 
tions resulted  in  excessive  tilt  rates  and  inadequate 
sensation  of  acceleration.  A residual  tilt  wash- in 
first-order  break  frequency  of  0.5  rad/sec  is  a good 
compromise  for  the  air-to-ground  flying  task. 


Pitch  Axis 


• Generally  adequate,  not  a subject  of  this  investigation. 
Yaw  Axis 

• Generally  adequate,  not  a subject  of  this  investigation. 
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SECTION  VI 

UTILIZING  THE  MEASURES  AND  CRITERIA  — 

LAMARS  OPTIMIZATION 

A.  INTRODUCTION 

The  results  obtained  from  the  simulator  experiment  described  in  the 
previous  section  are  used  to  further  refine  the  criteria  for  motion  fidelity- 
in  the  various  axes.  In  the  case  of  the  yaw  and  heave  axes,  this  simply 
amounts  to  inferring  appropriate  values  for  the  marginal  utility  parameters 
on  the  basis  of  the  experimental  results.  These  marginal  utility  parameters 
are  selected  to  cause  trends  in  the  criterion  values  consistent  with  pilot 
subjective  evaluations  elicited  as  the  experimental  variables  were  changed. 

For  the  pitch,  surge,  roll  and  sway  axes  it  turns  out,  the  experimental 
data  show  that  the  original  criteria  are  inadequate  for  defining  the 
subjective  optimum  reproduced  motion.  The  roll-sway  and  pitch-surge  axes 
must  be  considered  as  coupled  sets.  This  necessitates  development  of  a 
method  for  combining  the  criteria  for,  say,  the  roll  and  sway  axes  into  a 
single  criterion  for  the  coupled  roll-sway  set.  A plausible  combination, 
suggested  by  the  variation  of  the  measures  as  a function  of  the  experimental 
variables  and  by  the  experimental  results,  is  developed  in  parametric  form. 

Once  the  final  parameterized  forms  for  the  revised  criteria  are  chosen, 
the  measures  generated  in  the  parametric  study  are  used  to  evaluate  the 
criteria  for  each  axis.  Assumed  values  for  the  parameters  and  the  particular 
values  for  the  experimental  variables  used  in  the  simulator  experiments  are 
used  in  these  evaluations.  For  each  axis,  the  appropriate  criterion  (J) 
value s were  computed  as  functions  of  the  free  parameters  (a,  (3  and  X.).  These 
three  constants  are  then  chosen  such  that  trends  in  the  criterion  values  with 
changing  experimental  variables  match  the  pilot  subjective  evaluations  elici- 
ted during  the  simulator  experiment. 

Once  the  a,  0 and  X parameter  values  are  determined,  the  criteria  develop- 
ment is  complete,  and  the  resulting  criteria  are  ready  for  use  for  optimizing 
any  particular  motion  base  configuration  for  any  given  task.  In  this  study, 
the  LAMARS  motion  base  is  optimized  for  the  air-to-ground  scenario.  Also, 
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two  nonlinear  washout  features  added  to  the  basic  LAMARS  framework  are  con- 
sidered briefly  once  the  basic  LAMARS  drive  logic  is  optimized. 

Obviously,  the  method  described  in  this  chapter  is  based  upon  a rela- 
tively few  simulator  experimental  conditions  explored  in  this  limited  scope 
program.  The  significance  of  the  results  transcends  the  limitations 
imposed  by  this  modest  experimental  data  base,  however.  This  is  so  because 
the  main  contribution  is  the  demonstration  of  an  orderly  method  for  pre- 
dicting the  effectiveness  of  simulating  motion  for  a particular  task  on  a 
particular  motion  base,  and  an  orderly  procedure  for  optimizing  the  available 
parameters  of  the  simulation  drive  logic.  This  method  and  procedure  may  be 
applied  readily  by  simulation  engineers. 

The  following  sections  discuss  revision  of  the  criteria,  the  LAMARS 
optimization  and  the  investigation  of  the  Parrish  and  subliminal  nonlinear 
washout  features. 

B.  REVISED  MEASURES  AND  CRITERIA 

The  results  of  the  simulator  experiment  suggest  reexamination  of  the 
measures  and  criteria  developed  earlier  in  order  to  more  comprehensively 
reflect  motion  fidelity  requirements  for  the  two  sets  of  coupled  axes. 

This  necessity  arises  mainly  in  connection  with  the  translational  washouts 
used  in  the  surge  and  sway  axes.  As  an  example,  consider  the  equation  for 
the  original  sway  axis  criterion  computed  using  the  post-drive  logic  measures 
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This  criterion  treats  lateral  specific  force  reproduction  fidelity  vir- 
tually independently  of  roll  angular  velocity  reproduction  fidelity.  How- 
ever, it  is  clear  that  roll  axis  scaling  and  washout  rate  strongly  influence 
lateral  specific  force  fidelity  and  vice  versa  via  requirements  for  coordi- 
nation. Furthermore,  low  frequency  representations  of  lateral  specific 
force  achieved  via  residual  tilt  results  in  spurious  roll  angular  velocity 
cues.  These  two  facts  make  it  clear  that  a single  criterion  must  he  for- 
mulated for  the  roll  and  sway  axes  considered  together.  This  criterion 
must  also  combine  all  key  and  usually  competing  demands  upon  the  simulator 
for  motion  fidelity  in  these  two  axes.  Similarity  of  the  pitch  and  surge 
axes  to  the  roll  and  sway  axes  makes  the  above  remarks  appropriate  for  the 
pitch-surge  axes  as  well. 

Consider  the  following  criteria  for  the  roll-sway  and  pitch-surge  axes 
which  are  responsive  to  the  above  requirements. 

Jyp  = [l  — p]  Jy  + P([l  — X]  Jp  + X Jyy)  (l) 

JXq  = [1  - f3]  JX  + P([1  - X]  Jq  + X JXT)  (2) 

Each  criterion  is  a combination  of  three  subcriteria.  Each  subcriterion 
emphasizes  one  particular  aspect  of  motion  fidelity.  The  free  parameters, 

P and  X,  allow  us  to  tailor  the  relative  weights  given  the  subcriteria  in 
order  that  the  combination  correspond  with  subjective  pilot  evaluations. 

The  particular  form  of  the  right  hand  sides  for  Eq.  1 and  2 assures  that: 

If  Jy  = Jp  = Jyp  — 1 then  Jyp  = 1 

and 

If  Jy  = Jp  = Jyy  = 0 then  Jyp  = 0 

thereby  maintaining  the  original  limiting  cue  properties  for  these  new 
criteria.  Next  consider  the  subcriteria.  Jyy  is  new,  and  Jy  is  modified 
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from  its  original  definition.  Jy^  represents  motion  reproduction  error 
effects  of  the  lateral  specific  force  response  to  lateral  specific  force 
command  inputs  in  the  absence  of  roll  rate  commands  (i.e.,  with  Kp  = 0). 


JyT  = (l  - K2r)  + 


+ 2k(  1 - Kgp) 


Vpo  ~ ^poF^Kp-p 


K§R 


(3) 


^2R  Vpo  - (vpo  VpoF)Kp=p 


K2R 


Jy  represents  sway  motion  reproduction  error  effects  for  the  sway  axis  in 
response  to  coordination  and  lateral  specific  force  command  inputs.  The 
error  upon  which  Jy  is  based  is; 

€s  = vpo  “ VpoF  = vpo  - (y  - g<ftn)  = (g<ftn  + ^po)  “ y 

■ 0 ” Ky)(g9m  + Vpo)  + KyCgcpm  + vpo)  - y 

Remaining  steps  in  the  development  of  Jy  are  the  same  as  those  used  earlier 
in  Section  IV. 


Jy  = (i  - Ky)2  + k2 


+ 2K(l  - Ky) 


(Kytff&n  * *!*>)  “ y)  Kp*0 

(Xytg^m  + ^po^  ~ y)2  “ (*y + tpo]  - y)2p.o  j 

Kp(g <Pm  + ♦po)2 

K£(g<Pm  + tp©)2  J) 

Ky(g<rm  + Vpo)2-C(gfPm  + Vpo)y]Kp,ol 


Ky(g<ftn  + Vpo)^ 


c(g«tw  + y>)y'iiCp-o  - f(gTw  * y)y|| 


The  third  subcriterion,  Jp,  retains  its  original  definition.  Jp  represents 
motion  reproduction  error  effects  for  the  roll  axis  in  response  to  roll  rate 
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and  lateral  specific  force  command  inputs.  (Lateral  specific  force  command 
inputs  produce  a false  roll  rate  cue  response  because  residual  tilt  is  used 
to  provide  low-frequency  lateral  specific  force  reproduction.) 


Corresponding  definitions  apply  for  Jxt>  Jx>  and  Jq  respectively.  The 
re stilting  criteria,  Jyp  and  JXq,  now  embody  terms  representing  all  effects 
which  potentially  may  compromise  motion  fidelity  in  the  roll- sway  and  pitch- 
surge  axis  respectively. 


The  final  step  in  revision  of  the  criteria  is  to  limit  the  measures  used 
for  all  criteria  (other  than  for  the  heave  axis  with  the  g-suit  active)  to 
those  computed  from  post-drive  logic  variables.  The  main  reason  for  this 
decision  is  the  finding  that,  almost  without  exception,  post-perceptual 
processing  and  post-threshold  processing  measures  are  little  different  and 
are  frequently  larger  than  corresponding  post-drive  logic  quantities.  This 
disappointing  finding  is  contrary  to  a key  initial  premise.  Namely,  it  had 
been  presumed  at  the  outset  that  post-perceptual  processing  and  post- threshold 
processing  measures  would  be  significantly  smaller  than  corresponding  post- 
drive logic  measures.  This  turned  out  to  be  false. 


There  is  one  notable  exception  to  the  above  finding.  The  post-perceptual 
processing  measures  which  incorporate  perceptual  integration  of  normal  specific 
force  motion  and  g-suit  tactile  cues  show  marked  reduction  with  respect  to 
corresponding  post-drive  logic  measures  (which  do  not  include  g-suit  effects). 
However,  the  post-perceptual  processing  measures  without  g-suit  effects 
approximate  the  corresponding  post-drive  logic  measures.  This  result  is 
consistent  with  the  experimental  subjective  opinion  that  g-suit  cues,  when 
present,  dominate  normal  specific  force  cues  for  the  air-to-ground  flying 
task.  Post-perceptual  processing  measures  are  used  in  the  heave  axis  cri- 
terion when  the  g-suit  is  active. 


A summary  of  the  revised  criteria  and  the  subcriteria  used  in  Eqs.  1 and 
2 is  given  in  Table  26. 


Now  that  the  measures  and  criteria  are  in  their  final  form,  the  last  step 
in  the  criteria  development  process  can  proceed.  In  this  step  values  for 
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TABLE  26.  REVISED  LIST  OF  NEASURES 


the  free  parameters  a,  p,  and  X are  found  for  each  axis  (or  set  of  axes) 
such  that  the  criterion  values  for  a given  set  of  experimental  variables 
reflect  the  trends  in  subjective  pilot  rankings  obtained  when  the  same 
values  for  the  experimental  variables  were  used  in  the  simulator  experiment. 
This  step  is  carried  out  by  constructing  plots  of  the  various  criteria  and 
subcriteria  as  a function  of  the  principal  experimental  variable,  washout 
break  frequency,  and  with  a,  p and  X as  parameters,  a,  p and  X values 
are  then  selected  to  cause  the  respective  criterion  values  to  display  the 
same  optimum  and  trends  as  the  subjective  pilot  evaluations. 

1 . Heave  Axis 

Results  from  the  simulator  experiment  for  the  heave  axis  are  complicated 
by  the  use  of  the  g-suit  to  simulate  normal  specific  force  cues.  To  begin, 
compare  the  Jz  criterion  computed  with  the  post-drive  logic  measures  for  four 
different  values  of  a,  using  Fig.  28  . The  object  here  is  to  match  the  rela- 
tive values  of  Jz  across  the  range  of  break  frequencies  with  the  subjective 
rankings  of  the  pilots  during  the  simulator  experiment.  In  the  simulator, 
the  ratings  showed  a preference  for  the  low  break  frequencies  and  a dislike 
for  the  higher  break  frequencies  — particularly  those  greater  than  1.0. 

The  Jz  values  for  = 1.2  present  the  best  match  to  the  pilot  preferences. 

Since  the  heave  axis  is  uncoupled  from  the  others,  is  the  only  con- 
stant which  needs  to  be  determined  in  this  first  optimization  step.  It  is 
interesting,  however,  to  examine  the  effect  of  the  g-suit  on  the  criteria 

i 

values.  Figure  29  presents  a plot  of  JZg-suit  vs*  break  frequency.  Note 
here  that  Kz  = 1.0  for  all  values  of  the  break  frequency,  this  is  because 
the  gain  of  the  g-suit  perception  is  presumed  to  be  unity.  The  values  com- 
puted for  JZg-suit  indicate  that  the  g-suit  is  essential  for  realistic  simu- 
lation of  normal  specific  force  cues;  the  JZg-suit  vaiues  8116  at  least  an 
order  of  magnitude  less  than  the  Jzno  g-suit  values  shown  in  Fig.  50.  This 
observation  is  in  total  agreement  with  the  subjective  pilot  ratings  con- 
cerning the  usefulness  of  the  g-suit  cues. 


♦Tables  of  measures  used  to  compute  the  J values  throughout  this  chapter 
can  be  found  in  Appendix  C. 
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2.  Roll- Sway  Axes 


For  this  coupled  axis  set  the  procedure  is  two-fold  — first,  the  a 
values  for  each  axis  must  be  determined  separately.  Once  this  is  accomp- 
lished, the  J values  must  be  combined  to  determine  the  p and  X values.  At 
each  stage  it  is  necessary  to  match  the  shape  of  the  criterion  curve  with 
the  subjective  pilot  ranking  trends. 

Figure  31  is  a plot  of  Jp  vs.  break  frequency  for  several  values  of 
ap.  The  pilot  rankings  presented  in  the  previous  chapter  ordered  the  break 
frequency  settings  in  the  following  way:  O.525  (best),  O.175,  0.70,  1.05 
(worst),  op  = 1.2  value  best  approximates  the  relative  ordering  of  the 
Kp,  Bp  combinations  as  determined  by  the  simulator  pilots. 

The  J y criterion  must  be  evaluated  as  a function  of  ay  and  C2y>  since 
Bp  was  held  constant  in  t..e  experiments  while  C2y  was  varied.  The  values 
used  for  the  other  parameters  are  indicated  in  Fig.  32  — note  that  the  gain 
Ky  is  here  a function  of  Cpy  The  pilot  rankings  of  the  Cpy  values  were  as 
follows:  0.25  (best),  0.0625,  0.0  (worst).  Any  value  for  c&y  less  than  0.6 
orders  the  Cpy  values  properly  so  ay  = 0.5  was  chosen  somewhat  arbitrarily. 
Figure  33  presents  a plot  of  the  Jy  criterion  versus  Bp  and  ay  for  a value 
of  C2y  = 0.25;  this  plot  merely  serves  as  a check  of  the  behavior  of  Jy  vs. 
Bp  for  the  values  of  ay  used  in  Fig.  32.  These  values  of  ay  will  also  be 
used  in  the  Jyp  computation. 

Only  a single  value  of  Jyp  was  computed,  since  there  was  no  simulator 
data  upon  which  to  base  a selection  of  ayp.  The  best  guess  at  this  value 
is  ayp  = 1.0.  Then,  assuming  Bp  = 0.5,  KpR  = 0.5  and  C2y  = 0.25, 

Jyp  = (l  — ^2r)^  + (aYT^2R)^FDLl  ] + aypKgpO  — K2r)[FDL3] 

= 0.628 


Now  that  the  a' s have  been  selected  for  the  three  subcriteria,  the  appro- 
priate J' s can  be  combined  to  form  the  Jyp  criterion. 

Figure  3^  is  the  plot  of  Jyp  versus  Bp  for  three  values  of  pyp  and  Xyp* 
According  to  the  pilot  rankings  of  the  roll- sway  axis  set  fidelity,  the 
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appropriate  combination  of  Byp  811(1  *yp  woul<i  be  one  which  places  a "minima" 
in  the  JYp  curve  at  a Bp  of  approximately  0.5  rad/sec.  Prom  Fig.  3h,  the 
curve  which  best  fits  the  relative  rankings  of  the  simulator  pilots  and  con- 
tains a minima  of  Jyp  at  Bp  = 0.5  rad/sec  has  Pyp  “ 0.45  and  *-yp  “0.1.  Note 
that  many  other  combinations  of  Byp  and  Xyp  were  tried  before  selecting  these 
values.  For  the  sake  of  clarity,  only  three  values  are  plotted  in  Fig.  34. 
These  three  values  came  closest  to  matching  the  subjectively  optimum  experi- 
mental conditions. 

3.  Pitch-Surge  Axes 

i 

Very  little  data  is  available  from  the  piloted  simulation  experiments 
for  the  pitch-surge  axes.  For  this  reason,  it  was  assumed  that  the  a values 
for  Jy,  Jp  and  Jyip  would  be  the  same  as  those  for  Jy,  Jq  and  Jytj>,  respectively. 
To  examine  the  appropriateness  of  this  assumption,  the  Jq  criterion  was  plotted 
versus  B , for  several  values  of  oq,  as  shown  in  Fig.  35.  In  comparing  this 
plot  to  the  Jp  plot  in  Fig.  31 > no  significant  difference  is  found.  Thus, 
the  assumption  that  ap  = oq  seems  reasonable. 

As  was  the  case  for  Jyj,  only  a single  value  of  Jy»p  is  computed: 

For  ^XT  = 1 Bq  — 0«5>  Kjip  = 0.5 

JXT  = (1  - K2P)2  + (axTK2p)2[PDLl]  + ofcTK2p(l  - K2P)[FDL3] 

= 0.392 

Also,  a single  value  is  computed  for  Jy,  since  in  this  special  case  (LAMARS 
does  not  have  a surge  motion  degree  of  freedom),  Kx  = 0.  Thus,  Jy  = 1.0. 

The  assumption  that  cgc  = ay  = 0.5  is  still  asserted  but  cannot  be  validated 
in  this  study. 

Again,  in  the  absence  of  simulator  data,  it  is  assumed  that  BXq  = Byp  and 
XyQ  = Xyp«  Figure  36  is  a plot  of  JyQ  vs.  Bq  for  three  values  of  BxQ  and 
XyQ  — the  assumed  values  of  OyQ  = 0.45  and  XyQ  = 0.1  provide  a JyQ  curve 
with  a shallow  minimum  between  Bq  values  of  0.25  and  0.5  rad/sec. 
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4.  Ymr  Axis 

The  yaw  axis,  like  the  pitch- surge  axes,  has  little  simulator  data  upon 
which  to  base  a choice  of  or.  The  assumption  is  again  made  that  or  = ctp  = cxq 
with  the  notion  that  the  angular  velocity  axes  can  in  some  sense  be  treated 
similarly.  Figure  37  presents  a plot  of  JR  versus  Br  for  several  values 
of  or,  which  compares  favorably  with  the  plots  of  Jq  vs.  Bq  and  Jp  vs.  Bp, 
shown  in  Figs.  35  and  31,  respectively.  Since  the  yaw  axis,  like  the  heave 
axis,  is  uncoupled,  or  is  the  only  parameter  to  be  determined,  or  = 1.2  is 
chosen. 

Table  27  lists  the  values  of  a,  B and  X for  each  axis  or  set  of  axes 
which  have  been  selected  in  the  process  of  matching  the  subjective  optimum 
and  relative  rankings  of  the  simulator  experiments.  The  values  of  the  cri- 
teria produced  when  these  parameter  values  are  used  suggest  the  same  optimum 
and  relative  ranking  of  the  configurations  as  did  the  experiments. 


TABLE  27.  SUMMARY  OF  CRITERIA  CONSTANTS 


PARAMETER 


VALUE 


Marginal  Utility 
Parameters 


Coupled  Criteria 
Parameter 


Pitch- Surge 
Roll- Sway 
Heave 
Roll- Sway 
Pitch- Surge 
Yaw 

Pitch- Surge 
Roll- Sway 


Pitch- Surge 
Roll- Sway 
Pitch- Surge 
Roll- Sway 


Figure 
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Given  the  criteria  forms  and  a,  3 and  X values  of  the  preceding  sub- 
section, a flying  task  and  a set  of  motion  simulator  hardware  constraints, 
it  is  now  possible  to  optimize  the  drive  logic  for  that  motion  base.  Con- 
sider the  air-to-ground  flying  task  and  the  LAMARS  motion  base  and  drive 
logic  as  a specific  example. 

The  motion  base  hardware  constraints  are  satisfied  via  scaling  selection. 
Scaling  gains  are  chosen  by  comparing  the  peak  displacement  values  obtained 
at  each  value  of  washout  break  frequency  of  interest  with  the  hardware  con- 
straint for  each  axis.  Gains  are  then  computed  for  each  break  frequency  so 
that  the  peak  displacement  will  not  exceed  the  hardware  constraints.  This 
identifies  those  scaling  gain -break  frequency  parameter  combinations  which 
may  be  considered  further. 

For  each  parameter  combination,  then,  the  appropriate  J value  can  be 
calculated  and  compared  with  the  J values  for  all  other  parameter  combina- 
tions for  each  axis  or  coupled  sixes.  The  optimum  parameter  combination  for 
each  axis  corresponds  to  the  minimum  J value.  If  all  of  the  J values  for 
a given  axis  are  greater  than  or  equal  to  1.0,  this  indicates  that  simulated 
motion  in  that  axis  does  not  contribute  to  the  simulation  fidelity. 

The  following  subsections  discuss  the  LAMARS  optimization  for  each  axis 
or  coupled  axis  set:  heave,  roll-sway,  pitch-surge  and  yaw. 

1.  Heave  Axis 


Figure  58  is  a plot  of  vs.  break  frequency.  This  washout  break 
frequency  and  scaling  gain  are  the  only  parameters  varied  for  this  axis. 

The  g-suit  is  presumed  not  used  for  Fig.  58.  Kz  values  are  computed  such 
that  the  peak  displacement  never  exceeds  10  ft.  These  values  are  also  given 
in  the  figure.  Since  all  of  the  Jz  values  are  greater  than  1.0,  it  is  appar- 
ent that  the  perceived  normal  motion  fidelity  is  not  Improved  with  the  use 
of  simulator  motion.  Therefore,  the  optimum  drive  logic  parameters  for  the 
heave  axis  are  both  0.0,  as  indicated  in  the  figure. 


Figure  59  i*  again  the  heave  axis  criterion,  but  here  the  effects  of  the 


g-suit  upon  perceived  motion  fidelity  are  included. 
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are  as  indicated.  The  same  Kg  values  apply  for  each  break  frequency  as  in 
Fig.  38.  From  this  plot  it  is  obvious  that,  although  the  heave  motion  does 
not  contribute  positively  to  perceived  normal  axis  motion  fidelity,  the 
g-suit  itself  contributes  a great  deal.  This  is  indicated  by  the  small 
JZga  values  for  all  break  frequency  values,  whereas  in  Fig.  38,  all  values 
were  1.0  or  greater.  The  optimum  conditions  for  the  heave  axis  in  terms 
of  attaining  the  best  perceived  normal  axis  motion  fidelity  specify  using 
the  g-suit,  but  not  heave  axis  motion. 

2.  Roll-Sway  Axes 

Optimization  of  the  roll-sway  axes  is  a two-step  process.  First  the 
subcriteria  values  (jy,  Jyy,  Jp)  for  each  scaling  gain-break  frequency 
parameter  set  must  be  computed.  This  parameter  set  involves  two  scaling 
gains,  and  K2R  and  two  break  frequencies,  Bp  and  ^/Cpy  . A third  scaling 
gain,  Ky,  is  unity  in  accordance  with  the  results  of  the  simulation  experi- 
ments . Next  the  subcriteria  are  combined  into  the  single  criterion,  Jyp. 
The  minimum  value  of  Jyp  over  the  parameter  set  indicates  which  parameter 
set  is  optimum. 

Figures  ^0  through  k2  are  plots  of  the  respective  Jp,  Jy  and  Jyy  sub- 
criteria. Each  varies  with  the  two  break  frequency  parameters  Bp  and  C2y. 
The  Kp  scaling  gains  required  to  satisfy  the  simulator  roll  constraint  of 
±25  deg  and  sway  constraint  of  ±10  ft  are  listed  in  Fig.  40  as  a function 
of  the  break  frequency  parameters.  These  Kp  values  are  calculated  assuming 
a constant  ratio  K2r/%  =2.0  for  the  sake  of  simplicity.  These  values  of 
the  scaling  gain  are  computed  using: 

KPl  PDIS  (roll)  = 25  deg 

Kp2  PDIS  (sway)  = 10  ft 

Kp  = min  (Kp.j,  Kpp) 


where  the  H)IS  data  are  from  Table  C-11  through  C-13  wherein  K^/Kp  = 2.0. 

The  three  criteria,  Jp,  Jy  and  Jyy  are  combined  to  form  Jyp,  which  is 
plotted  versus  Bp  and  C2y  in  Fig.  ^3.  Now  the  optimum  drive  logic  parameters 
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Figure  4l . LAMARS  Optimization  — Jy  vs.  C^y  and  Bp 


parameters  can  be  identified  simply  by  choosing  the  set  of  parameters  which 
results  in  the  smallest  value  of  Jyp.  This  point  is  identified  on  the  plot 
as  the  point  which  corresponds  to  % = 1 .0,  Kp  = 0.06,  KgR  = 0.12,  Cgy  = 0.0025 
and  Bp  = 0.5.  These  are  then  the  optimum  roll-sway  drive  logic  parameters  for 
the  assumed  form  of  the  LAMARS  drive  logic  and  hardware  constraints,  and  for 
the  air-to-ground  scenario. 

This  optimization  does  not  include  optimization  with  respect  to  the 
K2r/%>  ra*io.  Simple  expansion  of  the  parametric  study  to  include  this 
effect,  if  desired,  is  possible. 

3.  Pitch-Surge  Axes 

Optimization  of  the  pitch- surge  axes  in  principle  at  least  is  similar 
to  that  for  the  roll- sway  axes.  However,  the  fact  that  there  is  no  surge 
degree  of  freedom  for  the  LAMARS  motion  simulator  requires  the  Kx  scaling 
gain  to  be  zero  and,  in  turn,  makes  the  C2x  break  frequency  parameter  super- 
fluous. Consequently,  optimization  of  the  pitch- surge  axis  drive  logic  is 
considerably  specialized  by  comparison  with  the  roll- sway  axis  procedure  for 
LAMARS.  Two  scaling  gains,  Kq  and  Kjip,  and  one  break  frequency,  Bq  are 
involved  in  optimization.  Since  only  three  parameters  are  involved,  it  is 
not  necessary  to  form  the  ratio  K2p/Kq}  Kgp  effects  can  be  included  directly. 

Figures  44  and  ^5  are  the  plots  of  Jq  and  Jy j,  respectively.  Since 
Kx  = 0.0,  Jy  is  always  unity.  The  values  computed  for  Kq  assume  a derotated 
pitch  authority  of  4.76  deg.  This  authority  is  computed  on  the  assumption 
that  the  heave  axis  would  use  its  full  authority*.  Since  the  optimum  Jz 


♦Calculation  of  the  approximate  pitch  and  yaw  authorities  remaining  after 
derotation  (required  to  accomnodate  LAMARS  simulator  cab  translation)  is  made 
as  follows: 

&(angular  motion)  > g - sin"’  «« 

= 25  - Sin-1  10/26  =■  2.38  deg 

The  available  authority  is  obtained  by  applying  a factor  of  2.0  to  the  2.38  deg 
in  order  to  account  for  the  correlation  between  heaving  said  pitching  motions, 
tod  between  swaying  and  yawing  motions.  The  factor  2.0  is  an  educated  guess 
as  to  an  appropriate  value,  and  may  require  revision  in  the  future.  The  result- 
ing estimate  of  available  authority  available  for  reproducing  angular  motions 
is  4.76  deg. 


indicated  that  the  heave  axis  should  be  stationary,  a full  25  deg  of  pitch 
authority  is  actually  available.  Thus,  the  values  for  Kp  shown  in  Fig.  44 
could  be  multiplied  by  a factor  of  approximately  5 when  the  heave  axis 
is  stationary. 

Figure  46  is  a plot  of  the  combined  criterion,  Jxq,  computed  using  the 
subcriteria  Jq,  Jx  and  Jxt*  The  minimum  value  of  JxQ  identifies  the  follow- 
ing optimum  pitch- surge  drive  logic  parameter  set:  Kgp  = 1.0,  Kx  = 0.0, 

Kq  = 0.1,  Bq  = 0.125  (rad/sec).  If  the  full  pitch  displacement  limits  are 
used  (in  the  absence  of  heave  motion),  the  effect  upon  the  Jq  component  of 
Jxq  is  considerable.  For  example,  if  K2p  =1.0  (which  minimizes  Jxt)#  then 
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The  optimum  set  of  parameters  is  then: 

K2P  = 1.0,  Kjj  = 0.0,  Kq  = 0.546,  Bq  = 0.0675  rad/sec 

corresponding  to  Jxq  = O.892.  Notice  that  the  decrease  in  Jxq  (from  0.950) 
with  the  increase  in  the  pitch  limits  from  4.76  deg  to  25  deg  is  significant. 
(Notice  also  that  absence  of  the  surge  degree  of  freedom  which  places  a lower 
bound  on  Jxq  of  0.550.) 

4.  Yaw  Axis 

Optimization  of  the  yaw  axis  proceeds  similarly  to  that  for  the  heave 
axis  without  g-suit  cues.  Figure  47,  a plot  of  Jr  versus  Bp,  illustrates 
the  ineffectiveness  (Jr  i 1.0)  of  yaw  axis  motion  for  other  than  very  low 
values  of  By.  Therefore,  the  optimized  parameters  for  the  yaw  axis  are  that 
pair  for  which  the  computed  value  of  Jr  is  least. 

Two  comments  are  in  order.  First,  the  optimum  parameters  for  the  yaw 
axis  drive  logic  may  be  affected  by  the  presence  of  dutch  roll  motions 
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Figure  46.  LAMARS  Optimization  — JXq 


LAMARS  Optimization  — Jr  vs 


which  are  not  represented  in  the  simple  math  model  used  to  generate  the 
air-to-ground  scenario  data  used  here.  This  point  may  be  gotten  around 
by  using  data  generated,  say,  in  a manned  simulation  of  the  task  or  in- 
flight. Second,  Fig.  47  implies  that  even  lower  values  of  Bq  may  produce 
lesser  values  of  Jr  than  does  the  lowest  value  of  Br  used  in  calculating 
Jr  here.  As  Br  ultimately  approaches  zero  it  can  be  shown  that  Ky  for  this 
air-to-ground  scenario  approaches  25/360  = 0.0694  and  Jr  therefore  approaches 
0.866.  Consequently,  further  reductions  in  the  yaw  axis  washout  break  fre- 
quency, Bj»  from  0.0675  rad/sec  will  result  in  only  slight  reduction  in  Jr 
from  0.900  while  the  reduction  in  Kr  is  approximately  by  a factor  of  2.5. 
Therefore,  Bp  = 0.0675  rad/sec  is  regarded  as  optimum. 

Table  28  presents  a summary  of  the  optimum  drive  logic  parameters  as 
determined  by  identifying  those  parameter  combinations  for  which  the  criterion 
value  for  each  axis  or  set  of  axes  was  a minimum. 

D.  NONLINEAR  WASHOUT  INVESTIGATION 

Once  the  linear  drive  logic  parameters  were  optimized,  a limited  investi- 
gation of  two  nonlinear  drive  logic  schemes  was  undertaken.  These  are  the 
Parrish  adaptive  gain  scheme,  described  in  Section  III  and  Refs.  7-9,  and  the 
subliminal  scheme  described  in  Section  III,  Appendix  B and  Refs.  10  and  11. 
Performance  is  compared  (in  terms  of  the  motion  fidelity  criterion)  for  one 
set  of  parameters  for  the  Parrish  configuration,  two  sets  of  parameters  for 
the  subliminal  scheme,  and  the  optimized  linear  scheme.  The  parameters  for 
each  washout  are  summarized  in  Table  29*  Note  that  parameters  are  chosen 
for  both  the  Parrish  and  the  subliminal  washouts  to  make  full  use  of  the 
available  authority  limits  of  the  LAMARS  motion  base.  Tables  of  the  motion 
fidelity  measure  data  sure  included  in  Appendix  C,  Table  C-18.  It  is  interest- 
ing to  observe  the  large  increases  in  Kp  or  K2R  which  are  afforded  by  the 
subliminal  schemes,  compared  to  the  linear  scheme. 

Several  other  sets  of  parameters  for  both  the  Parrish  scheme  and  the 
subliminal  scheme  were  investigated,  but  for  brevity  only  the  three  cases 
producing  the  smallest  criterion  values  are  presented  here. 

Figure  40  presents  a plot  of  Jp,  JY  and  Jyr  for  each  of  the  washout 
schemes.  These  subcriteria  are  the  constituents  of  Jyp,  and  it  is 


TABIE  26.  OPTIMUM  LAMARS  DRIVE 

LOGIC  PARAMETERS 

Pitch-Surge  Axes 

Kq 

0.1  » 

* 

Kgp 

1.0  J 

(*  1.*  Bq) 

0.175 

rad/sec 

«B5P  (■  Bq) 

0.125 

rad/sec 

to^p  (=  Bq) 

0.125 

rad/sec 

?3P 

0.7 

Roll-Sway  Axes 

*P 

0.06 

*2R 

0.1 

% 

1.0 

> 

-3- 

• 

II 

0.7 

rad/sec 

“>1R 

0.0 

rad/sec 

“2R 

0.0 

rad/sec 

“^R  (=  ®p) 

0.5 

rad/sec 

“*R  (-  Bp) 

0.5 

rad/sec 

Cly/(2^) 

0.7 

0.0025 

(rad/sec)2 

?3R 

0.7 

Yaw-Heave  Axes 

■ 

Kr 

0.176 

^‘*22 

0.0 

air  (-  Br) 

0.0675 

rad/sec 

C1Z 

0.0 

rad/sec 

C2Z 

0.0 

(rad/sec)2 

C3z  (-  Bg) 

0.0 

rad/sec 

cuz  (-  b|) 

0.0 

(rad/sec)2 

C5Z  Bg) 

0.0 

rad/sec 

*Gains  confuted  based  on  a 
If  the  heave  axis  is  not  used, 
the  full  23  deg,  and  the  pitch 
approximately  3. 

derotated  pitch  authority  of  4.76  deg. 
as  indicated,  then  the  pitch  authority  is 
gain,  Kq,  can  be  increased  by  a factor  of 

TABLE  29.  SUMMARY  OP  PARAMETERS  USED  IN  THE 
NONLINEAR  WASHOUT  INVESTIGATION 


/ 


LINEAR 

PARRISH 

SUBLIMINAL  1 

SUBLIMINAL 

*p 

0.06 

(Adaptive) 

0.12 

0.05 

0.1 

0.1 

0.2 

0.5 

*y 

1.0 

1.0 

1.0 

1.0 

"p 

0.7 

0.7 

0.7 

0.7 

“>1R 

0.0 

0.0 

0.0 

0.0 

<*>28. 

0.0 

0.0 

0.0 

0.0 

“5R 

0.5 

0.5 

0.5 

0.5 

“4R 

0.5 

0.5 

0.5 

0.5 

c2y 

0.0025 

0.0025 

0.0025 

0.0025 

*Pmin 

0.06 

Cq, 

1.0* 

cK 

0.1* 

C2y’ 

0.25 

0.25 

* These  parameters  vere  chosen  from  the  investigation  reported  in 
Ref.  U. 
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Wubout  Comparison 
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Interesting  to  compare  the  various  J values  among  the  washouts.  Figure  U9 
is  the  plot  of  Jyp  for  each  washout.  There  is  no  significant  difference 
among  the  schemes  on  the  basis  of  the  motion  fidelity  criterion  value  alone, 
although  the  linear  washout  has  a slightly  smaller  Jyp  value  than  any  of 
the  nonlinear  schemes.  Within  this  limited  context  (which  does  not  include 
subjective  pilot  evaluation),  it  would  appear  that  the  linear  washouts  used 
in  the  LAMARS  roll-away  drive  logic  cannot  be  improved  upon  by  either  of 
the  two  nonlinear  configurations  examined. 


Figure  49.  Washout  Comparison 
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SECTION  VII 

RESULTS  AN D CONCLUSIONS 


The  results  of  this  research  can  be  summarized  as  follows: 

1 . A method  for  abstracting  mission  scenarios  from  strip  chart  data 
or  verbal  descriptions  was  developed  and  utilized.  This  method,  imple- 
mented in  FORTRAN,  uses  approximate  time  histories  for  p,  q and  Vm  and 
linearized  kinematic  equations  to  compute  body  axis  angular  velocity 
and  specific  force  time  histories  which  approximate  actual  flight  val- 
ues for  a particular  scenario.  The  latter  quantities  provide  realistic 
and  relatively  accurate  inputs  for  evaluating  and  optimizing  motion 
base  drive  logic. 

2.  A new  nonlinear  washout  technique  (the  subliminal  washout  scheme) 
was  developed.  This  scheme  exploits  the  indifference  threshold  which 
is  generally  acknowledged  to  be  present  in  the  specific  force  sensory 
paths.  In  effect,  the  nonlinear  washout  scheme  operates  to  increase 
the  washout  rate  when  the  input  specific  force  is  subthreshold  and 
the  cab  is  in  an  off-center  position.  This  is  in  order  to  drive  the 
cab  to  its  center  position  at  faster  rate  than  it  would  be  for  a 
comparable  linear  scheme.  For  one  input,  a roll-in  to  a constant 

h g turn,  a 70  percent  reduction  in  the  peak  lateral  travel  of  the 
cab  was  demonstrated- for  the  subliminal  scheme  relative  to  a compar- 
able linear  scheme. 

5.  Criteria  for  discerning  fidelity  of  motion  reproduction  were 
developer..  These  criteria  are  composed  of  various  mean- square  error 
measures  which  reflect  differences  between  simulated  and  actual  in- 
flight motion.  The  criteria  can  be  used  to  answer  three  questions 
which  arise  for  every  moving  base  simulation  program: 

• Is  moving  base  simulation  appropriate  for  the  flying 
tasks  to  be  investigated? 

• How  can  motion  base  drive  logic  parameters  be  opti- 
mized for  the  flying  tasks  to  be  investigated? 

• Will  this  set  of  motion  drive  logic  parameters  result 
in  motion  within  simulator  capabilities  for  the  fly- 
ing tasks  to  be  investigated? 

The  impact  of  the  three  sets  of  variables  which  quantify  the  given 
elements  of  the  above  questions  can  be  investigated  individually  or 
simultaneously  via  these  motion  fidelity  measures  and  criteria:  flying 
task  variables,  motion  base  hardware  displacement,  velocity  and  accelera- 
tion constraints,  and  parameters  of  the  motion  base  drive  logic. 


-t. 


155 


Obviously,  with  such  a large  number  of  variables,  it  was  necessary 
to  limit  the  scope  of  the  parametric  study  of  this  research  investi- 
gation to  those  combinations  required  to  demonstrate  formalized,  orderly 
methods  for  identifying  answers  to  the  three  questions  posed  above. 

k.  The  utility  of  the  criteria  and  measures  was  demonstrated  in  an 
application  to  optimize  the  LAMARS  drive  logic  with  respect  to  its  hard- 
ware constraints  and  drive  logic  parameters.  An  air-to-ground  flying 
task  was  used.  The  optimized  parameters  were  identified  for  the  yaw, 
roll- sway  and  pitch- surge  axes.  Furthermore,  the  criterion  value  for 
the  heave  axis  indicated  that  motion  in  that  axis  does  not  add  to  motion 
fidelity.  The  heave  axis  criterion  did  indicate  the  effectiveness  of 
the  g-suit  for  presenting  normal  specific  force  cues  is  significant. 

In  addition,  a comparison  of  the  roll-sway  criterion  values  for  the 
optimized  linear  drive  logic  and  two  nonlinear  schemes  showed  that 
in  this  case,  no  particular  benefit  was  derived  from  using  the  non- 
linear techniques. 

It  is  interesting  to  compare  the  current  LAMARS  drive  logic  para- 
meters for  the  F-17  configuration  (Ref.  l)  with  the  optimized  drive 
logic  parameters  based  upon  the  motion  fidelity  criteria.  Such  a com- 
parison is  made  in  Table  30.  The  major  difference  between  the  two 
parameter  sets  is  in  the  heave  axis  — the  fidelity  criterion  for 
that  axis  indicated  that  motion  did  not  improve  the  simulation  fidel- 
ity. Also,  the  parameters  for  the  pitch-surge  axis  set  are  somewhat 
smaller  than  the  current  LAMARS  values.  This  is  due  to  the  very  low 
(4.76  deg)  pitch  authority  which  is  assumed  to  remain  after  satisfy- 
ing cab  derotation  requirements.  As  was  pointed  out  in  an  earlier 
section,  the  pitch  authority  in  the  absence  of  heave  axis  motion  is 
the  full  25  deg.  Thus,  the  optimized  values  correspond  to  this  full 
authority  since  heave  axis  motion  is  not  used  for  optimized  drive 
logic. 

In  conclusion,  the  main  goal  of  the  research  reported  herein  was 
to  establish  an  orderly,  relatively  simple  method  for  optimizing  the 
presentation  of  motion  cues  in  moving  base  simulations.  This  is 
accomplished  via  choice  of  the  drive  logic  parameters  for  a given 
simulator  and  flying  task.  The  method  developed  here  is  based  upon 
use  of  motion  fidelity  criteria.  It  has  been  applied  to  optimize 
the  LAMARS  motion  base  drive  logic  parameters  for  an  air-to-ground 
scenario.  The  software  supporting  the  optimization  process  is  docu- 
mented in  a separate  volume.  The  methods  and  procedures  of  this 
report  can  be  applied  for  any  given  motion  base  and  any  desired 
scenario.  Preliminary  results  with  this  optimization  procedure 
indicate  that  further  applications  and  research  in  this  area  are 
warranted.  Suggestions  for  directions  which  future  research  should 
take  are  found  in  the  next  and  final  section. 
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DRIVE  LOGIC  PARAMETERS 
PARAMETERS 


Optimized 
Linear  LAMARS 


0.55  (0.1)* 

1.0 

0.09U5  (0.175)* 

rad/sec 

0.0675  (0.125)* 

rad/sec 

0.0675  (0.125)* 

rad/sec 

0.7 

0.0 

rad/sec 

0.06 

0.1 

1.0 

0.7 

rad/sec 

0.0 

rad/sec 

0.0 

rad/sec 

0.0 

rad/sec 

0.5 

rad/sec 

0.5 

rad/sec 

0.7 

0.0025 

(rad/sec)2 

0.7 

0.176 

0.0 

0.0675 

rad/sec 

0.0 

rad/sec 

0.0 

rad/sec 

0.0 

(rad/sec)2 

0.0 

rad/sec 

0.0 

(rad/sec)2 

0.0 

rad/sec 

axe  reduced  to  accommodate  derot a- 
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SECTION  VIII 


SUGGESTIONS  FOR  FURTHER  RESEARCH 


The  procedure  developed  for  generating  scenario  data  and  optimizing 
motion  simulator  drive  logic  is  a research  tool  at  this  stage.  Figure  50 
provides  an  overview  of  this  tool,  outlining  the  major  steps  of  the  procedure. 
Steps  in  blocks  with  light  outlines  are  manual  engineering  procedures,  while 
steps  in  blocks  with  heavy  outlines  are  implemented  in  software.  Further 
development  of  this  tool  should  convert  it  to  a form  suitable  for  routine 
use  as  an  engineering  simulator  experiment  design  aid.  Three  general  areas 
deserve  further  attention: 

1 . Additional  applications  of  the  optimization  procedure  should 
be  undertaken  to  shed  light  on  its  strengths  and  weaknesses  and 
suggest  areas  which  might  be  changed  and  improved.  The  greatest 
benefit  would  be  derived  if  this  tool  was  used  in  a variety  of 
different  simulation  environments  (different  motion  bases,  dif- 
ferent scenarios  and  different  drive  logic  schemes). 

2.  Motion  fidelity  criteria  should  be  sharpened.  At  present 
the  four  motion  fidelity  criteria  (one  for  each  of  roll- sway, 
pitch-surge,  yaw  and  heave  motion  axes)  are  each  functions  of 
several  parameters  (a,  0,  X).  Values  for  these  parameters  have 
been  inferred  from  a rather  small  amount  of  experimental  evi- 
dence. Additional  critical  experiments  are  recommended  to 
strengthen  the  data  base  upon  which  the  a,  0 and  X parameter 
values  are  based,  and  perhaps  to  provide  revised  values  for 
those  parameters. 

3.  Further  applications  of  this  tool  would  benefit  from  in- 
creased automation  of  the  optimization  procedure.  This  auto- 
mation would  result  in  less  cumbersome  "trial  and  error"  choice 
of  drive  logic  parameters,  for  example,  if  the  motion  fidelity 
criteria  optimizations  were  performed  via  a computerized  search. 

With  these  three  areas  targeted,  the  following  immediate  research  plan 
is  suggested: 

Task  1 . Extend  Capabilities  of  Software 

1 . 1 Convert  existing  software  from  its  present  research  con- 
figuration to  an  engineering  simulator  experiment  design 
configuration. 


1.2  Add  "motion  fidelity  criteria  evaluation"  function  to 
existing  software. 

1.3  Add  automatic  "optimization  of  drive  logic  parameters" 
(parameter  search)  function  to  software. 

Task  2.  Expand  Flying  Task  Data  Base 

2.1  Construct  engineering  abstractions  for  additional  air- 
to-ground  and  air-to-air  flying  tasks. 


2.2  Revise  existing  software  to  facilitate  input  of  flying 
task  data  from  piloted  simulation  runs. 
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APPENDIX  A 


^ § 


inputs  resulting  from  an  engineering 
abstraction  of  the  flying  task 


vT  * 4vt 


(longitudinal  kinematic  acceleration) 


KINEMATIC  EQUATIONS  FOR  SCENARIO  GENERATION 

The  overall  description  of  the  method  used  to  generate  scenarios  for 
this  research  was  presented  in  a previous  section.  Hie  purpose  of  this 
Appendix  is  to  present  approximated  kinematic  equations,  derived  from  the 
exact  equations  of  motion.  These  are  used  to  compute  the  drive  logic  inputs 
corresponding  to  various  flying  tasks.  The  summary  of  assumptions  needed 
to  arrive  at  the  kinematic  equations  from  the  equations  of  motion  is  given 
below: 

• For  a given  small  angle  e,  sin  e = e,  tan  e = e,  cos  e = 1 

• $Q  and  e are  small  angles 

• By  and  0 are  identically  zero 

• a is  small 

where  e is  pitch  attitude,  a is  angle  of  attack  and  50  is  the  thrust  incli- 
nation angle. 

The  ten  basic  kinematic  equations  are  given  below.  The  additions  or 
deletions  to  this  list,  which  are  scenario-dependent,  are  discussed  in  the 
succeeding  paragraphs. 


CP  = P 


(bank  angle  rate/roll  rate) 


♦ = r cos  cp  + q sin  cp  (heading  rate) 


(yaw  rate;  from  requiring 


(longitudinal  acceleration) 


The  first  seven  equations  are  integrated,  at  each  time  interval,  resulting 
in  ten  variables.  Except  as  indicated  below,  all  initial  conditions  are 


For  Scenarios  1 and  2 (the  first  two  HQDT  scenarios)  q is  not  an  input 
Instead,  q is  computed  from  an  equation  which,  if  satisfied,  assures  that 
the  airplane  is  in  constant  altitude  flight. 


Maintenance  of  constant  altitude  while  in  a coordinated  turn  requires 
the  vertical  conponent  of  lift  to  equal  the  weight.  Let  load  factor,  Nz, 
be  defined  by: 


Solving  for  (l  +NZ)  and  differentiating  both  sides  of  the  resulting  equation 
gives 

Nz  = p sin  cp/cos2  cp 


where  p is  in  rad/sec  units.  Normal  acceleration,  load  factor,  and  pitch 

• 9 • 

rate  are  related  in  the  following  equation  wherein  w = aV^  + aV^  has  been 
assumed  to  be  zero. 


1 


Solving  for  q In  terms  of  Nz,  and  differentiating  (assuming  V<j>  is  constant) 
results  in; 


Substitution  for  Nz  results  in 


In  addition,  the  following  initial  condition  is  imposed 


For  Scenario  3 (the  HQDT4  scenario),  p and  q are  not  always  input  as 
discrete  values.  During  certain  intervals  the  roll  and  pitch  rates  required 
to  achieve  a given  load  factor  onset  rate  while  maintaining  constant  alti- 
tude flight,  must  be  computed  on  the  basis  of  the  appropriate  constraint 
equations.  The  following  equations  are  used  to  compute  inputs  to  the  kine- 
matic equations  during  these  intervals. 


If  the  previous  equation  for  Nz  is  solved  for  p said  then  differentiated 
assuming  Nz  is  constant  the  result  is: 


If  the  assumption  that  V?  is  constant  is  relaxed  in  obtaining  the  q 
equation,  the  result  is: 


q = (g/V^)(Nz  + p sin  9 - (Vt/Vit)(1  + Nz  - cos  9)]  rad/sec 


As  was  true  for  the  other,  HQDT  scenarios,  the  following  initial  condition 
is  imposed: 

VT(0)  = 560  ft/sec 

Scenario  4,  the  air-to-ground  task,  imposes  the  following  initial  condi- 
tions : 

VT(0)  = 420  ft/sec 
❖T(°)  = -5  ft/ sec2 

In  order  to  simulate  the  air-to-air  tasks  (Scenarios  5-8),  it  was  neces- 
sary to  remove  the  assumption  that  a,  the  angle  of  attack  rate  of  change,  ia 
small.  A Act  equation  was  added,  and  several  equations  were  changed  to  include 
the  effect  of  a.  In  addition,  it  was  necessary  to  make  one  new  assumption. 
This  assumption  and  the  resulting  changes  in  the  kinematic  equations  which 
are  affected,  sure  given  below: 

a = Oq  + Aa 

where  a©  is  the  trim  angle  of  attack  for  straight,  level,  upright,  constant 
speed  flight. 

The  kinematic  terms 


-pw/VT  = “P(Aa  + Qo) 

wq  - Vt  sin  (a©  + Aa)Ax  = VrCa©  + Aa)(q  - Aa,) 


w = AxVrji  + Vt(oo  + Aa) 


are  restored  to  the  Sy/V^,  a*  and  az  equations,  respectively. 

r = (g/Vj)  sin  cp  + p(Aa  + a©) 
ax  = Vj  + g0  + Vt(oo  + Ax)(q  - Ax) 
az  = -V^q  - g cos  cp  + AxV  + Vt(oo  + Aa) 


The  normal  acceleration  Is  alternatively  given  by 


/ ' ' 


az  = — g + Z"  V<j^x(Vi/Vt(0)] 


where  z*  represents  the  primary  aerodynamic  characteristic  of  the  aircraft 
via 


pSVT(°) 

Zw  = —an—  (°La  + Cp) 


If  az  is  eliminated  from  the  two  expressions  for  az  above,  and  the  result- 
ing expression  is  solved  for  Ao,  the  result  is 


g(  1 — cos  cp)  (— 2w)AoV t VT(oto  + Aff.) 

**  = * - --  - y2--v^o) 


vT 


VT 


where  Vt(<x0  + £&)  /Vt  is  usually  negligible  with  respect  to  the  other  terms 
of  the  equation. 


Also,  the  following  initial  conditions  or  constants  apply: 


°o  = 

4 deg 

0(0)  = 

4 deg 

vT(o)  = 

800  ft /sec 

-Zw  = 

1 

Thus,  for  each  scenario  there  exists  a set  of  ten  (or  eleven,  for  the 
air-to-air  tasks)  state  and  algebraic  equations,  and  the  assumption  that 
ay  = 0.  The  drive  logic  inputs  are  generated  from  these.  For  the  angu- 
lar inputs,  this  is  a trivial  matter: 


PA  = P 

<*A  " * 

rA  = r 
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The  specific  force  inputs  are  complicated  by  the  fact  that  the  pilot's 
location  in  the  aircraft  is  significantly  displaced  from  the  aircraft  c.g. 

In  order  to  transform  the  ax,  ay  and  quantities  to  UpQ,  Vpo,  WpQ  specific 
forces  as  sensed  by  the  pilot,  it  is  necessary  to  use  the  following  equa- 
tions. 


In  each  case,  the  following  numbers  were  used  for  the  "location"  constants 
(typical  for  the  A- 10  aircraft) : 


For  details  on  the  implementation  of  these  kinematic  equations  as  com- 
puter software,  refer  to  the  Software  User's  Guide  (Ref.  2). 


APPENDIX  B 


SUBLIMINAL  WASHOUT  SCHEME  INVESTIGATION 

This  appendix  presents  the  nonlinear  washout  scheme  designed  in  the 
course  of  the  LAMARS  optimization  research.  Figure  B-1  depicts  an  appli- 
cation of  this  subliminal  washout  scheme  to  a first-order  roll  axis  washout. 
This  concept  is  the  result  of  an  attempt  to  utilize  so-called  "indifference" 
thresholds  which  pilots  exhibit  under  normal  workload.  These  thresholds 
may  be  operative  for  both  angular  velocity  and  specific  force  perception. 

The  hypothesis  is  that  pilots  do  not  perceive  angular  velocities  and  specific 
forces  which  have  magnitudes  below  the  respective  indifference  thresholds. 

It  was  hoped  that  exploitation  of  this  particular  phenomenon  would  yield 
a design  which  reduced  simulator  travel  requirements. 

The  overall  design  goal  is  to  drive  the  simulator  cab  towards  its  zero 
position  more  rapidly  than  would  the  underlying  linear  washout  whenever 


p * Input  Angulor  Velocity 
ps  = Scaled  Angular  Velocity 
pT  - Indifference  Threshold  Level 
^ = Commanded  Cob  Roll  Rote 
+ s Commanded  Cab  Roll  Angle 


Figure  B-1.  Subliminal  Washout  Filter, 
First-Order  Case 


the  motion  stimulus  is  below  the  indifference  threshold  level.  This  is 
accomplished  by  using  the  two  nonlinear  functions  in  boxes  A and  B of 
Fig.  B-1.  The  input  to  the  function  in  A is  the  scaled  angular  velocity. 
This  function  produces  a weighting  factor  which  serves  as  a variable  feed- 
back gain  in  the  washout  circuit.  If  the  input  magnitude  is  larger  than 
the  indifference  threshold,  p^,  the  weighting  factor  is  zero.  If  the  input 
is  less  than  p>p,  the  weighting  factor  is  a fraction  less  than  1.0,  which 
is  a sinusoid-like  function  of  the  input  (for  the  form  of  the  weighting 
function  used  here). 

The  input  to  B,  a soft  saturation  nonlinear  function,  is  cab  roll 
angle  cp.  If  cp  is  large,  the  value  of  the  function  output  is  equal,  to 
the  indifference  threshold,  ±pr*  If  cp  is  small,  the  value  of  the  func- 
tion output  is  proportional  to  cp. 

The  outputs  from  boxes  A and  B are  then  multiplied  to  arrive  at  an 
incremental  washout  rate  command  signal.  The  particular  choice  of  func- 
tions assures  that  this  signal's  magnitude  never  exceeds  the  indifference 
threshold  level.  Also,  the  smoothness  of  those  functions  tends  to  prevent 
discontinuous  commanded  changes  in  the  washout  rate.  The  value  of  this 
incremental  washout  rate  command  signal  will  be  non- zero  whenever  the  cab 
roll  angle  is  non- zero  and  the  scaled  angular  velocity  is  below  the  indif- 
ference threshold  level.  This  signal  is  then  subtracted  from  the  scaled 
input  singular  velocity.  The  subtracted  component  provides  an  additional 
toward-center  angular  velocity  input  to  the  integrator  not  present  in  the 
underlying  linear  scheme.  Thus,  the  cab  is  driven  back  to  its  zero  posi- 
tion more  quickly  than  it  would  be  for  the  linear  scheme  during  intervals 
of  subthreshold  angular  velocity  inputs. 

The  assumed  form  of  the  indifference  threshold  is  shown  in  Fig.  B-2a. 
The  soft  saturation  function  of  box  B,  Fig.  B-1,  is  complementary  to  the 
indifference  threshold  function  and  is  shown  in  Fig.  B-2b.  The  nonlinear 
function  in  box  A (which  generates  the  weighting  factor)  approximates  the 
partial  derivative  of  the  soft  saturation  function  with  respect  to  its 
argument.  This  function  is  shown  in  Fig.  B-2c.  The  specific  choice  of 


vV, 


a)  Presumed  Form  of  Indifference  Threshold 


b)  Complementary  Function  For  Indifference  Threshold 


c)  Approximate  ' Partial  Derivative  of  Complementary  Function 


Figure  B-2.  Summary  of  Nonlinear  Functions 


the  indifference  function  is  actually  somewhat  arbitrary,  and  subject  only 
to  the  constraints  that  T(o)  = 0,  A(o)  = 1 and  B(±«)  = ±p^.  The  particular 
analytical  functions  used  here  are  also  indicated  for  each  part  of  Fig.  B-2. 

Preliminary  tests  of  this  subliminal  washout  concept  for  the  roll  axis 
showed  it  to  be  ineffective  for  reducing  simulator  maximum  roll  and  sway 
displacements.  There  was  some  reduction  in  simulator  motion  requirements, 
but  not  really  enough  to  warrant  further  investigation.  Thus,  all  sub- 
sequent study  centered  around  schemes  which  included  the  subliminal  functions 
in  the  specific  force  paths  only.  Figure  B-3  shows  an  application  of  the 
subliminal  washout  concept  within  the  lateral  specific  force  channel  of  a 
crossproduct  washout  configuration  for  the  roll- sway  axes. 

For  crossproduct  schemes  such  as  the  one  in  Fig.  B-3,  the  scaled  input 
specific  force  K2RVp0  and  the  recovered  specific  force,  Vp0p  are  always 
equal  for  Ciy  = Cpy  = C^y  = Cpy  = 0.  This  is  characteristic  of  the  cross- 
product  configuration  (Refs.  1,  4)  for  implementation  of  the  residual  tilt 
and  coordinating  crossfeed  paths.  Introduction  of  the  nonlinear  subliminal 
washout  paths  (C£y  > 0 and/or  C^y  > o)  results  in  intentional  miscoordina- 
tion  of  reproduced  roll- sway  motions.  The  indifference  threshold  phenomenon 
permits  this  deliberate  introduction  of  a spurious  specific  force  cue  due 
to  miscoordination.  As  long  as  the  degree  of  miscoordination  does  not  cause 
a spurious  lateral  specific  force  cue  which  exceeds  the  specific  force  indif- 
ference threshold  level,  the  pilot,  under  normal  workload,  should  not  detect 
the  miscoordination.  On  the  other  hand  this  deliberate  miscoordination 
reduces  the  peak  sway  displacement  of  the  simulator  substantially. 

Computer  simulation  results  for  the  washout  in  Fig.  B-3  are  presented 
in  Fig.  B-4.  The  input  to  the  simulation  corresponds  to  a roll-in  to  a 
constant  4 g turn.  The  inputs  sire  roll  rate,  p^,  and  lateral  specific  force, 
VpQ.  There  is  no  significant  reduction  in  the  peak  magnitudes  of  the  sway 
acceleration,  y,  or  the  sway  velocity,  y.  But  these  results  clearly  show 
a substsmtial  reduction  in  the  maximum  sway  displacement,  y.  The  subliminal 
washout  scheme  accounts  for  a 70  percent  reduction  in  y compared  to  the  linear 
case  (maximum  required  y equals  13  ft  for  linear  washout,  4 ft  for  subliminal 
washout) . 
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Figure  B-1+.  Comparison  of  Outputs  for  Linear  and 
Subliminal  Washout  Schemes 


In  order  to  accomplish  this  substantial  reduction  in  lateral  transla- 
tion, however,  a substantial  increase  in  recovered  specific  force,  Vp0p, 
which  is  spurious  is  observed.  This  is  due  to  the  miscoordination  deliber- 
ately introduced  to  increase  the  washout  rate.  Notice  that  this  component 
of  miscoordination  is  constrained  to  be  approximately  within  the  assumed 
indifference  threshold  level  of  0.1  g.  Under  normal  workload,  the  pilot 
should  not  be  able  to  detect  this  level  of  miscoordination. 
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The  computer  simulation  of  the  subliminal  washout  has  been  exercised 
for  a variety  of  inputs.  Significant  reductions  in  translational  require- 
ments have  been  observed.  On  the  basis  of  these  results,  the  following 
conclusions  are  drawn: 

1 . The  subliminal  washout  concepts,  as  implemented  in  the  trans- 
lational axes  of  the  crossproduct  scheme,  are  effective  in 
reducing  the  maximum  displacement  requirements  of  the  motion 
base. 

2.  The  subliminal  washout  scheme  is  most  effective  for  sub- 
indifference threshold  specific  force  inputs.  The  washout 
effectively  reduces  the  tinder  lying  linear  scheme  when 
inputs  exceed  this  threshold  by  a very  large  amount. 

3.  The  use  of  the  subliminal  washout  scheme  results  in  an 
increase  in  recovered  specific  force  which  is  spurious. 

This  spurious  motion  is  due  to  additional  miscoordination. 

The  nonlinear  implementation  insures  that  this  miscoordina- 
tion component  is  never  greater  than  the  assumed  indiffer- 
ence threshold  level.  Thus,  under  normal  workload  the  pilot 
should  be  unable  to  detect  this  false  cue. 
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APPENDIX  C 
TABLES  OF  MEASURES 


This  appendix  presents  the  raw  measures  data  generated  by  the  measures 
software.  The  following  Tables  are  included: 


Table 

C-1; 

Pitch  and  Surge  axis  measures  versus  scenario 

Table 

C-2: 

Roll  and  Sway  axis  measures  versus  scenario 

Table 

C-3: 

Yaw  and  Heave  axis  measures  versus  scenario 

Table 

C-k: 

Pitch  and  Surge  axis  measures  versus  break  frequency, 

Kjjp  = 0.0 

Table 

C-5: 

Pitch  and  Surge  axis  measures  versus  break  frequency, 
-K2p  =0.5 

Table 

C-6: 

Pitch  and  Surge  axis  measures  versus  break  frequency, 

K2p  = 1.0 

Table 

C-7: 

Roll  and  Sway  axis  measures  versus  break  frequency, 

K2R  = 0.0,  C2y  = 0.0025  and  Kp  = 0.5 

Table 

C-8: 

Roll  and  Sway  axis  measures  versus  break  frequency, 

K2R  = 0.0,  C2y  = 0.0225  and  Kp  = 0.5 

Table 

C-9: 

Roll  and  Sway  axis  measures  versus  break  frequency, 

K2R  = 0.0,  C2y  =0.25  and  Kp  = 0-5 

Table 

C-10: 

Roll  and  Sway  axis  measures  versus  break  frequency, 

K2R  = 0.5#  C2y  = 0.0225  and  Kp  = 0.5 

Table 

C-11 : 

Roll  and  Sway  axis  measures  versus  break  frequency, 
k2R  =1*0,  C2y  = 0. 0025  and  Kp  = 0. 5 

Table 

C-12: 

Roll  and  Sway  axis  measures  versus  break  frequency, 

*2R  = 1*°>  c2y  = 0.0225  and  Kp  = 0.5 

Table 

C-13: 

Roll  and  Sway  axis  measures  versus  break  frequency, 
k2R  = 1*0,  C2y  = 0.25  and  Kp  = 0.5 

Table 

C-l4: 

Yaw  and  Heave  axis  measures  versus  break  frequency 

Table 

C-15: 

Revised  measures,  Pitch-Surge  axes  versus  break  fre- 
quency and  K2p 

Table 

C-16: 

Revised  measures,  Roll-Sway  axes  versus  break  frequency 
and  C2y 

Table 

C-17: 

Revised  measures.  Yaw  and  Heave  axes  versus  break 
frequency 

Table 

C-1 8: 

Revised  measures,  Roll-Sway  axes  versus  washout 
configuration 
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TABLE  C-1.  PITCH  AND  SURGE  MEASURES  VERSUS  SCENARIO 


1 

2 

3 

4 

5 

6 

7 

8 

PDLA 

0.75E+02 

0.68E+02 

0.14E+03 

0.26E+02 

0.26E+02 

0.31E+03 

0.42E+02 

0.70E+02 

PDL1 

0.6263? 

0.69404 

0.96003 

1.08905 

0.91431 

0.38442 

0.18338 

0.35533 

PDL2 

0.03132 

0.02884 

0.07142 

0.05983 

-0.07251 

-0.12939 

-0.02069 

-0.04326 

PDL3 

0.73160 

0.72728 

0.97375 

0.94394 

0.91827 

0.35970 

0.41830 

0.29831 

PDL4 

0.01937 

0.01697 

0.03537 

-0.01619 

-0.09861 

-0.11135 

-0.02692 

0.03979 

PDLS 

2.17330 

3.99232 

0.47547 

1.87816 

1.33772 

59.49042 

3.96663 

11.13155 

PSENA 

0.28E+02 

0.26E+02 

0.19E+02 

0.92E+01 

0.35E+01 

0.20E+03 

0. 19E+02 

0.59E+02 

PSEN1 

0.70044 

0.49247 

1.10686 

1.04947 

0.71742 

0.17017 

0.17841 

0.26482 

PSEN2 

0.03407 

0.02190 

0.08566 

-0.02999 

0.03997 

-0.05393 

0.01669 

-0.01967 

PSEN3 

0.74986 

0.46449 

1.02832 

0.82517 

0.62734 

0.14103 

0.40500 

0.21558 

PSEN4 

0.01939 

0.01509 

0.03842 

-0.09485 

-0.23225 

-0.06569 

0.02438 

0.05784 

PSENS 

1.38857 

3.59124 

0.27638 

1.28337 

0.85097 

49.25410 

1.60659 

10.38558 

PNLA 

0. 15E+02 

0.13E+02 

0.66E+01 

0.32E+01 

0.83E+00 

0. 18E+03 

0. 12E+02 

0.45E+02 

PNL 1 

0.81950 

0.74655 

0.99783 

0.83605 

0.81143 

0.11708 

0.36988 

0.29627 

PNL2 

0.01256 

0.01271 

0.00302 

-0.04638 

-0.06634 

-0.00377 

0.10862 

0.02058 

PNL  3 

0.89725 

0.78820 

0.99880 

0.88522 

0.89440 

0.13425 

0.59638 

0.36651 

PNL4 

0.00762 

0.00766 

0.00151 

-0.07366 

-0.05702 

-0.04359 

0.08899 

0.07073 

PNLS 

0.08484 

0.53403 

0.00040 

0.13301 

0.01454 

40.79810 

0.31808 

4.98215 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

0.77E+01 

0.79E+01 

0.32E+01 

0.10E+02 

0.67E+01 

0.1  1E+03 

0.15E+02 

0.22E+02 

PDIS 

17.99040 

17.65545 

14.30437 

34.02669 

41.84064 

86.60170 

22.90169 

25.21867 

a)  Pitch  Axis  Measures  vs  Scenarios 


1 

2 

3 

4 

5 

6 

7 

8 

PDLA 

0.00044 

0.00049 

0.00304 

0.01721 

0.07357 

0.64214 

0.02893 

0.08249 

PDL1 

0.72949 

0.62285 

0.37423 

0.32217 

0.42617 

1.14220 

1.07048 

1.05091 

PDL2 

27.97088 

14.72705 

6.25722 

0.31440 

-0.01809 

-0.11487 

0.26158 

-0.19794 

PDL3 

0.84603 

0.75452 

0.59972 

0.30103 

0.43858 

1.00174 

1.01873 

1.01604 

PDL4 

4.23309 

2.13022 

1.89706 

-0.07105 

-0.04561 

-0.10365 

0.06017 

-0.13039 

PDLS 

0.00851 

0.00517 

0.00802 

0.02025 

0.04577 

0.14842 

0.00504 

0.00674 

PSENA 

0.00022 

0.00023 

0.00112 

0.00525 

0.02174 

0.38365 

0.01351 

0.03357 

PSEN1 

0.75564 

0.73141 

0.36139 

0.66433 

0.73779 

1.22880 

1.06437 

1.07506 

PSEN2 

28.41140 

14.60515 

5.95116 

0.67220 

-0.01640 

-0.13519 

0.24937 

-0.30478 

PSEN3 

0.86104 

0.81675 

0.58409 

0.56371 

0.65744 

1.06298 

1.01575 

1.02733 

PSEN4 

4.19120 

1.55125 

1 . 80800 

-0.08536 

-0.05411 

-0.09795 

0.06433 

-0.19719 

PSENS 

0.00433 

0.00267 

0.00283 

0.00725 

0.01119 

0.06274 

0.00207 

0.00369 

PNLA 

0.13E-07 

0.16E-07 

0.18E-05 

0.50E-03 

0.76E-02 

0.30E+00 

0.33E-02 

0.18E-01 

PNL1 

0.98982 

0.96470 

0.76247 

0.98460 

0.92633 

1.22936 

1.01016 

1.00033 

PNL2 

0.12E+05 

0.41E+04 

0.13E+02 

0.14E+01- 

-0.74E-01- 

-0.13E+00 

0. 15E+00- 

-0.19E+00 

PNL  3 

0.99481 

0.98184 

0.86988 

0.96338 

0.87962 

1.08587 

1.00505 

1.00015 

PNL4 

73.21522 

20.00813 

1.77258 

-0.10080 

-0.05868 

-0.09272 

0.07185 

-0.10343 

PNLS 

0.15E-03 

0.65E-04 

0.18E-04 

0.83E-03 

0.16E-02 

0.33E-01 

0.33E-04 

0.37E-03 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDIS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

b)  Surge  Axis  Measures  vs.  Scenarios 
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TABLE  C-2.  ROLL  AND  SWAY  MEASURES  VERSUS  SCENARIO 


1 

2 

3 

4 

5 

6 

7 

8 

PDLA 

0.17E+03 

0. 44E+03 

0.78E+02 

0.16E+03 

0.65E+02 

0.49E+04 

0.73E+03 

0.20E+03 

PDL1 

1.38120 

1.42730 

1.25512 

1.31958 

1.41428 

1.32469 

1.39301 

1.03192 

PDL2 

-0.02564 

0.01391 

0.03747 

0.22213 

0.02631 

-0.30410 

0.11461 

0.63262 

PDL3 

1.0694? 

1.13320 

0.96080 

0.98028 

1.04476 

0.97984 

1.00646 

0.78101 

PDL4 

0.03563 

0.03391 

0.10162 

0.17313 

0.09063 

-0.12203 

0.14208 

0.33483 

PDLS 

0.06199 

0.11664 

0.03271 

0.0947? 

0.02759 

3.73685 

0.38277 

0.2203? 

PSENA 

0.12E+03 

0.32E+03 

0.56E+02 

0. 14E+03 

0.53E+02 

0.38E+04 

0.61E+03 

0. 16E+03 

PSEN1 

1.50638 

1.49302 

1.31497 

1.33852 

1.44843 

1.20769 

1.37130 

1.02612 

PSEN2 

-0.01938 

0.03665 

0.05242 

0.24990 

0.03628 

-0.26908 

0.09437 

0.75415 

PSEN3 

1.08929 

1.14727 

0.93887 

0.96518 

1.03784 

0.90262 

0.99503 

0.72697 

PSEN4 

0.05724 

0.05351 

0.13570 

0.19767 

0.10771 

-0.10614 

0.13567 

0.40287 

PSENS 

0.05679 

0.10299 

0.03057 

0.08972 

0.02566 

3.31999 

0.31328 

0.20456 

PNLA 

0.87E+02 

0.27E+03 

0.42E+02 

0. 11E+03 

0. 35E+02 

0.37E+04 

0.55E+03 

0. 13E+03 

PNL1 

0.99951 

1.00302 

0.99742 

0.99770 

0.99939 

1.06521 

1.01440 

0.96801 

PNL2 

0.00161 

0.00029 

0.00388 

0.00878 

0.00174 

-0.20411 

0.00247 

0.04672 

PNL3 

0.99973 

1.00149 

0.99853 

0.99832 

0.99966 

0.95896 

1.00470 

0.98307 

PNL4 

0.00083 

0.00016 

0.00210 

0.00485 

0.00090 

-0.09478 

0.00334 

0.023S4 

PN  LS 

0.17E-05 

0.14E-04 

0.31E-05 

0.40E-04 

0.50E-06 

0.12E+01 

0.14E-02 

0.48E-03 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

0. 11E+02 

0. 12E+02 

0.13E+02 

0.25E+02 

0.10E+02 

0. 16E+03 

0.33E+02 

0.31E+02 

PDIS 

0.1 3E+02 

0.25E+02 

0.12E+02 

0.31E+02 

0. 12E+02 

0.11E+03 

0.21E+02 

0.30E+02 

a)  Roll  Axis  Measures  vs.  Scenarios 


1 

2 

3 

4 

5 

6 

7 

e 

PDLA 

0.00027 

0.00098 

0.00027 

0 . 00086 

0.00021 

0.21070 

0.00241 

0.00251 

PDL1 

0.00623 

0.00947 

0.00471 

0.00750 

0.00725 

0.00122 

0.00655 

0.00414 

PDL2 

0.04376 

0.05186 

0.01986 

0.02386 

0.02664 

-0.00034 

0.01716 

0.00901 

PDL3 

0.03755 

0.06113 

0.02539 

0.04468 

0.04324 

0.00669 

0.05267 

0.02167 

PDL4 

-0.01691 

0.05226 

-0.01924 

0.00734 

-0.02448 

-0.00848 

-0.04058 

0.00157 

PDLS 

0.11E-04 

0.33E-04 

0.11E-04 

0.32E-04 

0.85E-05 

0.85E-02 

0.96E-04 

0.97E-04 

PSENA 

0.00011 

0.00050 

0.00011 

0.00044 

0.00011 

0.04460 

0.00126 

0.00100 

PSEN1 

0.01017 

0.01337 

0.00724 

0.01011 

0.01004 

0.00243 

0.00774 

0.00630 

PSEN2 

0.07037 

0.06933 

0.03493 

0.03465 

0.03950 

-0.00058 

0.01702 

0.01605 

PSEN3 

0.06609 

0.09010 

0.04190 

0.06247 

0.06255 

0.01693 

0.0641? 

0.03645 

PSEN4 

-0.0256? 

0.07355 

-0.02216 

0.02013 

-0.02797 

-0.02370 

-0.02866 

0.01462 

PSENS 

0.45E-05 

0.15E-04 

0.45E-05 

0.15E-04 

0.41E-05 

0.18E-02 

0.48E-04 

0.37E-04 

PNLA 

0.11E-07 

0.11E-05 

0.44E-07 

0.28E-05 

0.14E-07 

0.30E-01 

0.93E-05 

0.45E-04 

PNL1 

0.92986 

0.93376 

0.92958 

0.93129 

0.92980 

0.36245 

0.93277 

0.92194 

PNL2 

-0.00361 

0.01182 

-0.00086 

0.00836 

-0.00796 

-0.01272 

-0.00672 

0.00253 

PNL3 

0.96429 

0.96631 

0.96414 

0.96503 

0.96426 

0.57275 

0.96580 

0.96018 

PNL4 

-0.00191 

0.00609 

-0.00047 

P.00431 

-0.00417 

-0.01229 

-0.00350 

0.00132 

PNLS 

0.63E-12 

0.33E-10 

0.24E-11 

0.11E-09 

0.90E-12 

0.28E-03 

0.54E-09 

0.27E-08 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.11E+02 

0.21E+02 

0.12E+02 

0.26E+02 

0.11E+02 

0.42E+03 

0.26E+02 

0.36E+02 

PVEL 

0.11E+02 

0.26E+02 

0.13E+02 

0.34E+02 

0.11E+02 

0.12E+03 

0.31E+02 

0.25E+02 

PDIS 

29.63115 

71.74157 

27.69105 

70.45265 

26.21994 

59.35406 

37.14861 

52.06571 

b)  Sway  Axis  Measures  vs.  Scenarios 


TABLE  C-3.  YAW  AND  HEAVE  MEASURES  VERSUS  SCENARIO 


1 

2 

3 

4 

5 

6 

7 

8 

PDLA 

0.60E+01 

0.61E+01 

0.80E+01 

0.53E+01 

0.57E+01 

0. AAE+03 

0.15E+02 

0. 15E+02 

PDL1 

0.79145 

0.78756 

0.91998 

0.95179 

0.91998 

0.25066 

0.A7009 

0.67853 

PDL2 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL3 

0.87577 

0.83831 

0.95271 

0.95178 

0.92563 

0.25663 

0.56520 

0.69159 

PDLA 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

PDLS 

0.05936 

0.16875 

0.02919 

0.06390 

0.09727 

81.98996 

1.23166 

1.10752 

PSENA 

0.19E+01 

0.42E+01 

0. 10E+01 

0. 16E+01 

0.11E+01 

0.A0E+03 

0. 12E+02 

0.76E+01 

PSEN1 

0.8526A 

0.80533 

0.92171 

0.87670 

0.67222 

0.2A966 

0.A6983 

0.A6171 

PSEN2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN3 

0.88767 

0.8A598 

0.92257 

0.87832 

0.67258 

0.25636 

0.577A6 

0.A6106 

PSENA 

0 • 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSENS 

0.036A5 

0.11785 

0.0195A 

0.0A893 

0.08829 

7A. 19031 

0.96381 

1.02355 

PNLA 

0.39E-01 

0. AAE+OO 

0. 16E-01 

0.77E-01 

0.38E-01 

0.37E+03 

0.52E+01 

0. A1E+01 

PNL1 

0.98990 

0.983A1 

0.98A22 

0.97030 

0.92177 

0.26059 

0.83033 

0.6A527 

PNL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL3 

0.99A90 

0.9915A 

0.99200 

0.98A78 

0.958A7 

0.29229 

0.90660 

0.77358 

PNLA 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNLS 

0.98E-06 

0.37E-0A 

0.86E-06 

0.1AE-0A 

O.A6ErOA 

0.62E+02 

0.22E-01 

0.10E+00 

QS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

1.02345 

1.8A3A6 

1.1A587 

2.10096 

2.51962 

76.72181 

5.A7180 

9.86203 

PDIS 

3.183A1 

3.15170 

3.2A885 

A. 12562 

A. 50007 

A7.79A32 

6.09939 

8.49522 

a)  Yaw  Axis  Measures  vs. 

Scenarios 

1 

2 

3 

4 

5 

6 

7 

8 

PDLA 

A. 39132 

3.93376 

9.01536 

1.10385 

2.22176 

10.95854 

1.95161 

4.76030 

PDL1 

1.13511 

1.27072 

1.01A26 

1.08A13 

1.00985 

1.12655 

0.67929 

1.57623 

PDL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL3 

1.05A80 

1.09117 

1.00682 

1.02876 

0.99876 

1.01628 

0.82061 

1.13016 

PDLA 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PDLS 

0.00252 

0.00782 

0.00013 

0.00066 

0.00062 

0.02317 

0.00167 

0.03384 

PSENA 

2.22AA6 

1.8A669 

3.201A0 

0.A8922 

0.68658 

5.4A358 

0.62359 

3.18802 

PSEN1 

1.18063 

1.39273 

1.01839 

1.13910 

1.07205 

1.23213 

0.75093 

1.60677 

PSEN2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

PSEN3 

1.07596 

1.13698 

1.00875 

1.05238 

1.02336 

1.07342 

0.86290 

1.16094 

PSENA 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PSENS 

0.001AA 

0.00A93 

0.00006 

0.00038 

0.00039 

0.01045 

0.00035 

0.020AA 

PNLA 

A. 52700 

A. 17595 

9. 0253 A 

1.155A1 

2.19557 

11.50537 

1.30719 

5.59484 

PNL1 

0.08683 

0.17A72 

0.02389 

0.07023 

0.0AA45 

0.11701 

0.08A69 

0.36638 

PNL2 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

PNL3 

0.203A1 

0.17851 

0.1AA07 

0.13A27 

0.1AAA8 

0.17943 

0 .42278 

0.24648 

PNLA 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNLS 

0.06681 

0.068A8 

0.1A919 

0.02038 

0.03697 

0.19322 

0.0099A 

0.09554 

GS1 

1.06651 

1.19282 

1.00386 

1.03A71 

1.00A66 

1.06264 

0.87321 

1.36615 

GS2 

1.02978 

1.08399 

1.00185 

1.01363 

1.00077 

1.02096 

0.93339 

1.13817 

GSS 

0.71E-03 

0.23E-02 

0.3AE-0A 

0.19E-03 

0.15E-03 

0.5AE-02 

0.19E-03 

0.01130 

PACC 

0.23E+02 

0.38E+02 

0.88E+01 

0.38E+02 

0.23E+02 

0. 15E+03 

0.37E+02 

96.41149 

PVEL 

0.29E+02 

0.A5E+02 

0.10E+02 

0.37E+02 

0.26E+02 

0.10E+03 

0.31E+02 

0.10E+03 

PDIS 

0.54E+02 

0.54E+02 

0. 25E+02 

0.62E+02 

0.A5E+02 

0.13E+03 

0.36E+02 

99.66941 

b)  Heave  Axis  Measures  vs.  Scenarios 


TABLE  C-4 


J 


i 


I 


PITCH  AND  SURGE  MEASURES  VERSUS  BREAK  FREQUENCY 

Kpp  = 0.0 


.0675 

.125 

.25 

.5 

1.0 

2.0 

PDLA 

25.60710 

25.60710 

25.60710 

25.60710 

25.60710 

25.60710 

PDL1 

0.58238 

0.80460 

1.16504 

1.20203 

1.09463 

1.03307 

PDL2 

0.00000 

O.OOOOC 

0.00000 

0.00000 

0.00000 

0.00000 

PDL3 

0.49643 

0.64237 

0.93280 

1.05512 

1.03705 

1.01463 

PDL4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDLS 

3.77401 

3.32806 

1.91691 

0.58764 

0.13143 

0.02441 

PSENA 

9.18581 

9.18581 

9.18581 

9.18581 

9.18581 

9.18581 

PSEN1 

0.15987 

0.49218 

1.16178 

1.35200 

1.19229 

1.07247 

PSEN2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN3 

0.09075 

0.29474 

0.78554 

1.07126 

1.07042 

1.03138 

PSEN4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSENS 

2.24677 

2.07303 

1.35653 

0.48104 

0.11815 

0.02230 

PNLA 

3.19468 

3.19468 

3.19468 

3.19468 

3.19468 

3.19468 

PNL1 

0.29299 

0.46650 

0.78778 

1.00181 

1.00512 

1.00107 

PNL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL3 

0.47081 

0.55040 

0.81790 

0.99801 

1.00234 

1.00053 

PNL4 

0.000 CO 

0 . 00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNLS 

0.28E+00 

0.29E+00 

0.12E+00 

0.46E-02 

0.34E-03 

0.10E-04 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000- 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

11.17123 

10.86064 

9.98160 

6.59412 

4.31896 

2.18626 

PDIS 

40.78298 

38.81647 

23.41467 

11.02253 

4.25406 

1.94151 

a)  Pitch  Axis  Measures  vs.  Break  Frequencies  (rod/sec) 


.0675 

.125 

.25 

.5 

1.0 

2.0 

PDLA 

0.01721 

0.01721 

0.01721 

0.01721 

0.01721 

0.01721 

PDL1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL2 

1.84815 

1.11629 

0.42077 

0.06658 

0.00652 

0.00052 

PDL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL4 

0.02261 

0.00913 

-0.08146 

-0.06634 

-0.02554 

-0.00663 

PDLS 

0.32E-01 

0.19E-01 

0.72E-02 

0.11E-02 

0.11E-03 

0.90E-05 

PSENA 

0.00525 

0.00525 

0.00525 

0.00525 

0.00525 

0.00525 

PSEN1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN2 

2.71456 

2.00411 

0.86713 

0.14935 

0.01491 

0.00114 

PSEN3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN4 

0.27864 

0.13130 

-0.09114 

-0.12497 

-0.05415 

-0.01431 

PSENS 

0.14E-01 

0.11E-01 

0.46E-02 

0.78E-03 

0.78E-04 

0.60E-05 

PNLA 

0.00050 

0.00050 

0.00050 

0.00050 

0.00Q50 

0 . 00050 

PNL1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL2 

0.50E+01 

0.26E+01 

0 • 43E+00 

0.S2E-O2 

0.15E-04 

0.11E-07 

PNL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL4 

0.64E+00 

0.36E+00 

0.27E-02- 

-0.23E-01- 

■0. 18E-02- 

-0.45E-04 

PNLS 

0.25E-02 

0.13E-02 

0.21E-03 

0.41E-05 

0.74E-08 

0.56E-11 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDIS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

b)  Surge  Axis  Measures  vs  Break  Frequencies  (rod /sec) 


179 


i 


1 


TABLE  C-5 


PITCH  AND  SURGE  MEASURES  VERSUS  SCENARIO 

K2P  = 0.5 


.0675 

.125 

.25 

.5 

1.0 

20 

PDLA 

25.60710 

25.60710 

25.60710 

25.60710 

25.60710 

25.60710 

°DL1 

0.58238 

0.80460 

1.16S04 

1.20203 

1.09463 

1.03307 

PDL2 

0.02788 

0.03552 

0.05757 

0.16119 

0.28006 

0.33848 

PDL3 

0.49643 

0.64237 

0.93280 

1.05512 

1.03705 

1.01463 

PDL4 

0.00623 

0.00024 

-0.00954 

0.02321 

0.07816 

0.10190 

POLS 

3.87278 

3.55240 

2.40768 

1.32235 

0.92362 

0.88666 

FSENA 

9.18581 

9.18581 

9.18581 

9.18581 

9.18581 

9.18581 

PSEN1 

0.15987 

0.49218 

1.16178 

1.35200 

1.19229 

1.07247 

PSEN2 

0.00036 

0.00220 

0.03911 

0.25779 

0.53792 

0.68060 

PSEN3 

0.09075 

0.29474 

0.78554 

1.07126 

1.07042 

1.03138 

PSEN4 

-0.00823 

-0.02505 

-0.05377 

0.00393 

0.12426 

0.17725 

PSENS 

2.28539 

2.19314 

1.69329 

1.05496 

0.78275 

0.77118 

PNLA 

3.19468 

3.19468 

3.19468 

3.19468 

3.19468 

3.19468 

PNL1 

0.29299 

0.46650 

0.78778 

1.00181 

1.00512 

1.00107 

PNL2 

-0.00649 

-0.01420 

-0.00737 

0.06551 

0.18484 

0.25358 

PNL3 

0.47081 

0.55040 

0.81790 

0.99801 

1.00234 

1.00053 

PNL4 

-0.00958 

-0.02711 

-0.05951 

-0.01521 

0.05584 

0.07938 

PNLS 

0.29074 

0.32403 

0.21055 

0.08123 

0.05877 

0.07573 

0S1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

11.34271 

11.30964 

11.46019 

9.68026 

7.91454 

9.02811 

PDIS 

42.11644 

44.01939 

31.98684 

23.17263 

21.10955 

20.85443 

0)  Pitch  Axis  Meosures 

us.  Break  Frequencies  ( 

(rod/sec) 

.0675 

.125 

.25 

.5 

1.0 

2.0 

PDLA 

0.01721 

0.01721 

0.01721 

0.01721 

0.01721 

0.01721 

PDL1 

0.54471 

0.38330 

0.30462 

0.22573 

0.10453 

0.03229 

PDL2 

8.43034 

5.12656 

1.90130 

0.31086 

0.02417 

0.00101 

PDL3 

0.57889 

0.38000 

0.27126 

0.17188 

0.07261 

0.02288 

PDL4 

0.04523 

0.01827 

-0.16291 

-0.13268 

-0.05109 

-0.01326 

PDLS 

6. 03755 

0.02458 

0.01286 

0.00627 

0.00467 

0.00436 

PSENA 

0.00525 

0.00525 

0.00525 

0.00525 

0.00525 

0.00525 

PSEN1 

0.74129 

0.67474 

0.66411 

0.53426 

0.24390 

0.07021 

PSEN2 

12.73880 

9.30336 

4.02552 

0.72150 

0.05737 

0.00205 

PSEN3 

0.72690 

0.61555 

0.53706 

0.36979 

0.14863 

0.04002 

PSEN4 

0.55729 

0.26261 

-0.18229 

-0.24995 

-0.10829 

-.0.02863 

PSENS 

0.01565 

0.01211 

0.00654 

0.00265 

0.00160 

0.00138 

PNLA 

0.00050 

0.00050 

0.00050 

0.00050 

0.00050 

0.00050 

PNL1 

0.99281 

0.98927 

0.98764 

0.94572 

0.76359 

0.53764 

PNL2 

27.16108 

16.63425 

4.88605 

0.36250 

-0.13047 

-0.06035 

PNL3 

0.99594 

0.99357 

0.99176 

0.96264 

0.84420 

0.71094 

PNL4 

1.21951 

0.50589 

-0.15706 

-0.23020 

-0.13887 

-0.04731 

PNLS 

0.31E-02 

0.20E-02 

0.65E-03 

0.11E-03 

0.28E-04 

0.19E-04 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

PDIS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

b)  Surge  Axis  Measures  us.  Break  Frequencies  (rod /sec) 
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TABLE  C-6 


PITCH  AND  SURGE  MEASURES  VERSUS  SCENARIO 
K2P  = 1.0 


.0675 

.125 

.25 

.5 

1.0 

2.0 

PDLA 

25.60710 

25.60710 

25.60710 

25.60710 

25.60710 

25.60710 

PDL1 

0.58238 

0.80460 

1.16504 

1.20203 

1.09463 

1.03307 

POL  2 

0.05846 

0.07709 

0.13621 

0.40432 

0.74707 

0.93538 

PDL3 

0.49643 

0.64237 

0.93280 

1.05512 

1.03705 

1.01463 

PDL4 

0.01246 

0.00048 

-0.01909 

0.04643 

0.15632 

0.20379 

PDLS 

3.98881 

3.31546 

3.03328 

2.58161 

2.91259 

3.40321 

PSENA 

9.18581 

9.18581 

9.18581 

9.18581 

9.18581 

9.18581 

PSEN1 

0.15987 

0.49218 

1.16178 

1.35200 

1.19229 

1.07247 

PSEN2 

0.00145 

0.00882 

0.11164 

0.69043 

1.51697 

1.98977 

PSEN3 

0.09075 

0.29474 

0.78554 

1.07126 

1.07042 

1.03138 

PSEN4 

-0.01646 

-0.05010 

-0.10754 

0.00787 

0.24852 

0.35450 

PSENS 

2.32570 

2.32341 

2.10681 

2.03044 

2.46040 

2.96353 

PNLA 

3.19468 

3.19468 

3.19468 

3.19468 

3.19468 

3.19468 

PNL1 

0.29299 

0.46650 

0.78778 

1.00181 

1.00512 

1.0C107 

PNL2 

-0.01259 

-0.02613 

0.01055 

0.37498 

1.21500 

1.89129 

PNL3 

0.47081 

0.55040 

0.81790 

0.99801 

1.00234 

1.00053 

PNL4 

-0.01937 

-0.05560 

-0.13136 

-0.05927 

0.16582 

0.31532 

PNLS 

0.30150 

0.36002 

0.33963 

0.39879 

0.70585 

1.00685 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

11.51418 

11.75864 

12.93877 

14.15045 

14.27637 

17.88819 

PDXS 

45.58288 

49.22908 

40 . 69039 

36.51659 

39.03415 

40.81094 

a)  Pitch  Axis  Measures  *s.  Break  Frequencies  (rad /sec) 


.0675 

.125 

.25 

.5 

1.0 

20 

PDLA 

0.01721 

0.01721 

0.01721 

0.01721 

0.01721 

0.01721 

PDL1 

0.54471 

0.38330 

0.30462 

0.22573 

0.10453 

0.03229 

PDL2 

0.24E+01 

0.14E+01 

0.53E+00 

0.89E-01 

0.56E-02- 

-0.19E-04 

PDL3 

0.57889 

0.38000 

0.27126 

0.17188 

0.07261 

0.02288 

PDL4 

0.02261 

0.00913 

-0.08146 

-0.06634 

-0.02554 

-0.00663 

PDLS 

0.04662 

0.03532 

0.02504 

0.01899 

0.01748 

0.01721 

PSENA 

0.00525 

0.00525 

0.00525 

0.00525 

0.00525 

0.00525 

PSEN1 

0.74129 

0.67474 

0.66411 

0.53426 

0.24390 

0.07021 

PSEN2 

3.65484 

2.64757 

1.14563 

0.21140 

0.01377 

-0.00011 

PSGN3 

0.72690 

0.61555 

0.53706 

0.36979 

0.14863 

0.04002 

PSEN4 

0.27864 

0.13130 

-0.09114 

-0.12497 

-0.05415 

-0.01431 

PSENS 

0.01779 

0.01486 

0.01008 

0.00660 

0.00562 

0.00535 

PNLA 

0.00050 

0.00050 

0.00050 

0.00050 

0.00050 

0.00050 

PNL1 

0.98246 

0.98219 

0.99727 

0.99779 

0.71798 

0.25056 

PNL2 

9.14106 

6.77975 

3.08539 

0.87414 

0.09893 

-0.00116 

PNI.3 

0.98405 

0.97474 

0.96753 

0.85732 

0.47550 

0.13738 

PNL4 

0.59223 

0.17472 

-0.13692 

-0.24333 

-0.15984 

-0.05162 

PNLS 

0.00400 

0.00325 

0.00172 

0.00083 

0.00060 

0.00054 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PD  IS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

b)  Surge  Axis  Measures  vs  Break  Frequencies  (rod /sec) 
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TABLE  C-7 

ROLL  AND  SWAY  MEASURES  VERSUS  BREAK  FREQUENCY 
k2R  = C2Y  = *0025  and  Kp  = 0.5 


.0675 


.125 


•25 


1.0 


2.0 


PDLA 

0.16E+03 

0.16E+03 

0.16E+03 

0. 16E+03 

0. 16E+03 

0.16E+03 

PDL1 

0.06404 

0.20797 

0.55343 

0.94366 

1.21167 

1.26063 

PDL2 

0.00000 

0 • 00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL3 

0.03543 

0.11937 

0.35261 

0.66463 

0.93627 

1.08140 

PDLA 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

POLS 

40.89412 

39.90810 

34.92462 

25.29759 

13.96307 

4.02845 

P SENA 

0.14E+03 

0.14E+03 

0. 14E+03 

0.14E+03 

0. 14E+03 

0.14E+03 

PSEN1 

0.04375 

0.14694 

0.44244 

0.88056 

1.22006 

1.29518 

PSEN2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN3 

0.02246 

0.08052 

0.27222 

0.59837 

0.91663 

1.09239 

PSEN4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSENS 

34.08286 

33.64171 

30.64212 

23.33425 

13.19895 

3.76732 

PNLA 

0. 1 1E+03 

0.  UE+03 

0. 11E+03 

0. 11E+03 

0.11E+03 

0. 11E+03 

PNL1 

0.04425 

0.14389 

0.43573 

0.83028 

1.07373 

1.13259 

PNL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL3 

0.13203 

0.19718 

0.37823 

0.66028 

0.90999 

1.04815 

PNLA 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNLS 

21.89133 

21.03117 

19.03984 

14.30231 

7.12027 

1.01812 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

55.80147 

56.66270 

48.01708 

42.05910 

32.62257 

15.77460 

PDIS 

85.77427 

87.7982 6 

83.26992 

52.10756 

22.10156 

10.66587 

a)  Roll  Axis  Measures  vs.  Break  Frequencies  (rad/sec) 


.0675 

.125 

.25 

• 5 

1.0 

2.0 

PDLA 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

PDL1 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL2 

12.53404 

8.28724 

2.93902 

0.46169 

0.04881 

0.00412 

PDL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL4 

-0.96844 

-0.73905 

-0.28230 

-0.06037 

-0.00681 

0.00612 

PDLS 

0.11E-01 

0.71E-02 

0.25E-02 

0.39E-03 

0.42E-04 

0.35E-05 

PSENA 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

PSEN1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN2 

11.64152 

8.78923 

3.49040 

0.60979 

0.06606 

0.00532 

PSEN3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN4 

-1.21305 

-0.98930 

-0.44138 

-0.10207 

0.00317 

0.01496 

PSENS 

0.51E-02 

0.38E-02 

0.15E-02 

0.27E-03 

0.29E-04 

0.23E-05 

PNLA 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

PNL1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL2 

0.38E+02 

0.33E+02 

0.25E+01 

0.53E-01 

0.11E-03 

0.94E-07 

PNL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL4 

-0.21E+01-0.21E+01-0. 14E+00- 

-0.30E-02 

0.17E-03 

0.31E-04 

PNLS 

0.11E-03 

0.91E-04 

0.69E-05 

0.15E-06 

0.32E-09 

0.26E-12 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

PACC 

53.82566 

47.30611 

48.32876 

27.41887 

11.68266 

6.28094 

PVEL 

0.17E+03 

0. 18E+03 

0.97E+02 

0.55E+02 

0.22E+02 

0.72E+01 

PDIS 

0.18E+04 

0.92E+03 

0.39E+03 

0.15E+03 

0.43E+02 

0. 11E+02 

b)  Sway  Axis  Measures  rs.  Break  Frequencies  (rad/sec) 
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TABIE  C-8 


ROLL  AND  SWAY  MEASURES  VERSUS  SCENARIO 


K2R  = 

0.5*  C2Y  = *0025,  and  Kp  = 0. 5 

.0675 

.125 

.25 

• 5 

1.0 

2.0 

PDLA 

0. 16E+03 

0.16E+03 

0.16E+03 

0.16E+03 

0. 16E+03 

0. 16E+03 

PDL1 

0.06404 

0.20797 

0.55343 

0.94366 

1.21167 

1.26063 

PDL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL3 

0.03S43 

0.11937 

0.35261 

0.66463 

0.93627 

1.08140 

PDL4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDLS 

40.89412 

39.90810 

34.92462 

25.29759 

13.96307 

4.02845 

PSENA 

0. 14E+03 

0.14E+03 

0.14E+03 

0. 14E+03 

0.14E+03 

0. 14E+03 

PSEN1 

0.04375 

0.14694 

0.44244 

0.88056 

1.22006 

1.29518 

PSEN2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN3 

0.02246 

0.08052 

0.27222 

0.59837 

0.91663 

1.09239 

PSEN4 

0.00000 

0.00000 

0.00000 

0.00000 

0 • 00000 

0.00000 

PSENS 

34.08286 

33.64171 

30.64212 

23.33425 

13.19895 

3.76732 

PNLA 

0.11E+03 

0.11E+03 

0.11E+03 

0.UE+03 

0.11E+03 

0.11E+03 

PNL1 

0.04425 

0.14389 

0.43573 

0.83028 

1.07373 

1.13259 

PNL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL3 

0.13203 

0.19718 

0.37823 

0.66028 

0.90999 

1.04815 

PNL4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNLS 

21.89133 

21.03117 

19.05984 

14.30231 

7.12027 

1.01812 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

55.80147 

56.66270 

48.01708 

42.05910 

32.62257 

15.77460 

PDIS 

85.77427 

87.79826 

83.26992 

52.10756 

22.10156 

10.66587 

a)  Roll  Axis  Measures  ws.  Break  frequencies  (rad/sec) 

.0675 

.125 

• 25 

• 5 

1.0 

2.0 

PDLA 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

PDL1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL2 

78.33399 

62.59290 

24.50791 

4.07210 

0.43979 

0.03723 

PDL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL4 

-2.24497 

-1.50476 

-0.37052 

-0.03195 

0.01814 

0.02822 

PDLS 

0.67E-01 

0.S3E-01 

0.21E-01 

0.35E-02 

0.38E-03 

0.32E-04 

PSENA 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

PSEN1 

0.00000 

0.00000 

0 • 00000 

0.00000 

0.00000 

0.00000 

PSEN2 

83.83268 

69.90791 

29.67314 

5.47468 

0.60242 

0.04851 

PSEN3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN4 

-3.01640 

-2.19215 

-0.73313 

-0.09376 

0.06814 

0.05917 

PSENS 

0.37E-01 

0.31E-01 

0.13E-01 

0.24E-02 

0.26E-03 

0.21E-04 

PNLA 

0.28E-05 

0.28E-05 

0.2BE-05 

0.28E-05 

0.28E-05 

0.28E-05 

PNL1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL2 

0.47E+04 

0. 36E+04 

0.90E+03 

0.38E+02 

0.69E-01 

0.55E-04 

PNL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL4 

-19.74601- 

-11.75227 

-0.06778 

0.04731 

0.01714 

0.00106 

PNLS 

0.13E-01 

0.10E-01 

0.25E-02 

0.11E-03 

0.19E-06 

0.15E-09 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

PACC 

52.76090 

56.02879 

43.55520 

21.44993 

12.65569 

6.79812 

PVEL 

0. 16E+03 

0. 13E+03 

0.11E+03 

0.60E+02 

0.21E+02 

0.68E+01 

PDIS 

0.73E+03 

0. 68E+03 

0.36E+03 

0.13E+03 

0.40E+02 

0.11E+02 

b)  Sway  Axis  Measures  rs.  Break  Frequencies  (rad/sec) 
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TABLE  C-9 


ROLL  AND  SWAY  MEASURES  VERSUS  SCENARIO 


K2R  = 

1.0,  C2y  - .0025,  and  Kj 

, = 0.5 

.0675 

.125 

.25 

•5 

t.O 

2.0 

PDLA 

O. 16E+03 

0.16E+03 

0.16E+03 

0. 16E+03 

0. 16E+03 

0.16E+03 

PDL1 

0.04404 

0.20797 

0.55343 

0.94366 

1.21167 

1.26063 

P0L2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL3 

0.03543 

0.11937 

0.35261 

0.66463 

0.93627 

1.08140 

PDL4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

POLS 

40.89412 

39.90810 

34.92462 

25.29759 

13.96307 

4.02845 

PSENA 

0.14E+03 

0.14E+03 

0.14E+03 

0 . 14E+03 

0. 14E+03 

0.14E+03 

PSENl 

0.04375 

0.14694 

0.44244 

0.88056 

1.22006 

1.29518 

PSEN2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 • 00000 

PSEN3 

0.02244 

0.08052 

0.27222 

0.59837 

0.91663 

1.09239 

PSEN4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSENS 

34.08284 

33.64171 

30.64212 

23.33425 

13.19895 

3.76732 

PNLA 

0. UE+03 

0.  UE+03 

0.  UE+03 

0.  UE+03 

0.  UE+03 

0.  UE+03 

PNL1 

0.04425 

0.14389 

0.43573 

0.83028 

1.07373 

1.13259 

PNL2 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PNL3 

0.13203 

0.19718 

0.37823 

0.66028 

0.90999 

1.04815 

PNL4 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PNLS 

21.89133 

21.03117 

19.05984 

14.30231 

7.12027 

1.01812 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0S2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

OSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

55.80147 

56.66270 

48.01708 

42.05910 

32.62257 

15.77460 

PDIS 

85.77427 

87.79826 

83.26992 

52.10756 

22.10156 

10.66587 

a)  Roll  Axis  Measures  vs.  Break  Frequencies  (rad/sec) 

.0675 

.125 

.25 

• 5 

1.0 

2.0 

POLA 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

PDL1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL2 

0.13E+03 

0.11E+03 

0 .62E+02 

0.19E+02 

0.30E+01 

0.30E+00 

PDL3 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PDL4 

-2.29020 

-1.18560 

-0.09902 

0.14211 

0.09097 

0.10069 

PDLS 

0.10960 

0.09541 

0.05300 

0.01653 

0.00256 

0.00025 

PSENA 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

PSEN1 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PSEN2 

0. 16E+03 

0. 14E+03 

0.84E+02 

0.28E+02 

0.42E+01 

0.39E+00 

PSEN3 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PSEN4 

-3.45136 

-2.04522 

-0.36981 

0.36459 

0.34690 

0.20934 

PSENS 

0.06809 

0.06083 

0.03683 

0.01232 

0.00185 

0.00017 

PNLA 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-0S 

0.28E-05 

PNL1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 • 00000 

PNL2 

0. 14E+05 

0. 12E+05 

0.59E+04 

0. 10E+04 

0.88E+01 

0.10E-01 

PNL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL4 

-12.83748 

-1.73929 

5.01623 

6.57397 

0.26171 

0.01518 

PNLS 

0.38E-01 

0.32E-01 

0.16E-01 

0.29E-02 

0.25E-04 

0.28E-07 

OS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

OS2 

0.00000 

0.00000 

0 • 00000 

0.00000 

0 • 00000 

0.00000 

GSS 

0.00000 

0.00000 

0 • 00000 

0.00000 

0.00000 

0.00000 

PACC 

26.57902 

23.12729 

21.70629 

19.97275 

12.70632 

6.24403 

PVEL 

56.81534 

60.44103 

49.64578 

25.94772 

12.70869 

4.65805 

PDIS 

0.16E+03 

0. 14E+03 

0. 12E+03 

0.70E+02 

0.24E+02 

0.69E+01 

b)  Seay  Axis  Measures  vs.  Break  Frequencies  (rad/sec) 
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TABLE  C-10 


ROLL  AND  SWAY  MEASURES  VERSUS 

SCENARIO 

K2r  = 

0.5>  C2Y  = .0225,  and  Kp  = 0.5 

.0675 

.125 

• 25 

• 5 

1.0 

2.0 

PDLA 

0.1 6E+03 

0.16E+03 

0.16E+03 

0. 16E+03 

0.16E+03 

0.16E+03 

PDL1 

0.06 AO A 

0.20797 

0.55343 

0.94366 

1.21167 

1.26063 

PDL2 

0.00038 

0.00207 

0.00860 

0.01730 

0.01024 

-0.04193 

PDL3 

0.035 A3 

0.11937 

0.35261 

0.66463 

0.93627 

1.08140 

PDLA 

0.0006A 

0.00250 

0.00942 

0.02325 

0.03523 

0.00700 

PDLS 

A0.857A3 

39.78760 

34.50283 

24.09550 

11.48383 

1.72578 

PSENA 

0.1AE+03 

0.14E+0 3 

0.14E+03 

0.14E+03 

0.14E+03 

0.14E+03 

PSEN1 

0.0A375 

0.1A694 

0.44244 

0.88056 

1.22006 

1.29518 

PSEN2 

0.00028 

0.00158 

0.00749 

0.01753 

0.01141 

-0.04732 

PSEN3 

0.022A6 

0.08052 

0.27222 

0.59837 

0.91663 

1.09239 

PSENA 

0.00058 

0.00229 

0.00919 

0.02502 

0.04007 

0.00791 

PSENS 

3 A. 05270 

33.53963 

30.27024 

22.22493 

10.85342 

1.61305 

PNLA 

0. 1 1E+03 

0.11E+03 

0.11E+03 

0.11E+03 

0.11E+03 

0.11E+03 

PNL1 

0.0AA25 

0.14389 

0.43573 

0.83028 

1.07373 

1.13259 

PNL2 

0.00036 

0.00193 

0.00847 

0.01608 

0.00942 

-0.04423 

PNL3 

0.13203 

0.19718 

0.37823 

0.66028 

0.90999 

1.04815 

PNLA 

0.00066 

0.00247 

0.00928 

0.02330 

0.03542 

-0.00741 

PNLS 

21.86A15 

20.94668 

18.77669 

13.44626 

5.39690 

0.19314 

GSt 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

55.78326 

56.56514 

47.50300 

41.55858 

29.52164 

10.22700 

PDXS 

85 . 7652A 

87.63050 

83.00689 

49.94592 

18.74533 

9.31515 

a)  Roll  Axis  Measures  vs.  Break  Frequencies  (rad/sec) 

.0675 

.125 

.25 

• 5 

1.0 

2.0 

PDLA 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

PULI 

0.13388 

0.1542A 

0.15556 

0.09845 

0.04766 

0.01888 

PDL2 

0.31E+03 

0.25E+03 

0.10E+03 

0.17E+02 

0.20E+01 

0. 19E+00 

PDL3 

0.07317 

0.08893 

0.12110 

0.13250 

0.12665 

0.10854 

PDL4 

-A.4899A 

-3.00953 

-0.74104 

-0.06390 

0.03628 

0.05644 

PDLS 

0.06930 

0.05573 

0.02222 

0.00394 

0.00058 

0.00019 

PSENA 

0.0004A 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

PSEN1 

0.16723 

0.19347 

0.20283 

0.13402 

0.06113 

0.02154 

PSEN2 

0.34E+03 

0.28E+03 

0.12E+03 

0.23E+02 

0.28E+01 

0.25E+00 

PSEN3 

0.09116 

0.11134 

0.15772 

0.18366 

0.16938 

0.12879 

PSENA 

-6.03280 

-4.38430 

-1.46627 

-0.18753 

0.13628 

0.11834 

PSENS 

0.0382A 

0.03208 

0.01388 

0.00270 

0.00035 

0.00009 

PNLA 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

PNL1 

0.71289 

0.73721 

0.77065 

0.76175 

0.71571 

0.67643 

PNL2 

0.20E+05 

0.16EA05 

0.42E+04 

0.24E+03 

0.13E+01 

0. 12E+00 

PNL3 

0.8A135 

0.85539 

0.87444 

0.87092 

0.84537 

0.82225 

PNLA 

-45.74344- 

-29.58179 

-2.27563 

-0.23843 

0.06402 

0.06728 

PNLS 

0.1AE-01 

0.1  IE-01 

0.29E-02 

0.17E-03 

0.81E-06 

0.13E-07 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

50.02A23 

53.44072 

42.96724 

21.40831 

12.77244 

8.12718 

PVEL 

0.16E+03 

0. 13E+03 

0.UE+03 

0.61E+02 

0.23E+02 

0.80E+01 

PDIS 

0.73E+03 

0.69E+03 

0.36E+03 

0.13E+03 

0.44E+02 

0.13E+02 

b)  Sway  Axis  Measures  ts.  Break  Frequencies  (rad/sec) 
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TABLE  C-11 


ROLL  AND  SWAY  MEASURES  VERSUS  SCENARIO 
K2R  = O.O,  C2Y  = 0.:,  and  Kp  = 0.5 


.0675 

.125 

•25 

.5 

1.0 

2.0 

y 

PDLA 

0.16E+03 

0. 16E+03 

0. 16E+03 

0.16E+03 

0. 16C+03 

0.16E+03 

PDL1 

0.06404 

0.20797 

0.55343 

0.94366 

1.21167 

1.26063 

i 

PDL2 

0.00077 

0.00416 

0.01742 

0.03621 

0.03016 

-0.04607 

PDL3 

0.03543 

0.11937 

0.35261 

0.66463 

0.93627 

1.08140 

i 

PDL4 

0.00127 

0.00499 

0.01885 

0.04650 

0.07045 

0.01399 

r 

PDLS 

40.82082 

39.66796 

34.09004 

22.95970 

9.40313 

0.97914 

PSENA 

0. 14E+03 

0. 14E+03 

0. 14E+03 

0. 14E+03 

0. 14E+03 

0.14E+03 

l 

PSEN1 

0.04375 

0.14694 

0.44244 

0.88056 

1.22006 

1.29518 

PSEN2 

0.00057 

0.00317 

0.01518 

0.03677 

0.03372 

-0.05190 

PSEN3 

0.02246 

0.08052 

0.27222 

0.59837 

0.91663 

1.09239 

L 

PSEN4 

0.00117 

0.00457 

0.01839 

0.05004 

0.08014 

0.01581 

1 

PSEN8 

34.02260 

33.43813 

29.90532 

21.17370 

8.87997 

0.91714 

PNLA 

0.  UE+03 

0.  UE+03 

0.  UE+03 

0.  UE+03 

0.  UE+03 

0.  UE+03 

PNL1 

0.04425 

0.14389 

0.43573 

0.83028 

1.07373 

1.13259 

j) 

PNL2 

0.00072 

0.00389 

0.01712 

0.03307 

0.02789 

-0.07058 

PNL3 

0.13203 

0.19718 

0.37823 

0.66028 

0.90999 

1.04815 

PNL4 

0.00133 

0.00495 

0.01855 

0.04618 

0.06958 

-0.01901 

PNL8 

21.83703 

20.86272 

18.49943 

12.63848 

3.99825 

0.10469 

i 

OS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

i 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

B 

PVEL 

55.76504 

56.46758 

46.98891 

41.05805 

26.72289 

8.58769 

PDIS 

85.75620 

87.46274 

82.74385 

47.78429 

18.30605 

9.33734 

a)  Roll  Axis  Measures  vs.  Break  Frequencies  (rad/sec) 


.0675 

• 125 

.25 

• 5 

1.0 

2.0 

PDLA 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

PDL1 

0.01630 

0.01819 

0.01745 

0.01078 

0.00522 

0.00208 

PDL2 

12.60749 

8.51663 

3.15231 

0.52832 

0.06196 

0.00641 

PDL3 

0.00956 

0.01481 

0.02934 

0.03886 

0.04014 

0.03564 

PDL4 

-0.96844 

-0.73905 

-0.28230 

-0.06037 

-0.00681 

0.00612 

PDLS 

0.01328 

0.00939 

0.00400 

0.00135 

0.00085 

0.00079 

PSENA 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

PSEN1 

0.01942 

0.02208 

0.02238 

0.01456 

0 . 00666 

0.00236 

PSEN2 

11.70503 

9.02788 

3.73847 

0.69734 

0.08480 

0.00857 

PSEN3 

0.01091 

0.01737 

0.03712 

0.05339 

0.05367 

0.04233 

PSEN4 

-1.21305 

-0.98930 

-0.44138 

-0.10207 

0.00317 

0.01496 

PSENS 

0.00662 

0.00525 

0.00244 

0.00079 

0.00043 

0.00039 

PNLA 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

PNL1 

0.01955 

0.02490 

0.03930 

0.04222 

0.02530 

0.01345 

PNL2 

0.12E+03 

0.  UE+03 

0.10E+02 

0.67E+00 

0.10E+00 

0.22E-01 

PNL3 

0.09081 

0.11454 

0.17182 

0.18659 

0.14820 

0.11185 

PNL4 

-7.91721 

-6.58732 

-1.39134 

0.11111 

0.13706 

0.05968 

PNLS 

0.37E-03 

0.35E-03 

0.38E-04 

0.31E-05 

0.16E-05 

0.19E-05 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 • 00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

51.23702 

45.25141 

48.24283 

28.31316 

14.79704 

10.28296 

PVEL 

0.17E+03 

0.18E+03 

0.10E+03 

0.S8E+02 

0.25E+02 

0.10E+02 

PDIS 

0.18E+04 

0.94E+03 

0.40E+03 

0.17E+03 

0.52E+02 

0.16E+02 

r 
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b)  Sway  Axis  Measures  vs.  Break  Frequencies  (rad/sec) 


TABLE  C-12 


ROLL  AND  SWAY  JEASURES 

VERSUS  SCENARIO 

K2r  = 

: 0-5,  C2Y  = °*5 y 

and  Kp 

= 0.5 

.0675 

• 125 

.25 

.5 

1.0 

2.0 

PDLA 

0.16E+03 

0.16E+03 

0. 16E+03 

0.16E+03 

0.16E-F03 

0. 16E+03 

PDL1 

0.06404 

0.20797 

0.55343 

0.94366 

1.21167 

1.26063 

PDL2 

0.00077 

0.00416 

0.01742 

0.03621 

0.03016 

-0.04607 

PDL3 

0.03543 

0.11937 

0.35261 

0.66463 

0.93627 

1.08140 

PDL4 

0.00127 

0.00499 

0.01885 

0.04650 

0.07045 

0.01399 

PDLS 

40.82082 

39.66796 

34.09004 

22.95970 

9.40313 

0.97914 

PSENA 

0. 14E+03 

0. 14E+03 

0. 14E+03 

0.  14E+03 

0. 14E+03 

0. 14E+03 

PSEN1 

0.04375 

0.14694 

0.44244 

0.88056 

1.22006 

1.29518 

PSEN2 

0.00057 

0.00317 

0.01518 

0.03677 

0.03372 

-0.05190 

PSEN3 

0.02246 

0.08052 

0.27222 

0.59837 

0.91663 

1.09239 

PSEN4 

0.00117 

0.00457 

0.01839 

0.05004 

0.08014 

0.01581 

PSENS 

34.02260 

33.43813 

29.90532 

21.17370 

8.87997 

0.91714 

PNLA 

0.11E+03 

0. 11E+03 

0.1 1E+03 

0.  UE+03 

0.  UE+03 

0.  UE+03 

PNL1 

0.04425 

0.14389 

0.43573 

0.83028 

1.07373 

1.13259 

PNL2 

0.00072 

0.00389 

0.01712 

0.03307 

0.02789 

-0.07058 

PNL3 

0.13203 

0.19718 

0.37823 

0.66028 

0.90999 

1.04815 

PNL4 

0.00133 

0.00495 

0.01855 

0.04618 

0.06958 

-0.01901 

PNLS 

21.83703 

20.86272 

18.49943 

12.63848 

3.99825 

0.10469 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

55.76504 

56.46758 

46.98891 

41.05805 

26.72289 

8.58769 

PDIS 

85.75620 

87.46274 

82.74385 

47.78429 

18.30605 

9.33734, 

a)  Roll  Axis  Measures  ts.  Break  Frequencies  (rad/sec) 

.0675 

• 125 

.25 

•5 

1.0 

2.0 

PDLA 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

PDL1 

0.13388 

0.15424 

0.15556 

0.09845 

0.04766 

0.01888 

PDL2 

78.87490 

64.35049 

26.29375 

4.66591 

0.55924 

0.05810 

PDL3 

0.07317 

0.08893 

0.12110 

0.13250 

0.12665 

0.10854 

PDL4 

-2.24497 

-1.50476 

-0.37052 

-0.03195 

0.01814 

0.02822 

PDLS 

0.07208 

0.05839 

0.02388 

0.00475 

0.00113 

0.00069 

PSENA 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

PSEN1 

0.16723 

0.19347 

0.20283 

0.13402 

0.06113 

0.02154 

FSEN2 

84.27433 

71.7588 7 

31.80060 

6.27024 

0.77469 

0.07833 

PSEN3 

0.09116 

0.11134 

0.15772 

0.18366 

0.16938 

0.12879 

PSEN4 

-3.01640 

-2.19215 

-0.73313 

-0.09376 

0.06814 

0.05917 

PSENS 

0.03998 

0.03377 

0.01496 

0.00316 

0.00060 

0.00032 

PNLA 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

PNL1 

0.22018 

0.27235 

0.34915 

0.29266 

0.17065 

0.09692 

PNL2 

0.54E+04 

0.43E+04 

0. 12E+04 

0.92E+02 

0.13E+01 

0. 14E+00 

PNL3 

0.37808 

0.43235 

0.50626 

0.49243 

0.39319 

0.30342 

PNL4 

-25.58068- 

■17.39147 

-2.15347 

0.04884 

0.26198 

0.15039 

PNLS 

0.15E-01 

0.12E-01 

0.34E-02 

0.26E-03 

0.33E-05 

0.91E-06 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0 • 00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

GSS 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PACC 

48.85858 

53.79541 

43.62885 

23.17445 

14.18067 

10.37257 

PVEL 

0.16E+03 

0.13E+03 

0. 1 1E+03 

0.62E+02 

0.2SE+02 

0.98E+01 

PDIS 

0.73E+03 

0.69E+03 

0.37E+03 

0.14E+03 

0.48E+02 

0.16E+02 

b)  Sway  Axis  Measures  rs.  Break  Frequencies  (rad/sec) 
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T.iBLE  C-13 

ROLL  AND  SWAY  NEASURES  VERSUS  SCENARIO 
KgR  = 1.0,  C2Y  = 0.5,  and  Kp  = 0.5 


.0675 

.125 

•25 

•5 

1.0 

2.0 

PDLA 

0. 16E+03 

0.16E+03 

0.16E+03 

0. 16E+03 

0.16E+03 

0.16E+03 

PDL1 

0.06404 

0.20797 

0.55343 

0.94366 

1.21167 

1.26063 

PDL2 

0.00077 

0.00416 

0.01742 

0.03621 

0.03016 

-O; 04607 

PDL3 

0.03SA3 

0.11937 

0.35261 

0.66463 

0.93627 

1.08140 

PDLA 

0.00127 

0.00499 

0.01885 

0.04650 

0.07045 

0.01399 

POLS 

AO. 82082 

39.66796 

34.09004 

22.95970 

9.40313 

0.97914 

PSENA 

0. 14E+03 

0. 14E+03 

0. 14E+03 

0. 14E+03 

0. 14E+03 

0.14E+03 

PSEN1 

0.04375 

0.14694 

0.44244 

0.88056 

1.22006 

1.29518 

PSEN2 

0.00057 

0.00317 

0.01518 

0.03677 

0.03372 

-0.05190 

PSEN3 

0.02246 

0.08052 

0.27222 

0.59837 

0.91663 

1.09239 

PSENA 

0.00117 

0.00457 

0.01839 

0.05004 

0.08014 

0.01581 

PSENS 

3A. 02260 

33.43813 

29,90532 

21.17370 

8.87997 

0.91714 

PNLA 

0.1 1E+03 

0.11E+03 

0.  UE+03 

0.  UE+03 

0.  UE+03 

0.  UE+03 

PNL1 

0.0AA25 

0.14389 

0.43573 

0.83028 

1.07373 

1.13259 

PNL2 

0.00072 

0.00389 

0.01712 

0.03307 

0.02789 

-0.07058 

PNL3 

0.13203 

0.19718 

0.37823 

0.66028 

0.90999 

1.04815 

PNLA 

0.00133 

0.00495 

0.01855 

0.04618 

0.06958 

-0.01901 

PNLS 

21.83703 

20.86272 

18.49943 

12.63848 

3.99825 

0.10469 

GS1 

0.00000 

0.00000 

0 . 00000 

0.00000 

O.COOOO 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

55.76504 

56.46758 

46.98891 

41.05805 

26.72289 

8.58769 

PDIS 

85.75620 

87.46274 

82.74385 

47.78429 

18.30605 

9.33734 

a)  Boll  Axis  Measures  rs.  Break  Frequencies  (rad/sec) 


.0675 

•125 

.25 

•5 

1.0 

2.0 

PDLA 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

0.00085 

PDL1 

0.58151 

0.63654 

0.72090 

0.67231 

0.42521 

0.18765 

PDL2 

0.13E+03 

0. UE+03 

0.65E+02 

0.21E+02 

0.36E+01 

0.43E+00 

PDL3 

0.37263 

0.39351 

0.41913 

0.41995 

0.38547 

0.32607 

PDLA 

-2.29020 

-1.18560 

-0.09902 

0.14211 

0.09097 

0.10069 

PDLS 

0.11265 

0.09868 

0.05612 

0.01856 

0.00345 

0.00065 

PSENA 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

0.00044 

PSEN1 

0.77512 

0.85079 

0.98157 

0.92471 

0.53516 

0.20753 

PSEN2 

0.15E+03 

0.14E+03 

0.87E+02 

0.31E+02 

0.S1E+01 

0.59E+00 

PSEN3 

0.50260 

0.53212 

0.57781 

0.58743 

0.49624 

0.36670 

PSEN4 

-3.45136 

-2.04522 

-0.36981 

0.36459 

0.34690 

0.20934 

PSENS 

0.07012 

0.06305 

0.03898 

0.01355 

0.00218 

0.00028 

PNLA 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

0.28E-05 

PNL1 

1.08436 

1.14903 

1.04723 

0.78052 

0.54066 

0.40748 

PNL2 

0.14E+05 

0.12E+05 

0.64E+04 

0.12E+04 

0.19E+02 

0.26E+00 

PNL3 

0.84317 

0.85745 

0.86095 

0.84389 

0.71683 

0.62937 

PNL4 

-14.90417 

-3.91391 

3.52785 

4.37840 

0.14775 

0.16286 

PNLS 

0.39E-01 

0.34E-01 

0.18E-01 

0.34E-02 

0.54E-04 

0.22E-06 

GS1 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

24.49761 

21.74777 

21.37760 

19.47116 

15.87156 

9.51661 

PVEL 

54.30840 

58.02134 

48.68298 

27.54009 

14.45196 

6.83523 

PDIS 

0. 16Ef03 

0.14E+03 

0.12E+03 

0.75E+02 

0.28E+02 

0.10E+02 

b)  Sway  Axis  Measures  vs.  Break  Frequencies  (rad/sec) 
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TABUS  C-lU 


YAW  AMD  HEAVE  AXIS  NEASURES  VERSUS  BREAM  FREQUENCY 

.0675 

.125 

.25 

•5 

1.0 

2.0 

PDLA 

5.29930 

5.29930 

5.29930 

5.29930 

5.29930 

5.29930 

PDH 

0.65466 

0.93944 

1.18190 

1.14679 

1.07502 

1.03102 

PDL2 

0.00000 

0.00000 

0 . 00000 

0 . 00000 

0.00000 

0 • 00000 

PDL3 

0.54160 

0.74246 

0.98077 

1.03396 

1.02612 

1.01335 

PDL4 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

POLS 

0.75709 

0.60215 

0.29193 

0.10448 

0.03017 

0.00572 

PSENA 

1.63021 

1.63021 

1.63021 

1.63021 

1.63021 

1.63021 

PSEN1 

0.21707 

0.59527 

1 . 1S282 

1.26259 

1.17356 

1.08087 

PSEN2 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PSEN3 

0.13114 

0.38600 

0.82608 

1.02317 

1.05268 

1.03382 

PSENA 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSENS 

0.38913 

0.33553 

0.20405 

0.08813 

0.02779 

0.00539 

PNLA 

0.07667 

0.07667 

0.07667 

0.07667 

0.07667 

0.07667 

PNL1 

0.50633 

0.66803 

0.88223 

0.97459 

1.00178 

1.00108 

PNL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL3 

0.70278 

0.79854 

0.93250 

0.98675 

1.00086 

1.00054 

PNLA 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNLS 

0.19E-02 

0.14E-02 

0.33E-03 

0.21E-04 

0.13E-05 

0.26E-07 

GS1 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

GS2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0 . 00000 

GSS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

3.87879 

4.36360 

3*66914 

2.31082 

1.44214 

0.72981 

PDIS 

26.99819 

14.55148 

7.93515 

3.92446 

1.84893 

0.69747 

a)  Yaw  Axis  Measures  vs.  Break  Frequencies  (red./ sec) 

| 

.0675 

.125 

• 25 

•5 

1.0 

2.0 

PDLA 

1.10385 

1.10385 

1.10385 

1.10385 

1.10385 

1.10385 

PDL1 

0.83764 

1.21807 

1.55947 

1 . 18993 

1.06558 

1.01853 

PDL2 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

PDL3 

0.63075 

0.85105 

1.12942 

1.05578 

1.02024 

1.00642 

PDLA 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

' PDLS 

0.01431 

0.01281 

0.00747 

0.00195 

0.00062 

0.00014 

PSENA 

0.48922 

0.48922 

0.48922 

0.48922 

0.48922 

0.48922 

PSEN1 

0.65404 

1.19463 

1.83484 

1.32172 

1.10502 

1.02095 

PSEN2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSEN3 

0.42263 

0.72214 

1.18275 

1.10095 

1 . 03699 

1.00848 

PSEN4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PSENS 

0.89E-02 

0.83E-02 

0.52E-02 

0.13E-02 

0.34E-03 

0.44E-04 

PNLA 

1.15541 

1.15541 

1.15541 

1.15541 

1.15541 

1.15541 

PNL1 

0.02100 

0.04448 

0.08403 

0.07643 

0.07223 

0.05933 

PNL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNL3 

0.16689 

0.11625 

0.08720 

0.13090 

0.13625 

0.14512 

PNLA 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PNLS 

0.02217 

0.02412 

0.02404 

0.02097 

0.02022 

0.01978 

GS1 

1.01189 

1.11171 

1.15905 

1.05259 

1.03069 

1.00686 

GS2 

0.98101 

1.03165 

1.06070 

1.01700 

1.01165 

1.00278 

GSS 

0.13E-02 

0.13E-02 

0.98E-03 

0.48E-03 

0.19E-03 

0.34E-04 

PACC 

0.96E+02 

0.10E+03 

0.87E+02 

0.56E+02 

0.32E+02 

0.26E+02 

PVEL 

0.28E+03 

0.23E+03 

0.19E+03 

0.58E+02 

0.30E+02 

0. 10E+02 

PDIS 

0.20E+04 

0.14E+04 

0.61E+03 

0.16E+03 

0.40E+02 

0.94E+01 

b)  Heave  Axis  Measures  vs.  Break  Frequencies  (rad /sec) 
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TABLE  C-15.  . REVISED  MEASURES,  PITCH-SURGE  AXES  VERSUS 
BREAK  FREQUENCY  AND  Kgp 


BF 

.0675 

•125 

•25 

•5 

1.0 

2.0 

PDLl 

0.58238 

0.80460 

1.16504 

1.20203 

1.09463 

1.03307 

PDL2 

0.00000 

0.00000 

0 • 00000 

0.00000 

0.00000 

0.00000 

PDL3 

PDL4 

0.49643 

0.00000 

0.64237 

0.00000 

0.93280 

0.00000 

1.05512 

0.00000 

1.03705 

0.00000 

1.01463 

0.00000 

Jq  measures. 

PVEL 

11.17123 

10.86064 

9.98160 

6.59412 

4.31896 

2.18626 

k2P  =0.0 

PDIS 

40.78298 

38.81647 

23.41467 

11.02253 

4.25406 

1.94151 

PDLl 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

0.00000 

P0L2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

Jx  measures. 

POL  3 

0.00000 

0 . 00000 

0.00000 

0.00000 

0.00000 

0.00000 

H0L4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

Kgp  ~ 0.0 

PACO 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PVEL 

0.00000 

0.00000 

0.00000 

0 . 00000 

0.00000 

0.00000 

- 

PDIS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDLl 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

JXT  measures. 

PDL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

Kgp  — 0.0 

°DL1 

0.58238 

0.80460 

1.16504 

1.20203 

1.09463 

1.03307 

PDL2 

0.02788 

0.03552 

0.05757 

0.16119 

0.28006 

0.33848 

Jq  measures. 

PDL3 

0.49643 

0.64237 

0.93280 

1.05512 

1.03705 

1.01463 

PDL4 

0.00623 

0.00024 

-0.00954 

0.02321 

0.07816 

0.10190 

Kgp  = 0.5 

PVEL 

11.34271 

11.30964 

11.46019 

9.68026 

7.91454 

9.02811 

PDIS 

42.11644 

44,01939 

31.98684 

23.17263 

21.10955 

20.85443 

1 

PDLl 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

Jx  measures. 

PDL4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

I 

PACC 

PVEL 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

K2P  =0.5 

PDIS 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDLl 

0.54471 

0.38330 

0.30462 

0.22573 

0.10453 

0.03229 

Jxr  measures, 

PDL3 

0.57889 

0.38000 

0.27126 

0.17188 

0.07261 

0.02288 

Kgp  = 0.5 

PDLl 

PDL2 

0.58238 

0.05846 

0.80460 

0.07709 

1 . 16504 
0.13621 

1.20203 

0.40432 

1.09463 

0.74707 

1.03307 

0.93538 

Jq  measures, 

Kgp  = 1.0 

PDL3 

PDL4 

0.49643 

0.01246 

0.64237 

0.00048 

0.93280 

-0.01909 

1.05512 

0.04643 

1.03705 

0.15632 

1.01463 

0.20379 

1 

i 

PVEL 

11.51418 

11.75864 

12.93877 

14.15045 

14.27637 

17.88819 

PDIS 

45.58288 

49.22908 

40.69039 

36.51659 

39.03415 

40.81094 

PDLl 

PDL2 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

I 

PDL3 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

Jy  measures. 

PDL4 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PACC 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

Kgp  = 1.0 

1 

PVEL 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

PDIS 

0.00000 

0.00000 

0.00000 

0.00000 

0 • 00000 

0.00000 

. 

PDLl 

0.54471 

0.38330 

0.30462 

0.22573 

0.10453 

0.03229 

Jxr  measures. 

| 

PDL3 

0.57889 

0.38000 

0.27126 

0.17188 

0.07261 

0.02288 

Kgp  * 1.0 

; 
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TABLE  C-l6.  REVISED  MEASURES,  ROLL-SWAY  AXES  VERSUS  BREAK 

BREAK  FREQUENCY  AND  Cgy 


.0675  *125  .25  .5 


PDLl 

C. 06404 

0.20797 

0.55343 

0.94366 

PDL2 

0.00077 

0.00416 

0.01742 

0.03621 

PDL3 

0.03543 

0.11937 

0.35261 

0.66463 

PDL4 

0.00127 

0.00499 

0.01885 

0.04650 

PVEL 

55.76504 

56.46758 

46.98891 

41.05805 

PDIS 

85.75620 

87.46274 

82.74385 

47.78429 

PPM 

0.00017 

0.00022 

0.00036 

0.00061 

prit.2 

0. 13064 

0.10071 

0.06550 

0.03001 

PO-J 

0.9E.949 

0 . 98749 

0.97636 

0.93235 

PDI-4 

-0.91526 

-0.92514 

-0.94044 

-0.91701 

PACC 

31 .23702 

45.25141 

48.24283 

28.31316 

PVEL 

0 . 1 7E+03 

0.18E+03 

0.10E+03 

0.58E+02 

PP  IS 

0.1 BE +04 

0.94E+03 

0.40E+03 

0.17E+03 

PDLl 

0.01630 

0.01819 

0.0174S 

0.01078 

PDL3 

0.00956 

0.01481 

0.02934 

0.03886 

PDLl 

0.06404 

0.20797 

0.55343 

0.94366 

PDL2 

0.00077 

0.00416 

0.01742 

0.03621 

PDL3 

0.03543 

0.11937 

0.35261 

0.66463 

PDL4 

0.00127 

0.00499 

0.01885 

0.04650 

PVEL 

55.76504 

56.46758 

46.98891 

41.05805 

PDIS 

85.75620 

87.46274 

82.74385 

47.78429 

!.  r , 

0.</0:.3v 

0.0013V 

0.00323 

0.00550 

r 

0.01730 

0.76097 

0.54631 

0.26503 

1 

C.  79017 

0.98840 

0.97797 

0.93499 

-0.43101 

-0.49975 

-0.65545 

-0.78325 

PACC 

48.85858 

53.79541 

43.62885 

23.17445 

PVEL 

0.16E+03 

0.13E+03 

0.11E+03 

0.62E+02 

PDIS 

0.73E+03 

0.69E+03 

0.37E+03 

0.14E+03 

1.0 

2.0 

1.21167 

1.26063 

Jp  measures, 

0.03016 

-0.04607 

0.93627 

1.08140 

K^r  = 1 .0, 

0.07045 

0.01399 

26.72289 

8.58769 

Cgy  = . 0025 

18.30605 

9.33734 

0.00104 

0.00138 

Jy  measures, 

0.01231 

0.00425 

0.76758 

0.44862 

& 

w 

11 

• 

0 

-0.76099 

-0.44598 

14.79704 

10.28296 

Cgy  — . 0025 

0.25E+02 

0.10E+02 

0.52E+02 

0.16E+02 

0.00522 

0.00208 

Jy«r  measures. 

0.04014 

0.03564 

K2R  = 1 • 0 
Cgy  = . 0025 

1.21167 

1.26063 

Jp  measures. 

0.03016 

-0.04607 

0.93627 

1.08140 

•s 

0 

. 

II 

K 

0.07045 

0.01399 

26.72289 

8.58769 

Cgy  = .0225 

18.30605 

9.33734 

0.009A7 

0.01252 

Jy  measures. 

0.11113 

0.03852 

0.77185 

-0.70723 

0.45397 

-0.42770 

KgR  = 1 .0, 

14.18067 

10.37257 

CgY  = .0225 

0.25E+02 

0.98E+01 

0.48E+02 

0.16E+02 

PDLl 

PDL3 

0.13388 

0.07317 

0.15424 

0.08893 

0.15556 

0.12110 

0.09845 

0.13250 

0.04766 

0.12665 

0.01888 

0.10854 

Jyp  measures 
KgR  = 1.0, 

CgY  = • 0225 

PDLl 

PDL2 

PDL3 

PDL4 

PVEL 

PDIS 

0.06404 

0.00077 

0.03543 

0.00127 

55.76504 

85.75620 

0.20797 

0.00416 

0.11937 

0.00499 

56.46758 

87.46274 

0.55343 

0.01742 

0.35261 

0.01885 

46.98891 

82.74385 

0.94366 

0.03621 

0.66463 

0.04650 

41.05805 

47.78429 

1.21167 

0.03016 

0.93627 

0.07045 

26.72289 

18.30605 

1.26063 

-0.04607 

1.08140 

0.0 1399 
8.58769 
9.33734 

Jp  measures, 
KgR  = 1 .0, 
Cgy  = .25 

PDLl 
PDI.2 
PDI.3 
POL  ♦ 
PACC 
PVEL 
PDIS 

0.00603 

1.30986 

0.99337 

0.09124 

24.49761 

54.30840 

0.16E+03 

0.00753 

1.32757 

0.99241 

0.09040 

21.74777 

58.02134 

0.14E+03 

0.01498 
1.34221 
0.98613 
0.05361 
21.37760 
48.68298 
0. 12E+03 

0.03819 
t. 20246 
0.95597 
-0.09672 
19.47116 
27.54009 
0.7SE+02 

0.08450 

0.70853 

0.81586 

-0.31008 

15.87156 

14.45196 

0.28E+02 

0.12441 

0.28572 

0.51537 

-0.27295 

9.51661 

6.83523 

0.10E+02 

Jy  measures, 

k2R  = 1 *°> 
C2Y  = ,25 

PDLl 

PDL3 

0.58151 

0.37263 

0.63654 

0.39351 

0.72090 

0.41913 

0.67231+'  0.42521 
0.41995/  0.38547 

0.18765 

0.32607 

Jytji  measures 
K2r  = 1.0, 
Cgy  = -25 
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TABIE  C-17.  REVISED  NEASURES,  YAW  AMD  HEAVE  AXES 
VERSUS  BREAK  FREQUENCY 


0.63466  0.93944 

0.34 160  0.74246 

3.87879  4.36360 

26.99819  14.53148 


PDL1 

PDL3 

PVEL 

PDIS 


1.18190  1.14679 
0.98077  1.03396 
3.66914  2.31082 
7.93515  3.92446 


1.07502  1.03102 
1.02612  1.01335 
1.44214  0.72981 
1.84893  0.69747 


PDL1 

PDL3 

PACC 

PVEL 

PDIS 


0.83764  1.21807  1.53947  1.18993  1.06558  1.01853 

0.63075  0.85105  1.12942  1.05578  1.02024  1.00642 

0.96E+02  0. 10E+03  0.87E+02  0.56E+02  0.32E+02  0.26E+02 

0.28E+03  0.23E+03  0.19E+03  0.58E+02  0.30E+02  0.10E+02 

0.20E+04  0.14E+04  0.61E+03  0.16E+03  0.40E+02  0.94E+01 


0.08403  0.07643  0.07223  0.05933 

0.08720  0.13090  0.13625  0.14512 


0.02100  0.04448 

0.16689  0.11625 


K)Ll 


1.11171 

1.03165 


1.01189 

0.98101 


1.05259 

1.01700 


1.15905 

1.06070 


1.00686 

1.00278 


1.03069 

1.01165 


TABIE  C-l8.  REVISED  M2ASUHES,  ROLL-SWAY  AXES  VERSUS 
WASHOUT  CONFIGURATION 


Washout 


0.94366  0.94366 

0.02973  0.02973 

0.66463  0.66463 

0.03875  0.03875 

41.22489  41.22489 
48.50484  48.50484 


0.94366 
0 .245^)6 
0.66463 
0.23249 
37.05383 
3 2.50349 


PDL1 

PDL2 

PDL3 

PDL4 

PVEL 

PDIS 


0.57632  0.57632 

-0.00394  -0.00394 

0.93013  0.93013 

-0.23146  -0.23146 
28.09876  28.09876 
57.85407  57.85408 
0.16E+03  0. 16E+03 


PDL1 

PDL2 

PDL3 

PDL4 

PACC 

PVEL 

PDIS 


0.73261 

-0.00663 

47.22521 

50.11210 

0.12E+03 


rn  0.66065 
0.40566 


0.01078 

0.03886 


0.01078 

0.03886 


PDL1 

PDL3 
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APPENDIX  D 

COMPARISON  METHOD  FOR  PHASE  CHARACTERISTICS  OF 
DIFFERENT  DRIVE  LOGIC  CIRCUITS 

Parametric  study  of  drive  logic  response  characteristics  can  be  extremely 
complex  if  many  parameters  are  involved.  It  is  also  difficult  to  compare 
the  characteristics  of  drive  logic  for  different  simulators  which  may  not 
only  have  different  parameter  definitions,  but  may  also  have  different  drive 
logic  forms.  Both  of  these  circumstances  motivate  development  of  simple 
metrics  which  can  serve  to  quantify  basic  dynamic  performance  of  drive 
logic  circuits  independently  of  their  complexity  or  form.  The  metrics 
proposed  here  focus  upon  the  effective  phase  lead  introduced  by  washout 
circuits  at  high  frequencies,  and  the  effective  phase  lag  introduced  by 
residual  tilt  circuits  at  low  frequencies.  The  matter  of  drive  logic 
scaling  gain  constitutes  another  key  characteristic  in  itself.  Its  sig- 
nificance is  obvious,  and  it  will  not  be  considered  further  here. 

The  method  for  calculating  effective  leading  and  lagging  phase  charac- 
teristics is  based  upon  Pade  approximants  for  exponential  functions.  A table 
of  these  approximants  is  given  below. 

Use  of  these  approximations  is  illustrated  with  some  simple  examples. 
EXAMPLE  #1 

Compare  the  drive  logic  for  the  roll  axes  of  two  simulators. 

Simulator  #1* 

Roll  washout  transfer  function  is: 

s s 

s + cqp  s + cap 


♦The  Vought  Corp.  LAMBS  simulator  roll  axis  drive  logic  has  this  form. 
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TABLE  D-l.  APPROXIMATIONS  FOR  DRIVE  LOGIC 
TRANSFER  FUNCTION  FACTORS 


/ 


For  I s I < i/t  or  l/v  (l.e.,  at  low  frequencies) 


1 + TS 


2 2 

1 + tjs  + v s 


= e 


■TS 


-TjS 


, , . TS 

1 + ts  = e 


. , x 2 2 . tis 

1 + ns  + v s = e ' 


For  Isl  > a or  y (i.e..  At  high  frequencies) 


- 0 + V-) 


+ a/s)  1 = e'01/8 


= ( 1 + P/s  + 72/s2)  1 « e‘p/s 


2 2 
s + ps  + 7 


s + a 


= 1 


*a/s  i e^8 


* ? * ^ = !♦<>/>♦  r2/.2  > *s/s 


19U 


vvfc- 


• . Jo.r^T  ..tAfl 


'*"'**  >:  \y'\ . y } '■  ' * ~ 


The  effective  lead  at  high  frequencies  is  obtained  from 

s g a e-a*p/g  ."Ofc/8  . + «&)/■ 

8 + ffi(p  8 + <Up 


using  the  fifth  entry  in  Table  D-1 . The  effective  lead 
frequency  parameter  is  (<ity  + a>p)  rad/sec. 

Simulator  #2* 

Roll  washout  transfer  function  is: 

s5 

(s  + 0.65)[s2  + 2(0.T07)0.65s  + (0.65)2] 

. e-0.65/s  e-2(0. 707)0.65/8  = e-1 .57/s 

The  effective  lead  frequency  parameter  is  1-57  rad/sec. 
INTERPRETATION 

In  order  that  the  roll  washout  for  Simulator  #1  introduce  less 
phase  (lead)  distortion  than  Simulator  #2  it  is  necessary  that 
((Up  + (Up)  < 1-57  rad/sec.  In  order  that  their  phase  distortions 
be  comparable,  (cup  + (Up)  = 1.57  rad/sec.  Furthermore,  notice  that 
alternative  selections  for  the  city  and  oip  parameters  of  Simulator  #1 
satisfying 

a*p  + a^>  = 1.57  rad/sec 

are  possible.  Still  further,  all  alternatives  satisfying  the  above 
equation  are  equivalent  in  terms  of  the  phase  distortion  metric 
used  here. 

The  usefulness  of  this  approach  for  reducing  the  number  of  parameters 
which  must  be  varied  in  a parametric  study  is  illustrated  by  the  last  point 
above.  By  way  of  further  illustration,  consider  equivalent  break  frequency 
combinations  for  the  first  and  second  order  transfer  function  factors  for 
Simulator  #2.  Any  combination  of  break  frequencies  satisfying 

ojfo  + 2(0.707)ooso  = 1*57 


♦The  Air  Force  LAMARS  simulator  roll  axis  drive  logic  for  the  F-17  is 
reported  to  have  this  transfer  function  (Ref.  l). 
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will  have  equivalent  phase  distortion  metrics.  Further,  if  the  factor 
(s  + 0.3)(s  + 0.l)/(s  + 0.2)2  is  included  in  the  transfer  function  for 
Simulator  #2,  the  phase  distortion  metric  will  be  unchanged. 

EXAMPLE  #2 

Evaluate  the  high  and  low  frequency  phase  distortion  in  lateral,  specific 
force  reproduction  for  the  nominal  configuration  of  the  Vought  Corp.  LAMBS 
drive  logic  with  and  without  addition  of  a lateral,  residual  tilt  feature. 

NOMINAL  LAMBS 

Lateral  specific  force  transfer  function  is: 

2 2 

0.20  -r 5 -5 

s + 2(0.707)2.0  s + (2.0)  s + 0.45s  + 0.1 


RESIDUAL  TILT  ADDITION 

Residual  tilt  transfer  function  is: 


Phase  lead  distortion  at  high  frequencies  without  residual  tilt  is 
obtained  from  the  first  transfer  function.  The  effective  lead  fre- 
quency parameter  is: 

2(0.707)2.0  +0.45  = 3.26  rad/sec 

The  phase  lag  distortion  at  low  frequencies  is  asymptotic  to  —360  deg 
in  the  absence  of  residual  tilt.  However,  the  amplitude  ratio  at  low 
frequencies  is  very  small  (proportional  to  frequency  to  the  fourth 
power). 

Phase  distortion  with  residual  tilt  is  obtained  by  adding  the  two 
transfer  functions  and  factoring  the  resulting  numerator  polynomials. 


fj 

i 


! 


i 
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Phase  lead  distortion  at  high  frequencies  with  residual  tilt  results 
in  the  effective  lead  frequency  parameter: 

-l.4i3  - 2(0.713)0.32  - 2(-0.03l)l.89  + 0.5  + 2(0.707)2.0  + 0.45 

= 2.03  rad/sec 

Phase  lag  distortion  at  low  frequencies  with  residual  tilt  results  in 
the  effective  delay  parameter: 


-1  - ilQi.p.3)  , giz9-03i ) + _L  + 2(0.707)  + 
1.413  0.32  I.89  0.5  2.0 


= 2.08  sec 

The  effective  phase  lag  is  asymptotic  to  0 degrees. 

INTERPRETATION 

Addition  of  residual  tilt  reduces  leading  phase  distortion  at  high 

frequencies  significantly. 

Addition  of  residual  tilt  reduces  lagging  phase  distortion  at  low 

frequencies  dramatically. 

As  a general  rule,  small  values  of  the  effective  lead  frequency  para- 
meter, and  of  the  effective  delay  parameter  are  desirable  from  the  point  of 
view  of  minimizing  phase  distortion  in  simulated  motion  cues.  Values  between 
1.0  and  2.0  rad/sec  for  effective  lead  frequency  and  between  0.5  and  1.0  sec 
for  effective  delay  are  to  be  considered  "good".  Certain  cautions  must  be 
presented  dong  with  this  information,  however. 

• "Good"  values  of  the  effective  lead  frequency  or  effective 
delay  parameters  obtained  by  means  of  extreme  motion  scaling 
in  the  drive  logic  are  counter-productive. 

• Given  motion  scaling  and  simulation  tasks  will  place  lower 
bounds  on  the  effective  lead  frequency  parameters  based  on 
the  available  simulator  travels. 


f 
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Spurious  motion  cue  generation  will  result  in  bounds  on  the 
effective  lead  frequency  and  effective  delay  parameters.  A 
lover  bound  on  the  effective  delay  parameters  is  set  by  the 
tolerable  level  of  spurious  pitch  and  roll  rate  cues.  These 
cues  result  from  the  residual  tilt  feature.  A lower  bound 
on  the  effective  lead  frequency  parameters  in  the  pitch  and 
roll  axis  drive  logic  and  an  upper  bound  on  the  lead  fre- 
quency parameters  in  the  longitudinal  and  lateral  transla- 
tional drive  logic  is  set  by  the  tolerable  level  of  spurious 
longitudinal  and  lateral  specific  force  cues.  These  cues 


are  the  result  of  imperfect  simulator  coordination. 


